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1.0  INTRODUCTION 


This  manuscript  represents  the  final  report  for 
Geophysics  Laboratory  Contract  No.  F19628-86-K-0036; 
Investigation  of  Structure  and  Transport  Processes  in  the  High 
Latitude  Ionosphere  Plasma  Using  Aircraft  and  Ground  Based  Digital 
lonosondes. 

The  project  addressed  three  topics  identified  by  a  GL 
Program  Research  and  Development  Announcement: 

i)  Specification  of  structure,  dynamic  behavior  and  transport  of 
ionization  at  high  latitudes. 

ii)  Measurements  and  long  term  data  base  development. 

iii)  Analysis  and  interpretation. 

The  University  of  Lowell  Center  for  Atmospheric 
Research  (ULCAR)  has  designed  and  built  an  advanced  ionospheric 
sounder,  the  Digisonde.  Several  of  these  have  been  purchased  and 
deployed  by  the  Air  Weather  Service  (AWS)  for  real  time  ionospheric 
data  acquisition.  The  Geophysics  Laboratory  (GL)  uses  several  high 
latitude  systems  for  scientific  studies.  ULCAR  contracted  to  use  these 
instruments  to  address  three  specific  aspects  of  the  above  topics: 

i)  Determination  of  ionospheric  convection  patterns. 

ii)  Observations  of  the  E-region  trough. 

iii)  Definition  of  the  source  and  structure  of  polar  cap  arcs  and  F- 
layer  patches. 

For  the  purposes  of  this  report,  the  research  has  been 
broken  down  into  three  different  categories:  analysis  techniques 
(Section  2.0);  observations  and  campaigns  (Section  3.0);  and 
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interpretation  of  data  (Section  4.0).  Section  5.0  lists  the  publications 
and  presentations  that  resulted  from  this  research. 

The  high  latitude  ionosphere  is  a  complex  region  because 
of  the  two  forms  of  solar  energy  transfer:  electromagnetic  radiation 
and  solar  wind  particles.  The  coupling  of  the  magnetic  fields  of  the 
earth  and  the  sun  (or  the  solar  wind)  leads  to  magnetospheric 
electric  fields  that  control  the  dynamics  in  the  magnetosphere.  These 
electric  fields  arc  also  mapped  along  magnetic  field  lines  into  the 
polar  cap  ionosphere  producing  an  E  x  B  drift  that  drives  the  plasma 
across  the  polar  cap  [Dungey,  1961;  Cauffman  and  Gurnett,  1972; 
Heppner,  1977].  Magnetospheric  accelerated  particles  impacting  onto 
the  atmosphere  produce  the  aurora  and  the  auroral  oval  surrounding 
the  polar  cap.  The  size,  or  diameter,  of  the  polar  cap  and  the  cross 
polar  cap  potential  are  direct  measures  of  the  energy  transfer  from 
the  solar  wind  into  the  ionosphere.  F  region  plasma  enhancements, 
or  patches,  observed  to  move  in  the  antisunwara  direction  across  the 
polar  cap  [Buchau  et  al.,  1983;  Weber  et  al.,  1984]  have  been 
explained  as  the  result  of  variations  in  the  diameter  of  the  polar  cap 
[Anderson  et  al.,  1988].  When  the  polar  cap  convection  pattern 
expands  so  that  sunlit  regions  at  noon  magnetic  local  time  (MLT) 
become  part  of  the  pattern,  then  higher  density  plasma  is  carried 
from  the  sunlit  ionosphere  into  the  dark  polar  cap,  with  its  lower 
electron  density. 

The  study  of  the  convection  pattern  is  one  of  the  major 
tasks  of  GL's  Digisonde  network  (Figure  1).  In  Figure  2,  the  locations 
of  the  four  high  latitude  Digisondes  are  shown  in  a  corrected 
geomagnetic  latitude  versus  geomagnetic  local  time  (CGLT)  system 
[Hakura,  1965;  Whalen,  1970].  The  locations  of  the  Millstone  Hill  and 
Else  AT  incoherent  scatter  radars  are  also  shown  for  reference.  The 
initial  observations  were  conducted  at  Qaanaaq,  Greenland  (87° 
CGLAT),  Goose  Bay,  Labrador  (65°  CGLAT),  and  Argcntia, 
Newfoundland  (57°  CGLAT).  The  station  locations  are  marked  by  Q, 
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G,  and  A  for  02  universal  time  (UT),  when  Qaanaaq  is  at  midnight 
magnetic  local  time.  The  plasma  convection  pattern  is  monitored  at 
its  center  by  Qaanaaq,  and  in  the  auroral  oval  and  trough  regions  by 
the  Goose  Bay  and  Argentia  stations,  respectively.  A  fourth  high 
latitude  Digisonde  was  installed  at  Sondre  Stromfjord,  Greenland  (77° 
CGLAT)  in  June  1989,  and  the  first  observations  are  discussed  in  this 
report.  The  Sondrestrom  Digisonde  fills  the  gap  between  Qaanaaq 
and  Goose  Bay.  The  Sondrestrom  station  is  important  because  it 
monitors  the  plasma  flow  in  the  dayside  throat  region  and  enables 
the  development  of  ionospheric  irregularities  to  be  investigated  as 
they  cross  the  polar  cap. 

For  the  modeling  and  forecasting  of  the  high  latitude 

ionosphere  the  plasma  convection  is  an  important  driver,  and  one  of 
the  exciting  results  of  the  work  reported  here  is  the  validation  of  the 
Digisonde  drift  measurements  as  a  monitor  of  the  polar  cap 

convection.  The  characteristic  flow  patterns  as  a  function  of  By,  for 

Bz  negative,  had  been  established  by  statistical  evaluation  of  satellite 
plasma  flow  measurements  [Heppner  and  Maynard,  1987],  Using 
similar  statistical  techniques,  we  have  evaluated  the  Digisonde  drift 
data  from  the  central  polar  cap  station  at  Qaanaaq,  and  found 

excellent  agreement  with  the  Heppner-Maynard  models. 

Understanding  of  the  polar  ionization  distribution 
involves  the  entire  polar  cap/oval/suboval  region,  including  the 
formation  of  the  night  and  daytime  F  region  trough  [Whalen,  1989, 
and  references  therein],  Whalen  concluded  from  his  analysis  of 
ionosonde  data  that  the  trough  in  the  daytime  F  layer  is  a  regular 
feature,  controlled  by  the  magnetospheric-ionospheric  convection. 
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2.0  ANALYSIS  AND  MEASURING  TECHNIQUES 


2.1  Electron  Densities 

The  labor  intensive  ionogram  scaling  and  evaluation 
required  by  conventional  ionosondes  has  been  replaced  by  automatic 
techniques.  The  Automatic  Real  Time  Ionospheric  Scaling  Technique 
(ARTIST)  developed  by  ULCAR  is  a  PC-based  system  which 
determines  values  for  the  important  parameters  from  an  ionogram. 
The  operation  and  testing  of  ARTIST  has  been  described  by  Gamache 
et  al.  [1988]  associated  with  this  contract,  and  by  Reinisch  et  al., 
1988a. 


Inversion  of  ionograms  to  yield  electron  density  profiles 
is  also  performed  automatically  by  ARTIST.  The  current  software 
includes  two  inversion  techniques.  These  are  the  ULCAR  modified 
Chebyshev  polynomial  method  and  the  Titheridge  POLAN  method 
[Titheridge,  1985].  Both  are  described  in  Scientific  Report  No.  2,  and 
by  Reinisch  et  al.  [1988b].  A  third  representation  of  the  F-region 
electron  density  profile  in  terms  of  LAY  functions  has  also  been 
investigated  [Bossy  et  al.,  1988].  McNamara  et  al.  [1987]  compiled 
statistics  on  the  comparison  of  hmF2  values  derived  from  the  POLAN 
model  with  those  obtained  by  the  simple  methods  of  Dudeney  [1983] 
and  Bilitza  et  al.  [1979]. 


2.2  Plasma  Drift 


A  very  important  aspect  of  the  Digisonde  is  its  ability  to 
measure  ionospheric  plasma  drift  velocities,  from  which  electric 
fields  may  be  derived.  Until  early  1986  drift  measurements  had 
been  made  by  manually  selecting  frequencies  and  setting  sampling 
gates,  thereby  requiring  continuous  manning  of  the  Digisonde  during 
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drift  measurements.  During  the  time  covered  by  this  contract,  the 
drift  measurements  have  been  automated.  The  drift  analysis 
programs  (running  on  GL's  Cyber  computer)  ’^ave  also  been 
significantly  improved.  This  work  is  summarized  .  w. 


2.2.1  Automation  of  Drift  Measurements 

The  automatic  drift  mode  operations  of  the  Digisonde  are 
based  on  ARTIST  trace-identification  software.  The  drift 
measurements  use  phase  coherent  pulse  transmission  on  two  HF 
frequencies  which  are  selected  to  sample  desired  heights  in  the  F- 
region.  During  normal  measurement  campaigns,  ionograms  are  taken 
in  5 -minute  intervals  to  monitor  the  dynamic  environment.  The 
time  between  ionograms  is  used  to  collect  a  series  (nominally  twelve) 
of  drift  measurements  (cases)  of  typically  10-s  each. 

The  commonly  used  drift  mode  in  the  Digisonde  256 
multiplexes  four  sounding  frequencies  and  samples  the  echoes  at  one 
range  gate  for  each  of  the  four  frequencies.  Automation  of  the  drift 
operation  requires  selection  of  appropriate  frequencies  FI,  F2,  F3  and 
F4,  corresponding  range  gates  HI,  H2,  H3  and  H4,  and  receiver  gain 
settings  Gl,  G2,  G3  and  G4.  The  "AUTODRIFT"  software  residing  in 
the  ARTIST  PC/AT  selects  frequencies  and  range  settings  from  the 
autoscaled  ionograms.  To  allow  for  height  (or  range)  changes  during 
the  drift  observation  intervals,  a  range  tracking  algorithm  is  included 
in  AUTODRIFT.  The  initial  settings  are  FI  =  F2  and  F3  =  F4;  and  H2  = 
HI  +  10  km  and  H4  =  H3  +  10  km.  At  the  end  of  each  drift  case 
(nominally  10  sec)  the  amplitudes  at  HL  (larger  range)  and  HS 
(smaller  range),  i.e.  AL  and  AS  are  compared  and  the  range  gates  are 
corrected  for  the  next  drift  case  in  the  following  way: 

If  AS  >  AL  ->  decrease  KL  by  20  km 

If  AS  <  AL  ->  increase  HS  by  20  km. 
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The  two  sounding  frequencies  are  selected  at  the  end  of  the 
ionogram,  as  soon  as  the  autoscaling  is  completed,  and  then  remain 
constant  during  the  drift  interval.  The  frequency  selection  criteria 
are  simple:  one,  find  interference-free  frequencies,  and  two,  choose 
frequencies  that  are  reflected  at  the  bottom  of  the  F  layer.  The 
reason  for  the  first  criterion  is  clear;  good  signal-to-noise  ratio.  The 
rationale  for  the  second  criterion  is  the  Digisonde  drift  analysis 
technique  [Reinisch  et  al.,  1987]  which  is  summarized  in  the 
following  paragraphs. 

A  wide  beam  transmit  antenna,  fed  by  10  kW  of  pulsed 
HF  signal,  illuminates  a  large  area  of  the  F  region  above  the 
ionosonde.  All  areas  in  the  F  region  where  the  plasma  frequency  is 
equal  to  the  sounding  frequency  and  where  the  area  surface  vector  is 
parallel  to  the  incident  wave  vector  generate  an  echo.  The  reflection 
point  is  called  a  source  [Dozois,  1983].  The  Doppler  shift  of  the  source 
signal  is  determined  by  the  location  and  the  velocity  of  the  source. 
The  Digisonde  determines  the  Doppler  shifts  by  forming  the  discrete 
Fourier  transforms  [Bibl  and  Reinisch,  1978]  of  the  signals  received 
on  four  or  seven  spaced  receiving  antennas.  To  find  the  locations  of 
the  sources,  two  assumptions  are  made:  one,  that  one  source 

contributes  to  a  given  Doppler  line,  and  two,  that  the  incidence  angle 
calculated  by  using  the  measured  spectral  phases  at  the  spaced 

antennas  specifies  the  source  location.  This  latter  assumption  is 
justified  when  refraction  in  the  ionosphere  is  small,  which  is  the  case 

for  signals  reflected  at  the  bottom  of  the  F  layer  where  the  vertical 

plasma  density  gradient  is  steep.  It  is  for  this  reason  that 
AUTODRIFT  selects  frequencies  in  the  flat  part  of  the  F  trace  in  the 
ionogram.  The  discussion  on  refraction  effects  and  on  the 
source/Doppler  relation  is  continued  in  Section  2.2.2. 

The  method  of  selecting  the  drift  sounding  frequencies  is 
as  follows.  The  AUTODRIFT  algorithm  divides  the  frequency  interval 
foF2  -  foFl  (or  foF2  -  foE,  or  foF2  -  fminF,  depending  on  ionospheric 
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conditions),  obtained  from  the  last  ionogram,  into  four  equal 
intervals.  From  the  two  central  intervals,  one  frequency  with  the 
best  signal-to-noise  ratio  (as  measured  in  the  ionogram)  is 
determined  for  each  interval.  The  virtual  echo  heights  for  these  two 
frequencies  are  obtained  from  the  ionogram  autoscaling  data.  The 
initial  receiver  gain  control  is  set  to  -20  dB  (G  =  2).  For  each  of  the 
two  frequencies,  the  gain  is  adjusted  at  the  end  of  each  drift  case 
until  the  larger  signal  amplitude  (of  the  two  range  gates)  is  within  6 
to  12  dB  of  the  maximum  allowable  digital  signal. 

The  first  automatic  drift  measurements  were  made  at 
Argentia,  Newfoundland,  in  October  1986.  Since  then,  several 
adjustments  have  been  made  to  the  AUTODRIFT  program,  which  is 
now  operating  reliably  at  Argentia,  Goose  Bay,  Sondrestrom,  and 
Qaanaaq. 


To  measure  in  the  AUTODRIFT  mode,  Digisonde  preface 
parameter  P3  [Bibl  et  al.,  1981]  must  be  >  8.  Different  automatic 
operational  modes  can  be  selected  by  the  appropriate  setting  of  the 
8,  4,  2  and  1  bits  of  preface  parameter  P7: 


Bit 

On 

Off 

8 

automatic  frequency  selec¬ 
tion  from  ionogram 

manual  frequency  selection 

4 

automatic  gain  control 

manual  gain  control 

2 

automatic  range  tracking 

manual  range  setting 

1 

select  sounding  frequen¬ 
cies  close  to  foF2  from 
ionoeram 

select  center  frequencies 
from  ionogram 
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For  routine  drift  observations  the  1-bit  is  off  and  the  two  sounding 
frequencies  are  selected  close  to  the  center  of  the  F2  trace  as 
discussed  before.  When  the  1-bit  is  on,  the  sounding  frequencies  are 
selected  from  the  interval  0.7  foF2  <  f  <  0.9  foF2. 


2.2.2  Echo  Source  Determination 

This  section  discusses  the  assumption  mentioned  in  the 
previous  section,  that  the  spectral  (Doppler)  analysis  uniquely 
identifies  the  echo  sources,  i.e.  that  each  Doppler  line  corresponds  to 
one  source.  This  is,  of  course,  not  generally  true  since  even  for  a 
horizontal  reflection  plane  that  moves  with  velocity 

V  =  (vx,  Vy,  Vz)  (1) 

there  exists  a  line  of  potential  source  points  that  have  identical 
Doppler  frequency 


d 


TC 


—  r»v 


A,r 


(2) 


where  X  is  the  wave  length  and  r  the  position  vector  of  the  source. 
The  loci  for  which  the  Doppler  frequency  or  the  radical  velocity  Vr  is 
constant  is  given  by 


r 

— .  v  =  Vf  =  const, 
r 

The  position  vector  r  is 
r  =  (x,  y,  H), 


(3) 
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where  H  is  the  height  of  the  reflecting  plane.  By  appropriate  choice 
of  the  X,  y-directions,  equation  (1)  becomes 


V  =  (vx,  0,  Vz), 

and  equation  (3)  can  be  written  as 
(xvx  +  zvz)2  =  r2  vr^, 
or  for  Vx  ^  0: 

(x  +  xq)^  -  Ay2  =  AH2 


(4) 

(5) 

(6) 


where 


A  = 


I 

(Vr/  Vx)  _ 

I 


This  is  a  set  of  hyper  das  shifted  along  the  x-axis  by 


2 

l-(Vr/Vx)  (7) 

Figure  3  shows  the  lines  of  equal  Doppler  velocity  for  different 
values  of  vr/vx  for  a  reflection  height  of  H  =  250  km,  and  a  vertical 
velocity  of  vz  =0.1  vx.  The  zero  Doppler  line  is  the  straight  line 
X  =  -  XQ.  When  the  vertical  velocity  vz  is  zero,  the  zero  Doppler  line 
goes  through  zero  (zenith). 

The  probability  that  iso-Doppler  lines  are  populated  with 
several  sources  is  small,  since  the  vector  normal  to  the  reflecting 
surface  must  be  parallel  to  the  wave  vector  of  the  incident  wave. 
However,  when  it  happens  that  two  sources  with  the  same  Doppler 
exist,  like  sources  SI  and  S2  for  vr/vx  =  0.5  in  Figure  3,  then  an  error 


1 1 


Lines  of  «<^al  Doppler  velocity.  Reflecting 
•lone  is  at  250  km  altitude  moving  with  velocity 
'**(Vx,0,V2).  It  is  assumed  thot  vi  is  10»  of  w. 


Figure  3.  Lines  of  Equal  Doppler  Velocity 
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in  the  calculation  of  the  source  location  may  occur.  Dozois  [1983]  has 
shown,  that  dependent  on  the  amplitudes  and  initial  phases  of  the 
two  sources,  the  apparent  Doppler  source  will  be  put  somewhere 
between  SI  and  S2,  like  point  A  in  Figure  3.  Fortunately,  the 
distance  of  A  from  its  Doppler  line  is  small  for  most  configurations, 
causing  only  small  errors  in  the  drift  velocity  calculation. 

Only  parts  of  the  iso-Doppler  lines  are  within  the 
Digisonde's  field  of  view  which  is  limited  by  the  fact  that  the 

ionosphere  is  illuminated  with  HF  pulses  of  66  |is  duration, 

corresponding  to  10  km  in  range.  For  a  perfect  plane  at  H  -  250  k.*! 
this  limits  the  field  of  view  to  a  circle  with  71  km  radius,  i.e.  a  zenith 
angle  of  16°.  Where  the  equidensity  surfaces  are  rough  with  height 

variations  of  Az  =  40  km,  the  radius  becomes  166  km,  corresponding 

to  a  zenith  angle  of  34°.  The  circles  with  R  =  75  km  and  R  =  175  km 
in  Figure  3  illustrate  the  fields  of  view  for  different  ionospheric 
roughness. 


2.2.3  The  Effects  of  Refraction 

To  estimate  the  effects  of  ionospheric  refraction  on  the 
drift  analysis,  a  simple  model  calculation  was  conducted,  that 

neglects  the  earth's  magnetic  field  and  horizontal  gradients  in 

ionization.  The  electron  density  N(h)  in  the  horizontally  stratified 
ionosphere  increases  linearly  with  height 

N(h)  =  a  (h-ho)  for  h  >  h© 

=  0  for  h  <  ho  (8) 

Reflections  occur  at  overdense  irregularities  on  a  point  A  on  the  ray 
path  such  that  u  =  0  at  this  point,  but  ixa  0  just  before  (see  Figure 
4).  In  our  model,  the  irregularities  move  together  with  the 
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ionosphere  with  the  velocity  v.  The  rate  of  change  of  the  phase 
P,  i.e.  the  Doppler  velocity  vr  is  [Kelso,  1964,  p.  278] 


v,  =  ^  =  2 


•[i(h,t)dS 


m  =  Unit  vector  in  the  ray  direction  at  point  v. 
v  =  Velocity  at  point  A. 

=  Index  of  refraction  at  point  A. 

P  =:  Phase  path. 

Integral  is  over  the  ray  path  S. 

Vertical  and  horizontal  velocities  are  considered  separately.  Fi 
assume  the  ionized  layer  is  drifting  upward  at  a  velocity  Vz.  T 
the  following  relationships  hold: 

N(h,t)  =  a(h-ho-v,t)  f 


H%,t)  =  1  -  ^(h-  ho  -  v^t) 


j  Ito  forh>ho  +  v,t 

^^l(h,t)  =  /  ^2  2^l 

1  0  for  h  <h„  +  v.,t 
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ionosphere  with  the  velocity  v.  The  rate  of  change  of  the  phase  path 
P,  i.e.  the  Doppler  velocity  vr  is  [Kelso,  1964,  p.  278] 


m  =  Unit  vector  in  the  ray  direction  at  point  v. 

V  =  Velocity  at  point  A. 

|1A.  =  Index  of  refraction  at  point  A. 

P  =  Phase  path. 

Integral  is  over  the  ray  path  S. 

Vertical  and  horizontal  velocities  are  considered  separately.  First, 
assume  the  ionized  layer  is  drifting  upward  at  a  velocity  vz.  Then 
the  following  relationships  hold: 


N(h,t)  =  a(h  -  ho-  v^t) 

n%,t)  =  l  -^(h-h„-v,t) 
1 


(10) 


(11) 


f  ll(h,t)  = 

at 


Vz 
2  11 


for  h  >ho  +  Vj,t 
for  h  <ho  +  V2t 


(12) 


dS  = 


pdh 


2  2 
[L  -  sin  0 


0 


(13) 
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(14) 


d  h  =  —  2ndji 
Ka 


Using  (10)  -  (14)  in  the  integral  portion  of  equation  (9): 


—  p(h,t  dS=  — dh_^^ 
at  2  2  *^2^1)  j— 1  ^  2^ 

f  V  p  -  sin  00 


=  -  V 


>-Pa 

'^=1  ri  2^  ' 

V  p  -  sin  00 


Let  u2  =  p2  .  sin20o  and  the  integral  becomes: 


- - - — .-iR/ 

j  du  =  -  VzU=-Vz  V  p  -  sin  00^ 


a  ./  2  /  2  2 

—  p(h,t)dS  =  V  1  -  sin  0o  -  v  p^  -  sin 
0t 


Snell's  law  relates  the  angle  0A  to  the  index  of  refraction 
pA  and  the  zenith  angle  0o: 


PAsin0A  =  sin0o  or  Pa=-^ — ! 

sin0 
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— ^i(h,t)dS  =  V,  cose„ - 


cose, 


Add  to  this  the  dot  product: 


sine. 

p^m»v  = - 2.  V  2  cose.. 

sine.  ^ 


to  get  a  Doppler  velocity  of:  ^  =  2  v ,  cose 

dt  ^  0*  (21) 

This  means  that  for  this  case  the  true  velocity  is  correctly  obtained 
from  the  measured  radial  velocity: 

V  =  vz  =  vr/2coseo.  (22) 

Assume  now  that  the  ionized  layer  is  moving  in  the  +x 
direction  with  a  velocity  vx.  In  this  case  N  and  4  are  no  longer 

functions  of  time  (between  the  bottom  of  the  layer  and  the  position 
of  the  irregularity)  so: 

N(h)  =  a  (h  -  hg) 


Mh)  =  l  ~(h-  h„) 


■li(h)  =  0 
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(26) 


—  H(h)dS  =0 

at 


The  irregularity  moves  with  velocity  vx,  so: 


^Am*v  =^aVxCos(90-  0^)  =  Ha  sine^ 


=  V. 


sin9 

sin0 


-  sinSA  =  Vx  sin0Q. 


(27) 


By  the  principle  of  superposition,  if  the  ionospheric 


velocity  is 


V  =  Vx  ax  +  vz  az 
the  Doppler  velocity  becomes: 

^  =  2(v2Cos0o  +  VxSin0j 


(28) 


(29) 


or  (see  Figure  4) 


dP 

dt 


=  2  V  ‘ar. 


(30) 


Bullett  [1990]  has  generalized  this  result  to  nonlinear  profiles  and 
has  shown  that  the  measured  rate  of  change  of  the  phase  path  is  the 
dot  product  of  the  actual  velocity  and  the  unit  vector  ar  in  the 
direction  of  the  ray  at  the  ground.  These  calculations  show  that  at 
least  to  first  order  there  are  no  refraction  effects.  Inclusion  of  the 
magnetic  field  will  not  fundamentally  alter  these  results. 
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2.2.4  Drift  Data  Processing 

The  ability  of  the  Digisonde  to  make  drift  measurements 
automatically  is  an  exciting  achievement,  which  led  to  a  new 
difficulty:  how  to  handle  the  vast  quantity  of  data  generated  by 

these  measurements.  A  significant  amount  of  effort  has  been 
involved  in  mitigating  this  problem,  leading  to  the  development  of 
more  efficient  techniques  which  take  advantage  of  the  strengths  of 
the  GL  Cyber  computer  on  which  the  processing  is  currently  done. 

The  plasma  drifts  were  originally  computed  by  the 
Frequency-Wavenumber  Power  Density  (Scanning)  technique 
described  by  Dozois  [1983].  The  FWPD  was  calculated  for  a  41  x  41 
pixel  area  covering  the  Digisonde  viewing  area.  Only  certain  spectral 
lines  were  used,  selected  according  to  their  amplitude.  The  FWPD 
method  used  both  amplitude  and  phase  measurements  to  determine 
source  locations.  Originally,  the  FWPD  method  was  thought  to  be  able 
to  separate  multiple  sources  contributing  to  the  same  spectral  line. 
However,  Dozois  [1983]  determined  that  multiple  sources  could  not 
be  reliably  separated  by  this  method,  which  required  10,000  CPU 
seconds  for  the  processing  of  one  magnetic  tape  containing  one  day 
of  drift  measurements. 


Dozois  simulated  the  multiple  sources  case  and  found  that 
two  sources  can  be  separately  identified  when  their  angular 
separation  is  larger  than  the  6  dB  beamwidth  of  the  antenna  array, 
otherwise  the  FWPD  method  determines  one  apparent  source  located 
somewhere  between  the  two  sources.  The  exact  location  depends  on 
the  relative  amplitudes  and  phases  of  the  two  sources,  which  are 
unknown.  As  discussed  in  Section  2.2.3,  the  apparent  source  location 
is  generally  close  to  its  proper  iso-Doppler  line,  so  that  the  effect  on 
the  accuracy  of  the  drift  vector  calculation  is  small.  Nevertheless, 
being  unable  to  reliably  separate  multiple  sources  the  computer 


intensive  FWPD  method  was  discarded  in  favor  of  a  computationally 
more  economic  technique  that  calculates  the  source  location  using  the 
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phase  only  (PO),  assuming  that  there  is  only  one  source  (for  each 
Doppler). 


In  the  PO  method,  spectral  amplitudes  are  used  as 
criteria  for  data  selection/rejection,  while  source  locations  are 
computed  from  phases  only.  The  approach  used  is  as  follows.  For  a 
given  source  k  with  Doppler  dk,  the  phase  <j)ik  measured  at  antenna  i 
is 

<t>ik  =  k*ri  +  8  +  27mi,  (31) 

where 

k  =  iTtlX 

X,  =  wavelength 

5  =  initial  phase  of  the  source 

271  ni  =  phase  ambiguity. 

ri  is  the  horizontal  antenna-position  vector  pointing  from  the  central 
antenna  1  to  antenna  i.  The  seven  antenna  array  used  for  these 
measurements  is  shown  in  Figure  5.  (Antennas  1,  5,  6  and  7  are 
used  for  four-antenna  measurements.)  Only  source  locations  that 
produce  phase  differences  of  less  than  n  at  the  spaced  antennas  are 
considered  here.  These  sources  are  located  in  the  primary  region  for 
which  the  zenith  angle  is  smaller  than  Zmax»  where 

sin(ZmJ  =  -  (32) 

a 

where  a  is  the  antenna  separation  between  antennas  2-3,  3-4  or  2-4. 
For  these  sources  the  integer  ni  is  0  or  1:  0  if  k*ri  is  positive;  1  if 
k»ri  is  negative,  i.e.  the  (positive)  measured  phase  is  greater  than 
k*ri  by  2k.  Sources  outside  Zmax  are  "aliased"  into  the  primary 
zone,  by  the  sidelobes  of  the  array.  The  same  situation  existed  for 
the  FWPD  calculations. 
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(Magnetic) 


(Corrected 

Geomagnetic) 


The  relative  directions  of  geographic  north,  corrected  geomagnetic 
north  and  compass  north  (i.e.  as  measured  with  a  compass)  are 
indicated  for  key  Digisonde  locations.  Compass  north  was  taken  from 
navigation  charts  for  1982  in  each  case  except  Goose  Bay,  which  was 
for  1983.  The  antenna  arrays  are  aligned  with  compass  north  (#6  -> 
#5). 


Figure  5.  Seven  Antenna  Receiving  Array 


The  symmetry  of  the  four-  and  seven-antenna  arrays 
used  for  drift  measurements  allows  us  to  determine  which  values 
need  to  be  decremented  to  obtain  the  corrected  phase  k»ri  -h  5.  A 
least-squares  fit  procedure  is  then  used  to  determine  the 
components  of  the  vector  k.  The  least-squares  error  is 


,  N-l  N 

;=x  X 

i=l  j=i+l 


(k.r,-k.rj)-(,t,;k-4.iJ  /  [n(N-iV2] 


(33) 


where  N  is  the  number  of  antennas  (typically  7  or  4),  and  (j)'ik  and 
(j)'jk  are  the  corrected  phases;  by  working  with  phase  differences,  the 
initial  phase  5  cancels  out.  (The  denominator  in  the  above  equation 
is  the  total  number  of  terms  in  the  double  sum,  normalizing  ek  with 
regard  to  the  number  of  antennas  used.)  Taking  the  derivatives  of 
w.r.t.  kx  and  ky  and  setting  them  each  equal  to  zero  yields  the  two 
equations  needed  for  determining  the  unknowns  kx  and  ky.  Sources 
with  large  r.m.s.  errors  are  deleted  from  the  analysis. 

A  new  skymap  printing  program  MAPTILT  was 
developed  that  permits  different  types  of  skymaps  to  be  displayed. 
One  of  the  new  features  is  the  ability  to  plot  skymaps  for  each  of 
four  frequency-range  bin«  in  a  case,  side  by  side,  to  separate  the 
information  from  different  heights,  as  shown  in  Figure  6.  When 
these  are  stacked  on  a  page,  the  temporal  development  of  the 
skymaps  also  becomes  apparent.  Only  sources  whose  amplitudes  are 
within  10  dB  of  the  maximum  amplitude  for  each  subcase  are 
printed,  as  indicated  on  the  right  side  of  each  row.  This  type  of 
display  is  useful,  for  example,  in  data  consistency  checks  and  for 
ionospheric  tilt  studies. 
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■  The  calculation  of  the  three  dimensional  drift  velocity 
vector  from  the  Doppler  sources  has  remained  unchanged.  Assuming 
a  uniform  velocity  field  v,  the  vector  is  determined  by  a  least- 
squares  process  including  all  qualified  sources  s  =  1,  2,  S.  The 
error  to  be  minimized  is: 


2 

e 


7C 


IA3 


(34) 


where  ksx  and  ksy  specify  the  location  for  source  s  with  Doppler  ds, 
and  k^sz  =  (2nlX)‘^  -  k^sx  -  k^sy;  the  weights  Ag  are  (linear)  source 
amplitudes  obtained  from  a  combination  of  the  amplitudes  measured 
at  each  of  the  four  (or  seven)  antennas. 

The  PO  method  gives  a  higher  precision  for  the  source 
locations.  With  the  FWPD  scanning  method,  the  source  locations  were 
quantized:  in  order  to  avoid  a  prohibitive  number  of  calculations,  the 
map  representing  the  sky  above  the  Digisonde  station  was  divided 
into  a  41-by-41  set  of  pixels,  and  the  FWPD  was  calculated  only  at 
these  1681  pixel  locations.  Sources  between  these  locations  ended 
up  at  the  closest  pixel  in  the  map  array,  and  only  one  source  per 
pixel  could  be  kept  in  the  map  array  when  there  were  several 
sources  close  together.  The  analytical  PO  method  is  not  limited  to 
quantized  locations.  It  keeps  the  same  visual  map  representation  of 
the  source  locations  (coordinates  (X,Y)  =  0,  ±1,  ±2,  ...,  ±20),  but  for 
the  drift  velocity  calculations  the  coordinates  are  given  with  one 
decimal  place  instead  of  being  integers.  The  accuracy  of  our 
measurements  does  not  justify  keeping  more  than  one  decimal  place. 
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2.2.5  Special  Processing  Concerns 

The  new  drift  processing  software  has  achieved  two 
objectives:  (1)  reduction  in  CPU  time,  (2)  improved  data  validation. 
The  time  required  to  process  one  magnetic  tape  containing  the  5 
minute  data  for  one  day  is  approximately  500  CPU  seconds,  a  factor 
of  20  improvement  over  the  FWPD  method. 

Identification  of  "good"  sources  is  most  important  if  valid 
drift  velocities  are  to  be  derived.  .  The  measured  drift  data  consist  of 
four  or  seven  64-line  spectra,  one  for  each  receiving  antenna.  Each 
spectrum  is  thresholded  by  zeroing  those  spectral  lines  which  have 
an  amplitude  below  the  most-probable-amplitude  of  the  spectrum. 
Any  such  "noise"  lines  are  not  used  for  the  drift  velocity  calculation. 
The  remaining  spectral  lines  are  considered  potential  Doppler  sources 
unless  a  spectral  amplitude  varies  from  antenna  to  antenna  by  20  dB 
or  more.  The  latter  case  indicates  that  more  than  one  source  exists 
with  the  same  Doppler  shift,  and  therefore  the  spectral  line  is  not 
used. 

The  phases  of  the  surviving  set  of  S  spectral  lines  are 
27C-corrected,  and  then  checked  for  consistency.  Because  of  the 
geometry  of  the  antenna  array,  the  corrected  phases  must  fulfill  the 
conditions  [Bibl,  1983]: 

tit  f 

(inner  antennas) 

and  (35) 

t  t  I  I 

^0  (outer  antennas) 

These  equations  are  easily  verified  by  inspection  of  Figure  4: 


(ks*r2+6)+(ksT3+5)+(ksT4-^)-  3(ks*r  1+5)  =  ks«(r2+r3+r4) +  35  -  35 

=  0  (36) 

since  r2  +  r3  +  r4  =  0 
and  ri  =  0. 

Because  of  noise  and  round-out  errors,  these  sums  are  generally  net 
exactly  zero.  If  a  A  is  larger  than  80°,  the  spectral  line  is  not  used 
for  skymap  and  drift  calculations.  Testing  the  distribution  of  A's  is  a 
convenient  way  of  checking  the  proper  operation  of  the  receiving 
array.  A  defective  antenna  results  in  systematically  large  values  for 
A. 


2.2.6  Validation  of  New  Phase  Only  Analysis  Technique 

To  verify  that  the  new  analysis  program,  based  on  the  PO 
technique  is  performing  correctly,  three  days  of  data  from  Qaanaaq 
were  processed  twice  using  alternately  the  proven  FWPD  based 
program  and  the  new  PO  based  program.  Figure  7  shows  the  results 
of  this  test  for  drift  measurements  made  during  the  October- 
November  1986  Viking  campaign.  The  velocities  calculated  from,  the 
FWPD  skymaps  are  joined  by  a  dashed  line.  It  is  clear  from  this 
comparison  that  the  two  methods  of  calculating  the  source  positions 
are  virtually  equivalent. 
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3.0  OBSERVATIONS  AND  CAMPAIGNS 


3.1  Drifts 


Development  of  the  instruments  at  the  three  Digisonde 
sites  Argentia,  Goose  Bay,  and  Qaanaaq  has  progressed  at  different 
rates  for  a  number  of  reasons  including  previous  history, 
accessibility,  and  equipment  failure.  Consequently,  the  amount  and 
quality  of  the  data  obtained  from  each  site  has  also  varied.  Due  to 
the  large  quantities  of  data  collected  by  these  instruments,  and  the 
availability  of  different  modes  of  operation,  a  computerized 
inventory  technique  was  developed. 

Tables  1  to  3  in  Appendix  A  depict  the  drift  observing 
schedules  for  Argenfia  (operating  autodrift  since  October  1986), 
Goose  Bay  (autodrift  since  January  1988  when  the  DGS  256  replaced 
the  DGS  128PS),  and  Qaanaaq  /.ntodrift  since  October  1986)  for  the 
years  1986,  1987,  1988  and  198i>'.  Table  lA  is  also  shown  on  the 
next  page.  Each  line  on  the  table  includes  a  "time  line"  which  is 
divided  into  288  columns  (one  column  for  each  five  minutes  of  a  24- 
hour  period),  preceded  by  the  date  and  "WD"  for  the  world  days.  On 
the  time  line,  the  start  time  of  each  group  of  continuous  drift 
measurements  (rounded  to  5-minute  increments)  is  indicated  by  a 
tic  mark.  To  ease  reading  the  inventory  over  a  day,  the  tic  marks  are 
longer  on  the  hour  (3  dots  instead  of  2),  \>’ith  an  additional  fourth  dot 
every  6  hours.  The  "p"  following  some  of  the  time  lines  indicates  the 
days  whose  data  has  been  processed  into  drift  velocity  graphs.  All  of 
the  data  from  all  three  stations  has  been  processed,  except  for  about 
12  months  of  Argentia.  The  latter  was  omitted  because  of  a 
hardware  fault  at  the  station,  discovered  in  August  1989,  which  had 
started  about  one  year  earlier,  and  resulted  in  erroneous  values 
being  recorded  on  the  tapes.  Note  that  the  ionograms  recorded  in 
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this  period  are  correct.  The  hardware  problem  could  not  be  fixed 
until  March  1989,  when  routine  processing  of  the  data  resumed. 

In  general,  the  drift  mode  operated  at  each  station  on  the 
three  World  Days  of  every  month  as  specified  in  the  contract.  The 
sounding  schedule  for  the  World  Days  was  one  ionogram  every  five 
minutes,  with  drifts  measured  between  ionograms.  In  addition  to 
the  contractual  requirements,  drift  data  were  also  taken  ''n  a  15 
minute  schedule  on  many  other  days.  Thus  a  substantial  data  base 
of  drift  measurements  has  been  compiled,  parts  of  which  have  been 
used  in  the  research  leading  to  the  publications  listed  in  Sections  5.0 
and  Appendix  B. 


3.2  Oblique  Ionogram  Measurements 

The  objective  of  this  part  of  the  project  was  to  develop 
and  use  oblique  ionospheric  sounding  between  existing  Digisonde 
(DISS)  stations  in  order  to  increase  the  effective  density  of 
ionospheric  sensors.  This  is  clearly  necessary  to  monitor  the 
occurrence,  location  and  extent  of  the  F  layer  trough.  The  intent  was 
to  use  the  vertical  ionograms  from  Argentia  and  Goose  Bay  together 
with  oblique  ionograms  between  these  two  stations  for  the  trough 
monitoring.  The  Goose  Bay  Digisonde  256  (DISS)  was  not  installed 
until  1988,  therefore  the  initial  oblique  propagation  experiments 
were  conducted  between  Argentia  and  Millstone  Hill.  Sample  data 
for  the  Goose-Argentia  link  are  also  now  available  for  limited  periods 
and  they  show  that  oblique  monitoring  is  a  valuable  tool,  as 
discussed  below. 
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3.2.1  Oblique  Sounding  Between  Goose  Bay  and  Argentia 

In  a  feasibility  study,  good  echoes  were  obtained  at 
Argentia  during  most  of  the  times  of  day  (see  Figures  8a  to  8f)  in 
spite  of  strong  interference  from  Europe  at  many  frequency  bands. 
For  most  of  the  recordings,  both  sounders  transmitted  and  received 
simultaneously.  This  was  done  for  testing  purposes,  but  is  not 
generally  advisable  because  of  the  large  difference  in  signal 
amplitudes  between  the  vertical  and  oblique  signals.  At  low 
frequencies  the  automatic  receiver  gain  and  printer  gain  squelch  the 
oblique  signals  because  the  vertical  echoes  are  so  strong.  Conversely, 
the  vertical  signals  near  foF2  are  gain-reduced,  because  the  oblique 
signals  are  strong.  However,  the  variations  in  the  maximum 
observable  frequency  (MOF)  are  clearly  identified,  and  can  be  u.sed 
to  aid  in  monitoring  the  trough  region. 


3.2.2  Oblique  Sounding  Between  Argentia  and  Millstone 

In  conjunction  with  another  experiment,  oblique 
propagation  between  the  University  of  Lowell  Digisonde  256  at 
Millstone  Hill,  Massachusetts  and  the  Air  Weather  Service  Digisonde 
256  at  Argentia,  Newfoundland  was  carried  out.  This  path  was  used 
to  establish  the  network  characteristics  of  monitoring  the  reflection 
points  between  established  ionosonde  stations  as  proof  of  concept  for 
routine  oblique  propagation  experiments  to  define  mid-path 
ionospheric  characteristics. 

Data  exist  from  both  stations  on  day  87293  (October  20, 
1987)  from  0018  to  0348  UT  (Tables  2A  and  2B)  and  from  day 
87296  (October  23,  1987)  1033  UT  to  day  87297  (October  24,  1987) 
0433  UT  (Tables  3  and  4). 

At  least  two  time  periods  cover  oblique  ionograms 
recorded  simultaneously  on  both  stations  which  establishes  the 


31 


gure  8.  Reception  in  Argentia  of  Oblique  lonograms  from  Goose  Bay 


Figure  8.  Reception  in  Argentia  of  Oblique  lonograms  from  Goose  Bay  (Continued) 


Figure  8.  Reception  in  ^entia  of  Oblique  lonograms  from  Goose  Bay  (Continued) 
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Figure  8.  Reception  in  Argentia  of  Oblique  lonograms  from  Goose  Bay  (Continued) 


Table  2A.  Oblique  Propagation  (from  Argentia)  Received  at  Millstone 


Table  2B. 


TIME 

UT 

0018 

0033 

0048 

0103 

0118 

0133 

0148 

0203 

0318 

0333 

0348 


Day  293,  October  20,  1987 


TIME 

foF2 

IF  MUF 

2F  MUF 

UT 

(MHz) 

(MHz) 

(MHz) 

0018 

5.25 

11.20 

0033 

5.20 

10.80 

0048 

5.20 

10.60 

U9.30 

0103 

4.90 

10.50 

U3.30 

0118 

4.60 

9.50 

7.80 

0133 

4.50 

9.20 

0148 

4.10 

9.25 

0203 

3.90 

Note:  "U"  values  are  uncertain. 


Oblique  Propagation  (from  Millstone)  Received  at  Argentia 


foF2 

IE  MUF 

IF  MUF 

2FMUF 

3F  MUF 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

4.60 

11.25 

4.60 

11.30 

10.80 

4.45 

11.60 

10.60 

8.40 

6.90 

4.30 

10.50 

10.50 

8.40 

4.00 

11.60 

9.30 

7.80 

4.00 

11.50 

9.15 

7.80 

3.90 

11.60 

9.25 

3.60 

3.30 

10.70 

7.10 

5.70 

4.80 

3.10 

11.60 

6.80 

5.80 

11.60 

6.70 
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Table  3A. 


Oblique  Propagation  (from  Argentia)  Received  at  Millstone 


Day  296,  October  23,  1987 


TIME 

foF2 

lEMUF 

IF  MUF 

2FMUF 

UT 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

1029 

6.90 

1033 

4.10 

7.80 

5.05 

1044 

4.60 

8.50 

1059 

4.80 

7.30 

10.40 

6.25 

1103 

U5.00 

10.90 

7.00 

1114 

5.00 

1118 

>11.50 

1129 

5.10 

1148 

10.10 

13.00 

1159 

5.40 

1203 

1214 

5.50 

11.40 

12.65 

8.60 

1218 

1229 

5.50 

11.10 

13.80 

8.70 

1233 

11.70 

14.25 

9.05 

REMARKS 


Direct 

Backscatter 


Change  Tx 
Antenna 
Configura¬ 
tion 


Change  Tx 

Antenna 

Orientation 
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Table  3B.  Oblique  Propagation  (from  Millstone)  Received  at  Argentia 


Day  296,  October  23,  1987 


TIME 

foF2 

lEMUF 

IF  i'4UF 

2FMUF 

3F  MUF 

UT 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

1018 

2.35 

1029 

2.30 

7.850 

1048 

2.70 

6.2  to 
(X)9.300 

1103 

3.75 

(X)10.850 

7.000 

4.80 

1118 

4.00 

>11.500 

(0)7.200 

(X)7.600 

1148 

4.70 

9.10  13.075 

(0)7.870 

(X)8.250 

(X)6.95 

1203 

5.00 

12.650 

(0)8.200 

(X)8.650 

1218 

5.40 

13.800 

(0)8.300 

1233 

6.00 

14.250 

(X)9.000 

1248 

6.30 

15.000 

(X)9.450 

1303 

6.70 

(F3) 

1318 

6.80 

14.750 

U9.250 

7.75 

1333 

6.5/ 

7.20 

1348 

6.70 

14.750 

9.475 

1403 

6.80 

>16.000 

10.000 

+  Nose 

1418 

6.90 

»16.000 

10.950 

9.20 

1433 

7.10 

» 16.000 

11.600 

U9.40 

1448 

7.30 

>16.000 

11.475 

9.20 

REMARKS 

Direct 

Backscatter 

First 

Oblique 

Echoes 


Whip  only 
(poor) 


Large 

Rhombic 


Strong 

Direct 

Backscatter 
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Table  4A.  Oblique  Propagation  (from  Argentia)  Received  at  Millstone 


Day  296,  October  23,  1987 


TIME 

foF2 

lEMUF 

IF  MUF 

UT 

(MHz) 

(MHz) 

(MHz) 

2103 

7.20 

6.80 

16.30 

2110 

13.45 

2118 

7.10 

15.40 

2133 

6.90 

13.10 

2148 

6.50 

9.00 

12.40 

2203 

6.20 

12.50 

2218 

5.90 

12.45 

2233 

5.90 

10.30 

12.90 

2248 

6.20 

9.10 

13.30 

2303 

5.90 

12.50 

2318 

5.00 

11.10 

2333 

5.40 

10.70 

2348 

5.30 

10.25 

*Note: 

values 

are  uncertain 

2F  MUF 

3F  MUF 

(MHz) 

(MHz) 

REMARKS 

10.30 

8.20 

Loops 

9.20 

8.20 

8.60 

6.90 

Long  High 

(0)7.90 

*U6.7 

Ray 

E  +  F  Mode 

(0)7.80 

*U6.5 

Present 

(0)7.85 

(X)8.30 

(X)7.50 

(X)8.40 

Strong  E  +  F 

(X)8.05 

(0)6.40 

Mode 

E  +  F  Mode 

(X)7.70 

Present 
Long  High 

6.80 

(0)5.40 

*U5.5 

Ray,  E  +  F 
Mode 
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Table  4B.  Oblique  Propagation  (from  Argentia)  Received  at  Millstone 


Day  191  y  October  24,  1987 


TIME  foF2  lEMUF  IF  MUF  2F  MUF  3F  MUF 

UT  (MHz)  (MHz)  (MHz)  (MHz)  (MHz)  REMARKS 


0003 

5.20 

9.40 

0018 

5.10 

10.10 

0033 

4.90 

9.40 

0048 

4.70 

8.70 

0103 

4.70 

(X)9.20 

0118 

4.30 

(X)8.10 

0133 

4.00 

(X)8.10 

0148 

3.90 

(X)7.90 

0203 

3.90 

(X)7.60 

0218 

4.00 

(X)8.10 

0233 

3.90 

(X)7.70 

0248 

3.90 

7.30 

0303 

3.80 

7.10 

0318 

3.80 

(X)7.20 

0333 

3.70 

(X)7.10 

0348 

3.70 

(X)6.85 

0403 

3.80 

(X)6.85 

0418 

3.70 

(X)6.95 

0433 

3.60 

6.90  (X)6.90 

(.) 

(•) 

Spread  F 

(.) 

(X)4.90 

(0)5.60 

(X)4.90 

(0)5.40 

(0)5.20 

5.20 

5.20 

5.05 

4.90 

4.90 

Nose 

4.80 

Extension 

8.1 

Nose  to  7.6 

4.70 

(0)4.50 

Best  lono- 

(0)4.50 

(0)4.40 

(X)4.81 

(X)4.80 

(0)4.40 

gram. 

Good  2F 
Mode 

Note:  (  .  )  indicates  mode  is  present  but  MUF  can't  be  scaled. 


absolute  group  path  between  the  stations.  This  occurred  on  day 
87293  0018  UT  to  0203  UT  and  on  day  87296  from  1033  UT  to  1233 
m  (Table  5). 

As  expected,  these  periods  show  the  same  MOF  on  both 
records,  but  the  vertical  critical  frequencies  at  the  endpoint  stations 
were  often  quite  different.  As  Figure  9  shows,  the  transmission 
factor  varies  substantially  between  morning  and  evening  hours.  This 
is  explained  by  the  change  in  E-region  ionization  and  the  change  in 
the  profile  of  the  F-region.  But  the  variation  of  the  transmission 
factor  in  each  of  these  time  periods  is  due  to  the  change  in  local 
gradients  of  the  maximum  electron  density.  Thus  it  is  inappropriate 
to  linearly  average  foF2  between  stations  more  than  1500  km  apart. 


3.2.3  Oblique  Sounding  Between  Goose  Bay  and  Millstone 

Because  of  its  excellent  synchronization  stability  the  best 
and  most  reliable  oblique  propagation  data  were  obtained  from  the 
Goose  Bay  -  Millstone  link  (see  Figure  10).  Although  it  was  not 
considered  necessary  for  the  study  of  the  ionospheric  trough  this  link 
will  improve  the  validity  of  the  Argentia  data.  The  necessary  time 
transfer  by  more  than  one  hour  of  the  Argentia  data  must  be 
verified  as  feasible.  If  the  motion  of  the  trough  can  be  distinguished 
from  other  time-dependent  changes,  then  the  study  of  the  three 
station  net  is  sufficient.  Also  the  lateral  extension  of  the  trough  can 
be  investigated. 


3.2.4  Trough  Studies 

A  substantial  number  of  ionograms  from  the  four 
stations:  Goose  Bay,  Labrador;  Argentia,  Newfoundland;  Millstone 
Hill,  Massachusetts  and  Wallops  Island,  Virginia,  have  been 
presented  as  ionogram  surveys  for  the  study  of  the  F  region  trough 
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Table  5A.  Oblique  Propagation  Received  Simultaneously  at  Millstone  and  Argentia 
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Table  5B.  Oblique  Propagation  Received  Simultaneously  at  Millstone  and  Argentia  (Continued) 
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LOCAL  TIME  (ARGENTIA)  =  UT  -  3  hr  30  min 
LOCAL  TIME  (MILLSTONE  HILL)  =  UT  -  4  hr  46  min 


Figure  9.  Variability  of  the  Virtual  Transmission  Factor  (M-Fac'tor) 


/T.ssurnin§  i^incar  oraOicnts  in  lOx  r^r^cntici 
Day  87293 


JlTXlAi.O  VVrilV/ 


46 


Figure  10.  Goose  Bay  to  Millstone  Oblique  Bistatic  lonograms 


between  the  auroral  and  mid-latitude  ionosphere.  Table  6  gives  a 
schedule  of  World  Day  data  from  1988  and  1989  when  ionogram 

surveys  were  available  for  at  least  three  of  the  four  stations.  Many 

more  days  will  be  available  from  1989  and  in  the  future. 

Analysis  of  digital  data  has  the  advantage  of  better 

discrimination  of  the  vertical  echo  traces  and  of  the  ordinary 
component.  Compared  to  the  dense  station  net  [Whalen,  1989], 

available  during  the  International  Geophysical  Year  1958,  the  chain 
of  available  stations  is  very  limited.  To  improve  the  situation  an 
effort  will  be  made  for  automatic  bistatic  sounding  between  the 
ionosondes  to  monitor  the  midpoint  conditions.  Note  that  Goose  Bay 
data  generally  appear  to  be  necessary  to  draw  a  cross-section  of  the 
maximum  F-layer  critical  frequency,  foF2,  as  a  function  of  time  of 
day. 


One  example  of  trough  development  is  presented  here. 
At  the  beginning  of  March  1989  (see  Figures  11a  and  11b)  an 
intensive  measurement  campaign  led  to  good  data  from  the  most 
significant  stations.  On  Day  061  the  ionospheric  trough  clearly 
reached  Millstone  Hill,  MA  and  had  its  minimum  between  06  and  07 
UT,  while  before  that  time  first  Argentia  and  later  Millstone  show 
additional  echo  traces  of  reflection  from  the  trough  walls.  On  Day 
065  Argentia  shows  echo  traces  from  the  walls  of  the  trough  since 
0100  UT  all  the  time  until  1000  UT,  while  the  Millstone  ionosphere 
breaks  down  almost  completely  at  0700  UT. 


In  the  future  edited  ARTIST  data  will  be  used  for  surveys 
of  the  trough  onsets.  The  ADEP  [Zhang,  1989]  software  package 
(ARTIST  Data  Editing  and  Printing)  permits  printing  of  the  critical 
frequencies  foF2,  foE  and  foEs,  and  the  electron  density  contours  as 
function  of  time.  An  example  is  shown  in  Figure  12  for  6  March 
1989,  the  trough  minimum  occurs  at  23  UT  at  Goose  Bay,  04  UT  at 
Argentia  and  06  UT  at  Millstone. 
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Table  6.  Available  Survey  lonograms 


mSi 

Goose  Bay 

Argentia 

Millstone 

Wallops  Is. 

75 

X 

X 

0  Bad 

X 

76 

X 

X  Absorption 

X  08-10  Z 

Very  low  foF2 

X 

77 

X 

0 

X 

X 

138 

X  02-06  Z 

Aurora  E 

X  06Z 

Low  foF2 

X  06Z 

Low  foF2 

X 

139 

X  Aurora  E 

All  night 

X  06-07  Z 

Low  foF2 

X  08Z 

Low  foF2 

X 

140 

X  07-15  Z 
foF2  and  FI  F2 

0 

X  08Z 

foF2  missing 

X 

166 

X  Aurora  Es 

X  0435  Z 
foF2  missing 

X 

X 

167 

X 

X 

X 

X 

194 

X  830-930  Z 
Absorption 

Aurora  01-07  Z 

0 

X 

X 

195 

X  No  morning 
data 

0 

X 

X 

196 

X 

X 

X 

xs 

287 

X 

X  Day  missing 

X 

X 

313 

X 

Absorption 

X  09-10  Z 
foF2  2.5  MHz 

X  02-07  Z 
foF2  2.5  MHz 

X 

314 

X 

X 

X 

X 

315 

X 

X 

X 

X 

348 

X 

X 

0 

X 

im. 

61 

X 

X 

X 

65 

X 

X 

X 
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MOTION  OF  THOUGH  OVER  THREE  STATION  CHAIN  1989  061 

Figure  11a.  Motion  of  Trough  Over  Three  Station  Chain 


50 


Figure  11b.  Motion  of  Trough  Over  Three  Station  Chain 
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Figurf  12.  Trough  Onset  6  March  1989 


L2 _ Campaigns 

Campaign  studies  have  formed  an  integral  part  of  this 
contract.  The  AWS  Digisonde  network  is  an  important  high  latitude 
diagnostic  available  to  both  the  geophysics  and  radar  engineering 
communities.  Several  of  the  campaigns  involving  the  Digisondes  are 
summarized  below  (chronological  order). 

a)  Viking/GL  (October  28  -  November  15,  1986) 

This  campaign  focussed  on  the  electrodynamics  of  polar  cap 
sun-aligned  F-region  arcs,  and  velocity  shears  contiguous  with 
arcs  in  the  auroral  oval.  Diagnostics  included  the  Viking  and 
HILAT  satellites,  EISCAT  and  Sondrestrom  radars,  GL  All  Sky 
Imaging  Photometers  (ASIP's),  and  AWS  Digisondes.  The  GL 
Airborne  Ionospheric  Observatory  (AIO)  was  based  at  Thule 
early  on  in  the  campa^^n,  and  operated  a  Digisonde  together 
with  satellite  scintillation  receivers  (250  MHz),  ASIP  and 
spectrometers.  The  AIO  Digisonde  operated  in  a  2.5  minute 
sequence  during  ground  based  and  airborne  observations. 
From  28  October  to  3  November,  the  Digisondes  at  Qaanaaq  and 
Argentia  operated  in  a  5  min  ionogram/drift  sequence.  On  3 
November,  the  AIO  <raTisferi'ed  to  Andoya,  Norway,  and 
participated  in  joint  observations  with  EISCAT  in  Tromso. 
.Analysis  of  this  data  set  cortributed  to  a  paper  currently  in 
preparation  by  Lars  Blaii  [1989  private  communication]. 

b)  Slant  Bs  (July  26  -  August  7.  1987) 

A  new  direction-scanning  ionosonde  mode  and  the  drift  mode 
were  implemented  at  Qaanaaq  for  this  campaign,  which  was  a 
collaborative  effort  between  GL,  ULCAR  and  DMI.  In  this  mode, 
each  frequency  sounding  within  the  ionogram  is  repeated  6  (or 
12)  times  while  the  beam  of  the  receiving  antenna  array  scans 
through  6  (or  12)  azimuth  directions  in  60°  (or  30°)  increments. 
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Analysis  of  slant  Es  data  from  this  campaign  and  the  one  in 
1986  resulted  in  a  paper  by  Hoeg  et  al.  [1987]  that  was 
presented  at  the  Fall  (1987)  AGU  meeting  in  San  Francisco  (see 
Appendix  B). 

Simultaneously  with  the  slant  Es,  F  layer  patches  were  also 
studied.  The  patches  often  produced  an  interesting  signature 
in  the  ionograms,  consisting  of  a  "nose"  similar  to  that  observed 
at  oblique  sounding  (see  Figure  8).  The  peculiar  triangular 
shape  in  the  spread  F  X-echoes  on  the  ionograms  (Figures  13 
and  14)  is  caused  by  the  oblique  reflections  from  a  patch, 
forming  the  nose  that  is  typical  for  oblique  ionograms.  With 
plasma  frequencies  estimated  to  be  only  about  10%  higher  than 
the  background,  these  patches  return  signals  from  distances 
maybe  up  to  200  km,  resulting  in  a  very  flat  nose.  Often  the 
nose  of  the  0-trace  is  covered  by  X-echoes  and  is  therefore  not 
visible.  But  at  1456  UT  (Figure  13)  the  0-nose  is  also  visible. 

(c)  Structure  and  Dynamics  of  the  Auroral  Ionosphere  (January  8  - 
15,  1988) 

The  GL-AIO  flew  to  the  north  of  the  EISCAT  radar  and 
performed  250  MHz  phase  and  amplitude  scintillation 
measurements,  recorded  all-sky  auroral  images  at  6300  A  and 
4278  A  and  performed  ionospheric  sounding  using  a  Digisonde. 
These  measurements  were  coordinated  with  EISCAT 
observations,  and  DMSP-F8  overpasses  in  which  precipitating 
electron  and  ion  fluxes  were  measured,  together  with  in-situ 
thermal  plasma  parameters  at  850  km  altitude. 

Large  scale  ionospheric  features  such  as  F-layer  blobs  and  the 
trough  were  observed.  The  location  and  intensity  of  small  scale 
irregularities  was  mapped.  All  ionograms  for  the  six  AIO 
missions  on  January  8,  9,  15,  16,  17  and  18  have  been  scaled  in 
terms  of  foF2,  foE,  foEs,  and  receiver  gain  (showing  the  level  of 
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Figure  14.  Qaanaaq,  Greenland  Digital  lonogram  14  August  1987, 
1706  UT 
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absorption;  high  gain  #  =  high  absorption).  Table  7  lists  these 
parameters  as  a  function  of  UT  and  aircraft  location  in  terms  of 
the  corrected  geomagnetic  latitude  (CGLAT).  The  ionization  in 
the  F  region  rarely  went  above  1  x  10^  cm"^  (3  MHz  plasma 
frequency).  This  is  comparable  to  the  ionization  measured  in 
the  E  region.  The  critical  frequencies  for  the  auroral  E  are 
generally  just  below  3  MHz,  and  sporadic  E  layers  which  are 
present  most  of  the  time  vary  between  3  and  9  MHz.  In  this 
environment  it  is  difficult  to  derive  vertical  electron  density 
profiles  for  the  F  region. 

Figure  15  shows  three  AIO  ionograms  from  the  18  January 
flight  at  1809,  1816  and  1819  UT.  The  F  layer  is  absent  or 
blanketed  in  the  last  two  ionograms,  so  a  profile  inversion  was 
attempted  only  for  the  1809  ionogram.  The  figure  shows  three 
different  electron  density  profiles  computed  for  1809  UT  using 
different  assumptions  about  the  electron  distribution  in  the 
transition  (valley)  region  between  the  E  and  F  layers.  The 
different  valleys  shown  in  the  profile  plots  affect  the 
calculation  of  the  hmF2  values,  causing  them  to  vary  from  236 
km  for  the  valley  shown  in  the  top  panel  to  227  km  (middle) 
and  to  209  km  (bottom).  In  the  latter  case  the  transition 
region  was  assumed  to  be  filled  with  ionization. 


d)  CEDAR-HLPS  (February  6  -  24,  1988) 

The  CEDAR-High  Latitude  Plasma  Structures  (HLPS)  initiative 
has  as  its  objective  the  study  of  the  source,  evolution  and 
ultimate  fate  of  ionospheric  irregularities  at  high  latitudes. 
Digisondes  obtained  drift  data  and  ionograms  at  Qaanaaq,  Goose 
Bay,  and  Argentia.  Time  resolution  was  generally  15  minutes, 
but  for  two  core  periods  (11-12  Feb.  21-02  UT;  18/19  Feb.  21- 
02  UT)  it  was  five  minutes.  Ionograms  were  also  recorded  by 

me  lyigi&uiiuc  uu  uu^iu  me  paiis.cu  at  iiiuic.  lauie 


57 


8  lists  the  available  data  for  Goose  Bay,  Qaanaaq,  Millstone  Hill 
and  Argentia. 

Many  other  instruments  supported  this  campaign,  including 
ISRs,  and  ASIPs. 

e)  CEDAR-HLPS  (December  5  -  13,  1988) 

The  format  of  this  HLPS  campaign  was  very  similar  to  that 
described  above. 

Results  from  the  analysis  of  the  CEDAR-HLPS  data  have  been 
presented  at  the  NSF  CEDAR  meetings  in  Boulder  1988  and 
1989,  at  the  URSI  meeting  in  Boulder  in  1989  and  Spring  1989 
AGU  meeting  (see  Section  5.0). 

f)  Slant-Es  and  Drift  Measurements  at  Sondrestrom  (July  7-11, 
1989) 

Drift  mode  operations  were  in  effect  at  all  three  Air  Force 
Digisondes:  Qaanaaq,  Goose  Bay  and  Argentia.  Unfortunately, 
the  Sondrestrom  ISR  was  also  undergoing  maintenance/repairs 
and  only  1-2  hours  of  simultaneous  observation  was  possible. 
Nature  did  not  cooperate  and  the  expected  Slant-Es  traces  were 
not  observed  due  to  a  lack  of  significant  magnetic  activity. 

However,  this  trip  was  extremely  valuable  because  it  produced 
five  days  of  drift  data  from  the  new  AWS  sounder  at 
Sondrestrom:  These  data  are  discussed  in  Section  4.0  and  yield 
some  interesting  geophysics  in  addition  to  confirming  the 
correct  operation  of  the  new  sounder.  The  experiment  also 
revealed  a  problem  with  interference  to  communication  systems 
which  was  resolved  in  November  1989. 
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UIGISOIIDE  OI’linATlUlJ  IM  SUPPOIIT  OF  CEDAR  HLPS  l'J08 


Table  8.  Digisonde  Data  in  Support  o 

(Febru  ary  1988). 


4.0  SONDRESTROM  DRIFT  OBSERVATIONS 


The  drift  data  obtained  at  Sondrestrom  during  the  period 
89188  -  89192  are  illustrated  in  Figure  16.  In  the  midnight  sector, 
the  drifts  are  generally  close  to  antisunward,  as  expected  for  polar 
cap  flows  under  Bz  southward  conditions.  The  antisunward  direction 
is  indicated  by  an  oblique  line  in  the  frame  marked  Az.  Magnetic 
midnight  is  indicated  by  MM.  Systematic  offsets  from  straight 
antisunward,  as  seen  e.g.  00-05  UT  on  88-190,  is  likely  due  to  the 
influence  of  By  on  the  exact  flow  direction,  as  shown  for  the  Qaanaaq 
data  [Appendix  B,  Cannon  et  al.,  1989].  Around  06  UT,  the  station 
passes  through  the  convection  reversal  into  a  region  of  sunward 
(eastward)  auroral  flow.  This  remains  strongly  eastward  until  the 
noon  sector  when  the  direction  switches  through  northward  to 
strongly  westward.  In  the  late  afternoon  (18  UT)  the  convection 
reversal  is  crossed  again,  as  Sondrestrom  re-enters  the  polar  cap. 
These  changes  are  fairly  abrupt  in  the  data  of  Figure  16.  In  only  one 
case,  (Day  188,  18-24  UT)  there  is  no  clear  transition  from  auroral  to 
polar  cap  flow. 

The  drift  speeds  also  varied  dramatically  during  the  five 
day  period  plotted  in  Figure  16.  There  is  no  dependence  on  time  of 
day,  although  both  the  horizontal  and  vertical  speeds  measured  in 
the  polar  cap  tend  to  be  larger  than  those  in  the  auroral  return  flow. 
There  is  also  no  clear  dependence  on  magnetic  activity  as  defined  by 
Kp  (shown  along  the  lower  axis  of  the  figure). 

The  most  readily  available  corroborative  data  are  the 
drift  data  from  Qaanaaq  (Figure  17).  Qaanaaq  is  always  inside  the 
polar  cap,  and  the  nature  of  these  measurements  provides  a  strong 
indication  of  IMF  conditions  [see  Appendix  B,  Cannon  et  al.,  1989]. 
Again,  the  oblique  lines  plotted  in  the  Frame  marked  Azimuth 
indicate  the  antisunward  drift  direction.  The  preliminary  Kp  indices 
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Figure  16.  Sondrestrom  Drift  Data,  89188  to  898192 
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are  shown  along  the  lower  axis  of  the  figure,  and  vary  from  1  to  3. 
The  drift  directions  tend  to  be  closer  to  antisun  ward  for  Kp  =  3,  while 
there  is  much  more  variation  for  lower  Kp. 

Since  Sondrestrom  and  Qaanaaq  are  at  approximately  the 
same  corrected  magnetic  local  time  (Figure  2),  a  polar  plot  in 
corrected  magnetic  coordinates  provides  an  interesting  opportunity 
to  explore  the  large  scale  features  of  convection  for  that  magnetic 
longitude.  Figure  18  depicts  the  drift  velocities  observed  at 
Sondrestrom  and  Qaanaaq  as  a  function  of  CGLT  for  the  24  hour 
period  between  06  UT  Day  89191  and  0545  UT  Day  89192.  This  was 
one  of  the  most  active  intervals  of  the  entire  data  set,  according  to 
the  Kp  index,  and  data  from  Sondrestrom  and  Qaanaaq  were 
available  almost  continuously. 

Several  interesting  features  emerge  from  examination  of 
Figure  18  (note  that  UT  ~  CGLT  +  2  hours).  At  the  start  of  the 
interval,  the  Qaanaaq  drifts  diverge  from  the  antisunward  direction, 
becoming  almost  sunward  by  09  CGLT,  suggesting  that  the  IMF  Bz 
component  is  turning  northward  (Bz  >  0).  The  Sondrestrom  drifts  are 
very  variable  during  this  interval,  though  sunward  flows 
predominate  from  about  07  CGLT.  At  0915  CGLT,  the  Qaanaaq  flows 
change  al  3tly  to  antisunward,  with  high  speeds  around  600  ms'^, 
typical  ot  IMF  Bz  southward.  The  Sondrestrom  flow  also  becomes 
antisunward  before  10  CGLT,  suggesting  that  plasma  is  entering  the 
polar  cap  rapidly  from  the  dayside. 

At  Qaanaaq,  the  drift  direction  remains  antisunward  for 
the  entire  day,  until  the  last  vector  plotted  at  0545,  indicating  that 
the  IMF  Bz  component  remained  southward  until  0545.  The 
Sondrestrom  drift  variations  are  more  complex.  Using  the  result  of 
Cannon  et  al.  (1989)  which  is  consistent  with  Heppner  and  Maynard 
(1987)  for  the  influence  of  By  on  polar  cap  flows,  it  is  possible  to 
predict  the  sign  of  the  IMF  By  component  from  the  Qaanaaq  flow 
direction.  On  average.  By  positive  leads  to  deviations  of  36°  anti- 
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cloc?:‘vise  from  the  antisunward  direction:  By  negative  leads  to 
de-vi'Uions  of  -12°.  Figure  19  depicts  the  predicted  By  direction 
bas?d  on  the  Qaanaaq  data.  For  much  of  the  time,  the  predictions 
bdped  on  Sondrestrom  data  are  the  same.  For  example  at  03  UT  the 
flow  direction  changes  from  approximately  southward  to  strongly 
south  westward,  indicating  an  IMF  By  change  from  southward  to 
northward.  When  IMF  data  become  available  in  the  near  future, 
these  predictions  will  be  tested. 

Unfortunately  software  is  not  yet  available  for  plotting 
figures  such  as  Figure  18  on  a  routine  basis.  Since  these  plots 
promise  to  be  extremely  useful,  the  appropriate  software  will  be 
developed  before  other  days  are  examined  in  detail. 


5.0  SCffiNTIHC  RESULTS 


5.1  Published,  Papers 

The  research  resulted  in  nine  published  papers  as  listed: 

1.  Basu,  Sa.,  S.  Basu,  C.  E.  Valladares,  E.  J.  Weber,  J.  Buchau,  G.  J. 
Bishop  and  B.  W.  Reinisch,  "Coordinated  Observations  of  High 
Latitude  Ionospheric  Turbulence,"  Physics  of  Space  Plasmas 
(1988),  SPI  Conference  Proceedings  and  Reprint  Series,  Vol.  8, 
pp.  137-151,  1988. 

2.  Bossy,  L.,  R.  R.  Gamache  and  B.  W.  Reinisch,  "Lay  Functions  for 
F2  Profiles,"  Adv.  Space  Res.,  Vol.  8,  No.  4,  pp.  4(201)-4(204), 
1988. 

3.  Buchau,  J.,  B.  W.  Reii.'  :h,  D.  N.  Anderson,  E.  J.  Weber  and  C.  G. 
Dozois,  "Polar  cap  plasma  convection  measurements  and  their 
relevance  to  the  modeling  of  the  high-latitude  ionosphere," 
Radio  Sci.,  Vol.  23,  No.  4,  pp.  521-536,  July- August  1988. 

4.  Cannon,  P.  S.,  B.  W.  Reinisch,  J.  Buchau  and  T.  W.  Bullett,  "A 
Statistical  Validation  of  the  Digisonde  Technique  for  Polar  Cap 
F-Region  Convection  Measurements,"  submitted  to  JGR  1989. 

5.  McNamara,  L.  F.,  B.  W.  Reinisch  and  J.  S.  Tang,  "Values  of  hmF2 
Deduced  from  Automatically  Scaled  lonograms,"  Adv.  Space 
Res.,  Vol.  7,  No.  6,  pp.  (6)53-(6)56,  1987. 

6.  Reinisch,  B.  W.,  J.  Buchau  and  E.  J.  Weber,  "Digital  lonosonde 
Observations  of  the  Polar  Cap  F  Region  Convection,"  Physica 
Scripta,  Vol.  36,  pp.  372-377,  1987. 

7.  Reinisch,  B.  W.,  R.  R.  Gamache,  X.  Huang  and  L.  F.  McNamara, 
"Real  Time  Electron  Density  Profiles  from  lonograms,"  Adv. 
Space  Res.,  Vol.  8,  No.  4,  pp.  (4)63-(4)72,  1988. 
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8. 


Reinisch,  B.  W.,  J.  Buchau,  K.  Bibl  and  G.  S.  Sales,  "Multistation/ 
Multiparameter  Observations  with  a  Network  of  Digital 
lonosondes,"  AGARD  Conference  Proceedings  No.  441,  pp.  37-1 
-  37-9,  1988. 

9.  Reinisch,  B.  W.,  K.  Bibl,  G.  S.  Sales,  J.  Tang,  Z.-M.  Zhang,  T.  W. 
Bullett  and  D.  F.  Kitrosser,  "The  Digisonde  256  Ionospheric 
Sounder,"  SCOSTEP  WITS  Handbook  #2,  November  1989. 

Appendix  B  of  this  report  contains  copies  of  papers  3  and 
4;  other  papers  were  contained  in  Scientific  Reports  #1  and  #2. 


5.2  Scientific  Reports 

Two  Scientific  Reports  have  been  published: 

1.  Reinisch,  B.  W.  J.  Buchau,  E.  J.  Weber,  L.  F.  McNamara  and  J.  S. 
Tang,  "Ionospheric  Convection  and  Structure  Using  Ground- 
Based  Digital  lonosondes,"  Scientific  Report  No.  1,  AFGL-TR-88- 
0042,  ULRF-440/CAR,  February  1988.  ADA194246 

2.  Reinisch,  B.  W.,  J.  Buchau,  R.  R.  Gamache,  K.  Bibl,  G.  S.  Sales,  X. 

Huang  and  L.  F.  McNamara,  "Electron  Density  Profiles  and 
Plasma  Drift  Measurements  with  Digital  lonosondes,"  Scientific 
Report  No.  2,  AFGL-TR-88-0233,  ULRF-442/CAR,  ADA205193, 
September  1988.  ada205193 


5.3  Conference  Presentation.^ 

Presentations  were  made  at  a  number  of  conferences: 

1.  Buchau,  J.,  D.  N.  Anderson,  E.  J.  Weber,  B.  W.  Reinisch  and  C.  G. 
Dozois,  "Polar  Cap  Plasma  Convection  and  Real  Time  Modeling," 
Ionospheric  Effects  Symposium,  Springfield,  VA,  May  1987. 
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2.  Buchau,  J.,  B.  W.  Reinisch,  K.  Bibl  and  D.  F.  Kitrosser,  "New 
Research  Opportunities  with  a  Network  of  Digital  lonosondes," 
URSI  General  Assembly,  Tel  Aviv,  Israel,  1987. 

3.  Cannon,  P.  S.,  J.  Buchau,  B.  W.  Reinisch,  "Statistical  Validation  of 
Digisonde  Drift  Data  from  Qaanaaq,"  AGU,  San  Francisco,  CA, 
December  1989. 

4.  Cannon,  P.  S.,  J.  Buchau,  B.  W.  Reinisch,  G.  Crowley  and  R. 
Lepping,  "IMF  Controlled  Variations  in  the  Convection  Flow 
Direction  Near  the  Corrected  Geomagnetic  Pole,"  US  National 
URSI  Meeting,  Boulder,  CO,  January  1990. 

5.  Crowley,  G.  and  B.  W.  Reinisch,  "Qaanaaq  Drift  Observations  in 
Support  of  the  December  1988  CEDAR/HLPS  Campaign,"  CEDAR 
Meeting,  Boulder,  CO,  June  1989. 

6.  Crowley,  G.,  J.  Buchau  and  B.  W.  Reinisch,. "First  Results  from  the 
Digisonde  at  Sondrestrc.  US  National  URSI  Meeting,  Boulder, 
CO,  January  1990. 

7.  Hoeg,  P.,  J.  Buchau,  B.  W.  Reinisch  and  J.  K.  Olesen,  "Drift 
Measurements  of  E-Region  Irregularities  in  the  Polar  Cap,"  AGU, 
San  Francisco,  CA,  EOS  December  1987. 

8.  Reinisch,  B.  W.,  J.  Buchau  and  E.  J.  Weber,  "Polar  Cap  F  Region 
Plasma  Convection,"  URSI  General  Assembly,  Tel  Aviv,  Israel, 
1987. 

9.  Reinisch,  B.  W.  and  J.  Buchau,  "Qaanaaq  Drift  Observations  in 
Support  of  the  January  88  CEDAR/HLPS  Campaign,"  CEDAR 
Meeting,  Boulder,  CO,  January  1988. 

10.  Reinisch,  ■  .  W.,  J.  Buchau  and  P.  S.  Cannon,  "Digital  lonosonde 
Measurements  of  High  Latitude  Ionospheric  Convection,"  lAGA, 
Exeter,  UK,  July  1989. 

1 1 .  Reinisch,  B.  W.,  J.  Buchau  and  C.  G.  Dozois,"  Correlating  Polar  Cap 
Plasma  Drift  with  IMF,"  US  National  URSI  Meeting,  Boulder,  CO, 
January  1989. 
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12.  Reinisch,  B.  W.,  J.  Buchau  and  C.  G.  Dozois,  "Response  of  the  Polar 
Cap  Convection  to  Changes  in  the  IMF,"  AGU  Meeting,  Baltimore, 
MD,  May  1989. 

13.  Sales,  G.  S.,  B.  W.  Reinisch  and  C.  G.  Dozois,  "A  New  Measure  of 
Ionospheric  Roughness,'  URSI  General  Assembly,  Tel  Aviv, 
Israel,  1987. 

Papers  #4,  #7  and  #11  are  contained  in  Appendix  B. 


72 


6.0  REFERENCES 


Anderson,  D.  N.,  J.  Buchan  and  R.  A.  Heelis,  Origins  of  density 
enhancements  in  the  winter  polar  cap  ionosphere,  Radio  Sci.,  23,  No. 
4,  1988. 

Bibl,  K.  and  B.  W.  Reinisch,  The  Universal  Digital  lonosonde.  Radio  Sci., 
13,  p.  519,.1978. 

Bibl,  K  esolves  Multiple  Reflectors 

with  ^  ..iultiplexed  Closely  Spaced 

Frequei  .nmission  F  1983  Symposium, 

Louvain,  Belgium,  p^  •  ^,  June  1983. 

Bibl,  K.,  B.  W.  Reinisch  and  D.  F.  Kitrosser,  Digisonde  256  -  General 
description  of  the  compact  digital  ionospheric  sounder.  University  of 
Lowell  Center  for  Atmospheric  Research,  1981. 

Bilitza,  D.,  N.  M.  Sheik  and  R.  Eyfrig,  A  global  model  for  the  height  of 
the  F2  peak  using  M3000  values  from  the  CCIR  numerical  maps. 
Telecomm.  J.,  46,  pp.  549-553,  1979. 

Bossy,  L.,  R,  R.  Gamache  and  B.  W.  Reinisch,  Lay  Functions  for  F2 
Profiles,  Adv.  Space  Res.,  Vol.  8,  No.  4,  pp.  4(201)-4(204),  1988. 

Buchau,  J.,  B.  W.  Reinisch,  E.  J.  Weber  and  J.  G.  Moore,  Structure  and 
dynamics  of  the  winter  polar  cap  F  region.  Radio  Sci.,  18(6),  p.  995, 
1983. 

Buchau,  J.,  B.  W.  Reinisch,  D.  N.  Anderson,  E.  J.  Weber  and  C.  Dozois, 
Polar  cap  plasma  convection  measurements  and  their  relevance  to 
the  modeling  of  the  high-latitude  ionosphere.  Radio  Sci.,  23,  No.  4,  pp. 
521-536,  1988. 

Bullett,  T.  W.,  Thesis,  University  of  Lowell,  1990. 


73 


Cannon,  P.  S.,  B.  W.  Reinisch,  J.  Buchau  and  T.  W.  Bullett,  "A  Statistical 
Validation  of  the  Digisonde  Technique  for  Polar  Cap  F-Region 
Convection  Measurements,"  submitted  to  JGR  1989. 

Cauffman,  D.  P.  and  D.  A.  Gurnett,  Satellite  measurements  of  high- 
latitude  convection  field,  Space  Sci.  Rev.,  13,  p.  369,  1972. 

Dozois,  C.  G.,  A  high  frequency  radio  technique  for  measuring  plasma 
drifts  in  the  ionosphere.  Rep.  AFGL-TR-83-0202,  Air  Force  Geophys. 
Lab.,  Hanscom  AFB,  MA,  ADA140509,  1983. 

Dudeney,  J.  R.,  The  accuracy  of  simple  methods  for  determining  the 
height  of  the  maximum  electron  concentration  of  the  F2  layer  from 
scaled  ionospheric  characteristics,  J.  Atmos.  Terr.  Phys.,  45,  pp.  629- 
640,  1983. 

Dunge/,  J.  W.,  Interplanetary  magnetic  fie’d  and  the  auroral  zones, 
Phys.  Rev.  Lett.,  6,  pp.  47-48,  1961. 

Hakura,  Y.,  Tables  and  maps  of  geomagnetic  coordinates  corrected  by 
the  higher  order  spherical  harmonics  terms.  Rep.  lonos.  Space  Res. 
Jpn.,  19,  p.  121,  1965. 

Heppner,  J.  P.,  Empirical  models  of  high-latitude  electric  fields,  J. 
Geophys.  Res.,  82,  7,  pp.  1115-1125,  1977. 

Heppner,  J.  P.  and  N.  C.  Maynard,  Empirical  high-latitude  electric  field 
models,  J.  Geophys.  R'--.-.,  92,.A5,  pp.  4467-4489,  1987. 

Hoeg,  P.,  J.  Buchau,  B.  W.  Reinisch  and  J.  K.  Oleson,  Drift 
Measurements  of  E-Region  Irregularities  in  the  Polar  Cap,  AGU,  San 
Francisco,  CA,  EOS,  December  1987. 

Kelso,  J.  M.,  Radio  Ray  Propagation  in  the  Ionosphere,  McGraw-Hill, 
1964. 

Reinisch,  B.  W.,  New  Techniques  in  Ground-Based  Ionospheric 
Sounding  and  Studies,  Radio  Sci.,  Vol.  21,  No.  3,  pp.  331-341,  May- 
June  1986. 


Reinisch,  B.  W,,  J.  Buchau,  K.  Bibl  and  G.  S.  Sales,  Multistation/ 
Multiparameter  Observations  with  a  Network  of  Digital  lonosondes, 
AGARD  Conference  Proceedings  No.  441,  pp,  37-1  -  37-9,  1988b. 

Reinisch,  B.  W.,  J.  Buchau  and  E.  J.  Weber,  Digital  ionosonde 
observations  of  the  polar  cap  F  region  convection,  Phys.  Scr.,  36,  p. 
372,  1987. 

Reinisch,  B.  W.,  R.  R.  Gamache,  X.  Huang  and  L.  F.  McNamara,  Real 
Time  Electron  Density  Profiles  from  lonograms,  Adv.  Space  Res.,  Vol. 
8,  No.  4,  pp.  (4)63-(4)72,  1988a. 

Titheridge,  J.  E.,  lonogram  Analysis  with  the  Generalized  Program 
POLAN,  World  Data  Center  A  for  Solar  Terrestrial  Physics,  Rep.  UAG- 
93,  1985. 

Weber,  E.  J.,  J.  Buchau,  J.  G.  Moore,  J.  R.  Sharber,  R.  C.  Livingston,  J.  D. 
Winningham  and  B.  W.  Reinisch,  F  layer  ionization  patches  in  the 
polar  cap,  J.  Geophys.  Res.,  89(A3),  p.  1683,  1984. 

Whalen,  J.  A.,  Auroral  oval  plotter  and  nomograph  for  determining 
corrected  geomagnetic  local  time,  AFCRL-70-0422,  ADA713170, 
1970. 

Whalen,  J.  A.,  The  daytime  F  layer  trough  and  its  relation  to 
ionospheric-magnetospheric  convection,  J.  Geophys.  Res.,  94,  A12,  pp. 
17,169-17,184,  1989. 

Zhang,  Z.-M.,  ARTIST  Data  Editing  and  Printing,  ULCAR  Report,  1967. 


75 


APPENDIX  A 

Invent  of  Drift  Observations 
from 

Argentia,  Goose  Bay  and  Qaanaaq 


77 


1987 


10  11  le  13  lA  IS  18  17  18  19  ftO  ftl  ftfi  fi3 


17 

WD  870A8 

18 

UO  870A9 

19 

UO  87050 

ftO 

87031 

Oi 

A70Sfi 

ftft 

87033 

&3 

870SA 

S4 

87033 

83 

87036 

88 

87037 

3 

8706A 

& 

87063 

7 

87066 

8 

87067 

9 

87068 

10 

87069 

11 

87070 

18 

87071 

IS 

a707A 

16 

87073 

17 

87076 

lA 

87077 

19 

87078 

8 

87098 

3 

87093 

A 

8709A 

18 

87108 

13 

87103 

lA 

8710A 

IS 

87103 

17 

87107 

18 

87108 

19 

87109 

80 

87110 

87 

87117 

88 

87118 

89 

87119 

30 

87180 

1 

87181 

8 

87188 

3 

87183 

A 

871&A 

S 

87183 

6 

87186 

7 

87187 

8 

87188 

9 

87189 

10 

87130 

11 

87131 

18 

87138 

13 

87133 

lA 

8713A 

IS 

87133 

16 

87136 

81 

UD  871A1 

88 

871A8 

83 

871  A3 

8A 

871AA 

83 

87  IAS 

86 

871A6 

87 

871A7 

88 

871A8 

89 

871A9 

30 

87130 

31 

87131 

1 

87138 

8 

87133 

A 

87133 

3 

87136 

6 

87137 

7 

87138 

8 

87139 

9 

87160 

10 

87161 

It  t •  It  it  I « 


I  I  I  I  I  *  I  I  I  I  I  (  I  i 
lllltllllllll 


I  *  «  I  I  I  I  I  *  I  I  I  I  I  (  I  »  i 
I  I  I  I  *  I  I  I  *  I  I  I  I  i  I  I  I  i 


■  III!  I  I  I  I  I 


«  I  I 
•  1  • 
•  •  I  I 


•  .  t  • 
t  t  I  • 


I  I  i  • 
I  I  I  • 


till 
•  •  I  . 


HOURLY 


1/  Itt  19 


IS-^XNLTre 


lO-HlNUTC 


S-MXNUTC  HCAeuRCnCNTQ 


TABLE  lAi  ARQENTXA  DRIFT  SURVEY 


79 


aa  aaoaaaa  .  aaaaaaaaaa  aaaaaaaa  aaaaa  aaa  aaaa  aaaa  aaaaaaaaanaoaaaaaaar.  aaaaaaaaaaaaa  a&aoaaa 


1987  (cont. > 

O  1  ft  i  4  S  6  7  8  9  lO  11  Ifi  ta  14  IS  14  17  18  19  eO  ftl  ftfi  63  U.T. 


JUN  Ifi 

87183 

13 

87184 

14 

8718S 

IS 

87188 

18  UD 

87187 

17  UD 

87168 

18  WD 

87189 

19 

87170 

60 

87171 

61 

07176 

66 

87173 

63 

87174 

£4 

87173 

6S 

87178 

68 

87177 

£7 

87178 

68 

87179 

69 

87180 

JUL  7 

87188 

8 

87189 

9 

87190 

10 

87191 

11 

8719C 

16 

87193 

13 

87194 

14 

8719S 

63  WD 

87604 

£4 

8760S 

£S 

87608 

68 

87607 

67 

87608 

68 

a76or 

69 

87610 

30 

87611 

31 

87616 

6U0  1 

87613 

7 

07619 

8 

87660 

9 

87661 

10 

87666 

11 

87663 

le 

87664 

13 

87663 

14 

87668 

17 

87669 

18  UD 

87630 

19  UD 

87631 

60  UD 

87636 

61 

87633 

66 

67634 

63 

87633 

64 

87638 

6S 

67637 

68 

87630 

67 

87639 

68 

87640 

69 

07641 

30 

87646 

31 

87643 

8CP  1 

87644 

6 

87643 

3 

87848 

4 

87647 

S 

07648 

8 

87849 

7 

87630 

8 

87631 

9 

87636 

10 

876S3 

13 

87638 

14 

87637 

IS 

87630 

18 

07639 

17 

07680 

18 

87681 

19 

87686 

66  UD 

87683 

63  UD 

87688 

64  UD 

67687 

es 

87688 

68 

87689 

67 

87670 

68 

87671 

69 

07676 

OCT  13 

87688 

14 

87607 

IS 

87688 

18 

07689 

17 

07690 

18 

87691 

60  UD 

87693 

61  WD 

87694 

66  UD 

07693 

63 

07698 

I  •  •  •  I  >  I  I  I  .  ,  i  I  I  ,  .  I  ,  «  ,  I  t  , 
I  •  •  •  I  •  •  t  I  .  I  .  I  ■  .  ,  I  .  I  i  .  I 

I . I . Ill 

I  >  >  •  t  •  .  •  t  .  >  I  I  •  .  I  4  •  •  *  I  I  > 


.  1  .  . 
•  I  t  . 
.  I  .  . 

.  t  •  * 


(■til 

I  •  1  •  I 
I  •  *  •  I 


I 

I 

I 

I 


I 

I 

I 

I 


. . . . I 

I  •  •  t  I  •  I  I  I  I  •  I  I  •  t  I  I  i  •  •  I  I  «  •  I  (  t  t  I  I  •  »  I  »  »  t  I  •  ■  I  I  •  »  •  I  •  «  ■  t  ■  «  <  <  •  <  <  i  •  *  •  I  •  i  <  I  •  •  I  I  i  •  •  I 

I  •  •  •  I  •  •  t  I  •  •  •  I  •  I  i  I  •  •  •  I  •  t  •  I  •  •  »  I  *  <  ■  I  *  *  «  1  t  t  «  I  «  f  •  I  ■  «  *  I  ■  «  «  I  t  «  »  I  •  i  •  I  i  f  I  I  I  *  *  I  I  I  i  I 

1  *  •  •  I  I  •  I  I  ■  •  i  I  •  •  •  I  ■  •  ■  I  I  t  I  I  •  I  •  I  *  ■  •  I  ■  •  »  I  t  «  I  1  I  I  •  I  •  <  1  I  •  «  *  I  t  •  •  I  I  ■  I  I  I  I  •  I  i  •  »  I  •  •  i  I 

I  •  I  •  I  I  I  i  I  •  •  i  I  •  •  •  I  I  •  I  I  I  •  I  I  I  •  f  I  »  •  •  I  •  I  •  I  »  »  «  t  •  »  (  I 

I  i  •  «  I  *  i  I  I  •  i  •  I  I  i  I  i  I  I  I  I  I  I  •  j  •  t  *  I  ■  «  (  »  k  t  *  I  >  (  I  I  t  •  «  I  •  <  «  I  t  t  <  I  *  i  <  I  <  i  •  I  i  I  I  I  I  •  I  I  *  *  I  j 

I  •  I  •  I  *  I  i  I . I  *  •  •  I  I  t  «  I  »  «  *  I  »  *  t  I  k  •  I  I  *  «  ■  1  »  t  ■  I  »  t  1  I  1  i  *  I  *  t  4  I  f  I  •  I  I  •  I  I  I  •  •  I  I  ■  i  I 

I . I  I  4  4  I  •  4  i  t  I  •  4  I  k  k  •  t  •  k  »  I  «  k  k  I  k  *  k  I  k  I 


I  4  I  4  I 
I  4  4  4  I 
I  4  4  4  I 
I  4  4  4  I 


I  4  4  4  I 

t  4  4  4  I 

I  4  4  4  I 


klkkkl444lk44l444l4k4|k44l444lk44l44kl 

|4l4l44>l44ll4l4lll4lll4|4k4lkkklkkklkkklkkklk44|kkktkk4l4kkl444lt44l444|444l444ll44l444l 

|444l444l444lk44l4l4l4tk|4kklkkklkkklkk4lkkklktk|kkklk4kt444l444l444l444j44ll444l444l44ll 

I•••l*••l•44l444i444ll44jkkklkkkl4kklkkklkkklk44|k4klk4k14(4l444l444l444^44tl444ll44t444l 

|444l444l444l4l4l44ki444jk4klkkklk4kl4kklkkklk4k|kkklk44t444l444l444l444|4k4i4kkl444l444l 

I*•*l•••l•••l••■ll•4lk44i4k4lkkklkk4lkk4l4kklk44ikkki4kki•k4ll44lk44l444|4k4l4l4l444l44tl 


. I  4  4  4  I  4  .  4  4  4  1 

l•••l••<••••)*••i•i•l••4|4kklkkkl444l44llkk4lkkk|4k4lklk)444l444i444l444|444l444l4l4l444l444l 
I . I . .  I  k  4  4  I  4  4 


I . . . . 

I . I . . . I  4  •  4  >  4  4  .  4  4  4  1  k  4  4  1  k  4  k  I  4  4  4  I  4  4  4  . . I  4  4  4  |  4  4  4  I  4  4  4  I . I  4  4  4  I 

I  •  •  •  I . . . I  k  I  4  I . I  4  4  4  I  4  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  >  4  4  4  i  4  4  4  I 

I  4  4  •  I  4  4  I  4  4  4  I  4  I . I  4  4  4  |  4  4  4  I  4  4  4  (  4  4  4  I  k  4  4  I  k  k  4  I  k  k  k  {  4  4  k  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  |  4  4  4  I . .  4  4  4  I  4  4  4  I 

I  •  . . . . I  4  4 . I  4  4  4  I  4  4  4  (  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  f  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  »  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I  4  4  4  I 

1  4  4  4  I  4  4  4  t  4  4  4 


I  4  4  4  I 
i  4  4  4  I 

)  4  4  I  I 
I  4  4  4  I 
I  4  4  I  I 
I  4  4  4  i 


I  4  4  4  I 
I  4  4  4  I 


I  4  4  4  I 


I  4  4  4  I 


I  4  4  t  I 


I4l4l444|444l444l4f4l444l444l 

|4>4ll44l444l4f4l44«l44t|444l4t4l  44(4t4l44klktt|444l444l444l4f4l444l444|444l444l444l444l444l 

|414ltl4l44il444l44tl444|l4/l444l4k4l444t4k4lktl|4tklk4ll444l444l444l444l444l444l444l444l444l 

|4>ll4ttll4>l>4ll4>4l444|444l4l4l4(4l44kl444l444|f4tl444l444lf44l444l444|444l444l444i444l444l 

|4lll4t4lf4ltt44l444i444|44fl44fl444l4k4l44kl4k4|44ttk4*l444lf44l444l444|4|ll444l444l444l444l 


I  4  . . I 


>  4  I  4  I 


I  4  4  4  I  4  I  4  I 
I  4  4  4  I  t  4  4  I 
I  4  4  4  I  4  4  4 


I  4  4  4  I  4  4  4  I  4  t  4  I  4  4  4  I  4  f  4  1  4  4  4  I 


. . 

•  4  4  i  I . 

I  4  4  I  I  4  4  4  I  •  4  I  I  4  4  4  I  4  4  4  I  4  4  4  t 

*  4  4  ,  1  4  4  .  1  . I  4  4  4  I 


4l444l444|444l444l444l444l444l 

|4lll4tll444ll4il4t4itf4l444l444lf44l4#4P444|44k|4kk|4tlll44ll44l444ljl4t444l444l444|444l444l 

|444l4l4l444l44lll4llf44|444l444l444l444l444l4kf|t44i44kll44ll44|444l 


I  I  •  •  •  >  1  .  I  .  •  .  I  .  4  4  .  .  I  ,  4  I  .  4  4  ,  ,  I  ,  4  4  .  4  I  M  ,  4  4  I  4  4  4  4  4  1  .  4  4  4  «  t  k  4  4  4  4  I  4  4  4  4  4  » 
J  .  .1  4  4.  4  .rf  4  4  4  .  4  I  4  4  44 . I . .  4  4  4  4  4  I .  4  4  4  4  I  4  4  4  4  4  I  4  4.  4  .1 


.  1  4  4  4  4  .  1  . 


.4.|. 


P 

P 

D 

P 


LCOCNDf 


6  J  4 

I  i  I 

MOURUY 


S  8  7  a  9  10  11  le  13  lA  IS  18  17  18  19  60  ftl  fifi  63  U.T. 

•  I  I  4  4  i  4  . .  I  4  4  4  I  4  4  4  . .  I  4  4  4  4  .  I . I . . . . . . . . 

IS-HtNUTC  lO-HlNUTC  S-MINUTC  HCAdUftOCNTS 


TOBLE  IB I  PRQENTIP  DRIFT  SURVEY 


80 


aaaaaaaaaaaaaaaaao  aaaaaaa  aaaaaaaaaa  aaaaaaa  aaa.aaaaaaaaaa.aaaa(iaaaaaaa  aaaaaaa  aaaaaaaa  aaaaaa 


1987  (cont. > 

9  lo  it  le  13  14  IS  IS  17  1ft  19  eo  ei  ea 


OCT  ftft  ft7ft99 

fi7  ft7300  I  • 

ftft  ft7301  I  • 


NOV  ft4  WD  a73&ft 
ftS  WD  ft73ft9 
ftfi  WD  ft7330 


OCC  17  ft73Sl 
1ft  ft73SA 
19  ft73S3 
ftO  ft73S4 
fit  ft733S 
fift  WD  ft73Sft 
63  WD  87337 
fi4  WD  A733A 


19BS 

9  ID  11  1ft  13  14  IS  1ft  17  lA  19  ftO  fil 


83  U.T. 


ft 

ftft037 

i  • 

7 

•ft03A 

1  • 

ft 

Aft039 

1  • 

9 

ftft040 

1  • 

10 

Aft041 

11 

A*04ft 

1ft 

ftft043 

1  I 

13 

Aft044 

1  • 

1ft 

WO 

A*047 

17 

WD 

Aft04ft 

1ft 

WO 

AA049 

K— 

19 

AftOSO 

84 

AAOSS 

83 

AftOSft 

1  • 

ftft 

ftftOST 

i  • 

87 

•ftOSft 

i  « 

8ft 

ftftOS9 

1  * 

89 

AAOftO 

1  • 

1 

ftftOftl 

1  • 

ft 

ftftOftft 

1  • 

3 

Aftofta 

4 

A*0ft4 

1  « 

S 

AftOftS 

I  « 

ft 

ftftOftft 

1  • 

7 

**0ft7 

1  • 

ft 

ftftOftft 

1  * 

9 

ftftOftft 

1  • 

11 

AA07i 

18 

ftft078 

1  • 

13 

aA073 

1  • 

14 

AA074 

1  « 

IS 

WO 

ftftOTS 

1ft 

WD 

ftft076 

i  « 

17 

WD 

ftA077 

83 

AAOA3 

84 

ftftOA4 

1  * 

83 

AftOftS 

1  • 

ftft 

AftOftA 

1  • 

87 

ftaoa7 

1  i 

8A 

ftftOAA 

1  • 

89 

AA0A9 

1  < 

30 

Aft090 

1  * 

3 

Aft094 

4 

AftOftS 

1  • 

S 

AAOftft 

1  • 

ft 

ftftOft7 

1  • 

7 

AA09a 

1  • 

ft 

ftftOftft 

1  • 

9 

AAIOO 

1  • 

10 

AAlOl 

1  > 

11 

AAlOft 

1  • 

14 

AAIOS 

1  « 

3 

Aaifi4 

5 

AAlfift 

ft 

ftai&7 

7 

AAICA 

A 

AAlCft 

V 

AAiao 

to 

AA131 

11 

AA138 

18 

Aai33 

13 

Bai34 

14 

AAiaS 

IS 

AA13ft 

1ft 

ftftl37 

17 

WD 

Aai3A 

1ft 

WD 

Aai39 

19 

WD 

AA140 

a  9  10  11  1C  13  14  IS  1ft  17  la  19  ftO  Cl  CC  &3  U.T. 

>lM  tfcl.llPlIlillW.llllhl.ai.MlI  — NMII  »  . . 

IS-HZNUrC  lO^ZNUTC  S-NXNUTC  KCAftURCMCNTft 


TRBLE  1C I  RRQENTIft  DRIFT  SURVEY 


il 


81 


isaa 

9  iO  tl  Ifi  13  14  IS  16  17  la  19  90  91 


HAR 


7 

66036 

6 

66039 

9 

66040 

10 

66041 

11 

66049 

19 

66043 

13 

66044 

14 

6604S 

IS 

66046 

16 

WD 

66047 

17 

WD 

66046 

16 

WD 

66049 

19 

660S0 

90 

660S1 

91 

660S9 

14 

66074 

IS 

WD 

6607S 

16 

WD 

66076 

17 

WD 

66077 

16 

66076 

19 

WD 

66110 

64 

66115 

93 

66116 

96 

66117 

97 

66116 

16 

66137 

17 

WD 

66136 

16 

WD 

66139 

19 

WD 

66140 

90 

66141 

14 

WD 

66166 

13 

WD 

66167 

16 

WD 

66166 

11 

66193 

16 

WD 

66194 

13 

WO 

66193 

14 

WD 

66196 

IS 

66197 

66 

66939 

67 

66940 

66 

66941 

99 

66949 

30 

66943 

13 

WO 

66937 

14 

WD 

66936 

IS 

WD 

66939 

10 

66964 

11 

WD 

66965 

16 

WD 

13 

WD 

66967 

14 

6flgOB 

6 

WO 

66313 

9 

WD 

66314 

10 

WD 

66315 

11 

66316 

6 

66337 

3 

66336 

S 

66340 

6 

66341 

7 

66349 

a 

66343 

9 

66344 

10 

66345 

11 

66346 

16 

66347 

13 

WD 

66346 

14 

WD 

66349 

15 

WO 

663 SO 

16 

66331 

Mlttltt*** 


•  . . 

. . . . . 


‘I* 


. . 


.  . . . 


**I«***4*IM>»I< 


. . . Ill 


4  4  •  4  I  4  4  4  I  4  I  I  I  I  4  4 
4  4  4  4  4  4  i  4  4  4  4  4  i  4  4  4  4  4 
4  4  4  I  4  4  4  4  4  1  4  4  4  4  4  I  4  4  4  4  4 


JAN  16  69016 

17  UD  69017  i< 
16  WD  69016  !• 
19  WD  69019  t' 


1969 

9  10  11  16  13  lA  IS  16  17  16  19  60  ftl  66  63 


ISoMXNUTC 


16 
.  1.4 


13  14  IS  16  17  16  19  60  61  66  63 


lO-MlNUTC 


S-NXNUT6  nCMmjRSnCNTtt 


TABLE  2Ai  GOOSE  BAY  DRIFT  SURVEY 


84 


19B9  <cont. ) 


10  11  18  13  14  IS  16  17  16  19  80  81  88  83 


FC8  13  69044 

14  UD  6904S  I* 

15  WD  69046 


e 

3 

4 

s 

6 

7 

6 

9 

10 

11 


69061 

69066 

69063 


69066 

69067 
69066 

69069 

69070 


1 1  1 1 1  Btn  •  MM  ti  l  mm  Imm 

. 

1  I  1 

* 

M«Ml«(MllMMII«M*lliMMllMM 

1  ‘  ■  ‘1 

13  69076 

14  UD  69073 


16  UD  6907S 


17 

16 

19 

60 


69076 

69077 
69076 
69079 


69100 


11  UD  69101 
16  UD  69108 

13  UD  69103 

14  69104 


HAV  9  UD  69169 

10  UD  69130 

11  UD  69131 


Mli  i 

MMlM— 

j  .  .  . 

, _ 

i  '  '  ‘ 

j  .  . 

. 

M  M  • 

JUN  16  69163 

13  UD  69164 

14  UD  6916S 

15  UD  69166 

16  69167 


JUL  10  69191 

11  UD  69198 
16  UD  69193 

13  UD  69194 

14  69193 


LCOCNDl 


6  3 

I  I 

HOURLY 


6 

i » • 


7 

•  lit 


0  10 
1 1  I  •  f  1 1 

1SH4XNUTC 


11  16  13  14  IS  16  17  16 

I  I  I  I  I  I  I^MMi 

10-441NUTE 


19  60  61  66  83  U.  T. 

S<^tNUTC  MCAOURCHENTO 


TABLE  EBi 


GOOSE  BAY  DRIFT  SURVEY 


85 


aaaao.  aaaaaaaaaa  aaaaaaaa  aaaaa  aao.  aaaaa  aaaaa 


APPENDIX  B 

Copies  of  Publications 


95 


For  Presentation  at  the  U.S.  National  URSI  Meeting  at  Boulder,  CO, 

4-6  January  1989 


Correlating  Polar  Cap  Plasma  Drift  with  IMF 

Bodo  W.  Reinischl,  Jurgen  Buchau^,  and  Claude  G.  Dozois^ 

1  University  of  Lowell  Center  for  Atmospheric  Research, 

450  Aiken  Street,  Lowell,  Massachusetts  01854 

2 Air  Force  Geophysics  Laboratory,  Hanscom  Air  Force  Base, 
Bedford,  Massachusetts  01731 


Abstract 


More  than  one  year  of  polar  cap  drift  data  of  the  F  region  have  been 
analyzed.  The  measurements  were  made  at  Qaanaaq,  Greenland 
(87°N  CGL)  with  a  Digisonde  256  operated  by  the  Danish 
Meteorological  Institute.  Drift  velocities  between  150  and  900 
meters  per  second  were  observed.  For  several  days  in  1987  and  for 
February  1988,  interplanetary  magnetic  field  (IMF)  data  were 
available  and  the  direction  of  the  F  region  drift  was  compared  with 
the  north-south  component  Bz  of  the  IMF.  For  Bz  negative,  the  drift 
is  almost  always  antisunward;  for  Bz  positive,  it  is  irregular  and  often 
sunward.  The  response  time  of  the  ionosphere  to  sign  changes  in  Bz 
is  fast.  Although  fluctuations  in  Bz  and  the  drift  directions  make  it 
difficult  to  define  a  precise  time,  one  can  see  a  response  often  within 
less  than  30  minutes. 

Introduction 


Since  1986  F  region  drift  measurements  were  made  in  Qaanaaq, 
North  Greenland  (87°  N  CGL)  with  a  Digisonde  256,  operated  by  the 
Danish  Meteorological  Institute  in  cooperation  with  the  Air  Force 

vju/vpaiyoxwo  auu  tiic  uiixvoiMiy  Ol  i-»OWcil.  1  lie  general 


state  of  the  ionosphere  was  observed  with  ionograms  that  were 
made  every  5  or  15  minutes.  Between  the  ionograms  the  Digisonde 
measured  the  F  region  drift  using  two  different  sounding 
frequencies,  sampling  the  F  layer  at  heights  that  are  50  to  150  km 
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apart.  The  sequence  of  ionogram/drift  measurements  is  illustrated 
in  Figure  1,  which  shows  the  on-line  printout  of  an  ionogram 
followed  by  eight  consecutive  drift  measurements.  The  Digisonde 
drift  technique  has  been  previously  described  /I,  2/.  This  paper 
discusses  some  drift  results  for  periods  for  which  interplanetary 
magnetic  field  (IMF)  data  could  be  obtained  from  the  National  Space 
Science  Data  Center. 

The  Polar  Cap  Plasma  Convection 

The  existence  of  polar  plasma  convection  was  first  inferred  from 
satellite-borne  electric  field  measurements  /3/  and  plasma  drift  meter 
measurements  /4/.  Ground-based  observations  with  the  incoherent 
scatter  radar  (ISR)  installations  at  Chatanika  /5/  and  Millstone  Hill  /6/ 
also  determined  the  convection  pattern.  The  Digisonde  256  can  easily 
determine  the  plasma  drift  as  function  of  time  over  long  time  periods. 
This  allows  to  study  in  detail  the  effects  of  the  IMF  on  the  plasma 
convection  for  those  times  where  IMF  data  are  available. 

Figure  2  shows  Heelis's  model  /?/  of  the  convection  pattern  in  a 
CGL/CGLT  coordinate  system.  The  Qaanaaq  observatory  moves  around 
the  magnetic  north  pole  in  24  hours  on  the  circle  marked  Q.  The 
antisunward  direction  at  Qaanaaq  varies  linearly  with  time,  and  since 
the  plasma  flow  is  expected  to  be  antisunward  for  Bz  south  conditions  , 
our  data  plots  always  show  this  direction. 

Figure  3  shows  the  Qaanaaq  drift  data  for  03  -05  March  1987, 
together  with  the  z-component  Bz  of  the  IMF.  The  data  are 
presented  as  function  of  universal  time.  The  direction  of  the  plasma 
drift  is  given  in  corrected  geomagnetic  coordinates,  with  south  in  the 
center  of  the  panel,  east  is  below  and  west  above.  The  antisunward 
(asw)  direction  is  indicated  by  the  straight  line.  The  magnitude  of 
the  drift  vector,  shown  in  the  top  panel,  is  given  in  m/s.  For  large 
data  sets,  this  presentation  is  well  suited  for  studies  of  the 
correlation  between  the  orientation  of  the  IMF  and  the  direction  of 
the  plasma  flow.  It  is  not  possible  to  determine  the  entire  convection 
pattern  with  just  one  observing  station.  It  would  require  at  least  10 
sounders  strategically  distributed  to  determine  the  shape  of  the 
convection  pattern  as  function  of  time. 
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Case  Studies 


To  illustrate  the  relation  between  Bz  and  the  direction  of  plasma  flow 
in  the  polar  cap  we  briefly  discuss  the  data  shown  in  Figures  3  to  12, 
covering  a  number  of  days  in  March,  May  and  June  1987,  and  in 

February  1988.  The  three  days  of  drift  data  in  Figure  3  show  a 
typical  behavior;  the  drift  direction  is  largely  antisunward  (asw), 
with  large  deviations  from  asw  when  Bz  turns  positive. 
Determination  of  the  response  time  of  the  ionospheric  reaction  to 

sign  changes  in  Bz  is  only  possible  in  some  clearly  defined  cases.  At 
00  UT  on  March  4,  it  takes  at  least  30  min  (there  is  a  gap  in  the  Bz 
data)  for  the  plasma  flow  to  become  asw  after  Bz  turns  negative.  On 
March  5  at  9:30  UT,  the  response  time  is  about  2  hours. 

On  March  16  and  17  (Figure  4),  Bz  is  predominantly  positive  and  the 
plasma  flow  is  very  irregular.  At  8:00  and  9:30  UT  on  March  16 

(Figure  5),  Bz  goes  through  zero,  and  the  drift  direction  changes 

within  about  30  min.  On  March  18,  Bz  is  predominantly  positive  and 
the  drift  direction  is  approximately  sunward.  At  20:30  UT,  Bz  turns 
negative  and  the  drift  directions  becomes  asw  almost  immediately. 
On  March  22,  Bz  turns  negative  atl6:30  UT  (Figure  6),  and  it  takes 
about  two  hours  before  the  drift  turns  asw. 


The  May  data  in  Figure  7  show  an  excepional  situation.  At  18:30  UT 
on  May  20,  Bz  turns  negative  and  remains  negative  for  9  hours, 
however,  the  drift  does  not  turn  a^w.  During  the  entire  three  days 
the  drift  velocity  is  unusually  small,  generally  below  200  m/s.  The 
June  data  (Figure  8)  show  a  fairly  regular  asw  drift;  unfortunately, 
the  Bz  data  have  large  gaps.  An  ionospheric  response  time  of  about 
one  hour  is  seen  at  19  UT  on  June  18,  when  Bz  turns  positive. 


Figure  9  (Feburary  9,  1988)  shows  another  example  at  21  UT  where 
Bz  is  clearly  negative,  but  the  drift  deviates  significantly  from  asw. 
A  similar  behavior  is  seen  for  the  first  90  min  on  February  25 
(Figurel2).  On  February  10  (Figure  10)  the  drift  turns  sunward  at 
16:45  UT,  about  30  min  after  Bz  turns  positve.  A  very  clear  Bz 
control  can  be  seen  for  February  23:  Bz  oscillates  through  zero 
several  times  (Figure  11)  and  the  drift  changes  accordingly  form  asw 
to  sunward  and  back  to  asw  with  a  response  time  of  less  than  15 
min. 
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Conclusions 


For  the  first  time  it  has  been  possible  to  conduct  extended  ground- 
based  measurements  of  the  polar  cap  F  region  plasma  flow  and 
compare  the  response  of  the  drift  to  the  variations  in  Bz.  It  has  been 
verified  that  the  plasma  flows  asw  when  Bz  is  negative,  and 
generally  sunward  when  Bz  is  positive.  There  are  exceptions, 
however,  where  the  drift  is  not  asw  for  some  4  hours,  in  one  single 
case  even  9  hours,  even  though  Bz  is  negative.  But  these  events  are 
real  exceptions.  The  response  time  of  the  plasma  flow  to  sign 
changes  of  Bz  varies  from  about  5  min  to  2  hours.  This  explains  why 
there  are  so  widely  differing  statements  in  the  literature  about  the 
response  time.  It  was  originally-  thought  /8/  that  the  response  time 
is  2  hours  or  more,  while  EISCAT  observations  suggested  a  5  min  or 
faster  response  /9/. 
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ABSTRACT 

Results  of  ionospheric  drift  measurements  with  a  Digisonde  256  digital  ionospheric 
sounder  located  at  Qaanaaq,  Greenland  (87°N,  corrected  geomagnetic  latitude),  are 
presented.  Digisonde  drift  data  have  been  related  to  the  interplanetary  magnetic 
field  (IMF)  measured  by  the  IMPS  satellite  for  thirty-two  days  during  1986,  1987 
and  1988.  Extremely  good  statistical  agreement,  between  these  measurements  and 
convection  directions  derived  from  satellite  instrumentation  is  demonstrated  when 
the  IMF  z  component  is  negative.  For  Bz  negative  and  By  negative  the  daily 
average  convection  direction  is  centred  on  -12° (anti-clockwise)  from  the  anti- 
sunward  direction.  When  Bz  is  negative  but  By  is  positive  the  convection  direction 
is  centred  on  +36°. These  directions  differ  by  respectively  3®  and  9°  from  the 
models  of  Heppner  and  Maynard  (1987).  The  variation  about  these  values  is  ~  ± 
20°.The  excellent  agreement  between  the  Digisonde  measurements  and  models 
derived  from  satellite  measurements  demonstrates  the  utility  of  the  Digisonde  for 
making  ground  based  convection  measurements  in  the  polar  cap  F  region  when  Bz 
is  south.  The  convection  directions  under  conditions  of  positive  Bz  have  also  been 
examined  and  we  have  measured  three  types  of  temporal  variation  in  azimuth; 
namely  an  ordered  and  slowly  (OS)  varying  change  in  direction,  an  ordered  and 
quickly  (OQ)  varying  change  in  direction,  and  disordered  (D)  variations  in 
direction.  The  latter  are  believed  to  result  from  a  breakdown  of  the  analysis 
technique  due  to  velocity  shears  in  the  vicinity  of  polar  cap  arcs  and  we  estimate 
that  they  accou,  ‘  for  ~259o  of  the  measurements  when  Bz>0.  When  Bz  is  positive 
and  By  is  negative  our  small  subset  of  OS  measurements  supports  the  distorted  two 
cell  model  of  Heppner  and  Maynard  (1987).  The  remainder  of  the  measurements 
show  no  well  defined  daily  average  convection  direction  or  diurnal  variation. 
Likewise  for  Bz  positive  and  By  positive  no  well  defined  convection  direction  can 
be  discerned,  nor  can  any  diurnal  variation.  The  existence  of  OQ  variations  when 
Bz  is  >0  suggests  that  meaningful  average  statistical  convection  patterns  may  be 
much  harder  to  synthetize  than  similar  patterns  when  Bz<0. 
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1.  INTRODUCTION 

Magnetospheric  electric  fields  play  an  important  part  in  the  control  of  the  high 
latitude  ionosphere.  Ionospheric  convection  due  to  magneto-hydrodynamic 
coupling  between  the  solar  wind  and  the  geomagnetic  field  was  first  suggested  by 
Dungey  (1961)  who  proposed  a  two  cell  pattern  with  antisun  ward  flow  across  the 
polar  cap  when  the  interplanetary  magnetic  field  (IMF)  shows  a  southward  directed 
component,  (see  also  Axford  and  Hines,  1961).  This  simple  two  cell  convection 
pattern  has  been  further  refined,  extended  and  incorporated  in  models  of  the  high 
latitude  ionosphere  (eg  Sojka  et  al,1986,  Sojka  and  Schunk,1987,  Quegan  et 
al.,1982).  With  this  progress  has  come  an  improved  understanding  of  the  gross, 
and  some  of  the  detailed,  aspects  of  the  high  latitude  ionosphere. 

Experimental  information  to  act  as  the  basis  of  the  convection  models  and  as  a  seed 
for  the  time  dependent  high  latitude  ionospheric  models  stems  from  a  wide  variety 
of  sources  (Sojka  and  Schunk,  1986).  Satellite  convection  measurements  have  been 
made  using  double  probe  instmmentation  aboard  0G06  (eg  Heppner,  1972, 1977), 
instmmentation  on  board  Atmosphere  Explorer  C  (AE-C)  (eg  Heelis  et  al.,  1976) 
and  using  data  from  Dynamics  Explorer  2  (DE-2)  (eg  Heppner  and  Maynard, 
1987).  Other  in-situ  measurements  of  the  electric  field  have  been  made  by  balloon 
probes  (Mozer  and  Lucht,1974)  and  rocket  measurements  (eg  Maynard  and 
Johnstone,1974).  Various  ground  bas'.d  remote  sensing  techniques  have  also  been 
employed  including  Incoherent  Scatter  Radars  (ISR)  located  at  Chatanika  (Foster, 
1981;1983),  Millstone  Hill  (Evans  et  al,  1980)  and  in  Scandinavia  using  EISCAT 
(Folkestad  et  al,1983).  Coherent  HF  and  VHF  Radars  have  provided  additional 
tools  to  investigate  both  E-region  irregularities  (Greenwald,1978;  Jones  et  al,1985) 
and  F-region  irregularities  (Greenwald  et  all 985;  Ruohoniemi  et  al,  1987)  and 
their  associated  electric  fields. 

Whilst  these  data  have  been  used  as  input  to  theoretical  models  to  deduce  the 
temporal  evolution  of  the  ionosphere,  they  have  been  more  usually  employed, 
either  separately  or  together,  to  derive  average  high  latitude  ionospheric 
convection  patterns.  When  the  IMF  z  (northward)  component  is  negative  many  of 
these  patterns  exhibit  the  basically  similar  form  of  a  two  cell  convection  pattern. 
These  include  those  derived  by  Heppner  (1977)  and  later  by  Heppner  and  Maynard 
(1987)  using  pattern  recognition  of  the  cross  polar  cap  potential  distribution  in 
order  to  minimise  the  spatial  blurring  which  must  necessarily  be  introduced  by 
averaging.  They  also  include  the  mathematically  convenient  models  of  Heelis  et  al. 
(1982)  which  describe  the  major  features  of  the  two  cell  convection  pattern  without 
undue  complication.  While  the  two  cell  models  of  these  and  other  authors  (see 
Schunk  (1988)  and  references  therein)  concentrated  on  periods  when  the  IMF  z 
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component  is  southwards  some  models  have  examined  the  convection  pattern  for 
Bz  north.  Burke  et  al.  (1979)  suggested  a  four  cell  convection  pattern  and  Lassen 
(1979)  interpreted  the  polar  cap  auroral  arcs,  observed  during  magnetically  quiet 
periods,  as  being  the  result  of  a  shear  associated  with  a  third  central  polar  cap  cell. 
Potemra  et  al.  (1984)  also  proposed  three  and  four  cell  convection  patterns  with  the 
number  of  cells  and  their  size  and  location,  dependent  on  the  sign  of  By.  Some  of 
the  expected  ionospheric  features  of  these  multi-cell  patterns  have  been  modeled  by 
Sojka  and  Schunk  (1987).  In  contrast  to  the  suggestion  that  three  and  four  cell 
patterns  represent  the  average  convection  pattern  when  Bz  is  northward  Heppner 
and  Maynard  (1987)  argued  in  favour  of  a  distorted  two  cell  model.  The  distortion 
proposed  was  that  of  a  translational  stretching  of  the  evening  cell  focus  towards  and 
beyond  noon,  accompanied  by  a  rotational  twist  of  the  streamlines  around  the  foci. 
Progressively  greater  distortion  was  proposed  as  Bz  increased. 

Each  of  the  different  instruments,  used  in  the  derivation  of  the  models,  has 
advantages  and  disadvantages.  Satellite-borne  instrumentation  provides  good  spatial 
coverage  of  the  high  latitude  region  at  the  expense  of  a  long  return  time  to  any 
particular  geographic  region  -  typically  the  orbital  period  of  a  low  altitude  satellite 
is  ~90  minutes.  In  addition,  any  dedu.  'ons  regarding  the  cross  polar  cap  potential 
from  satellite  data  must  include  some  temporal  blurring  because  the  convection 
field  is  assumed  to  remain  stable  during  the  period  (>10  min)  that  the  satellite 
crosses  the  polar  cap.  Multistatic  incoherent  scatter  radars  provide  excellent 
temporal  resolution  at  the  expense  of  some  limitation  on  their  spatial  coverage. 
Monostatic  ISR  facilities  are  limited  to  making  line  of  sight  velocity  measurements 
unless  beam  swinging  techniques,  coupled  with  assumptions  of  spatial  uniformity 
and  temporal  stability  of  the  bulk  plasma  velocity,  are  made.  Other  practical 
constraints  also  severely  limit  the  operational  availability  of  ISR  facilities,  making 
synoptic  measurements  difficult.  Tbe  coherent  HF  and  VHF  radars  provide  good 
temporal  resolution,  their  spatial  coverage  is  significant,  (particularly  HF  radars), 
and  their  near  continuous  operation  makes  them  ideal  for  routine  studies.  The 
primary  disadvantage  of  coherent  radars  is  that  the  necessary  radar  wave  vector  to 
geomagnetic  field  orthogonality  condition  (Bates  and  Aroee,1970)  cannot  always  be 
met  in  the  desired  region  of  the  ionosphere.  Further  coverage  restrictions  are 
imposed  by  the  logistics  of  locating  the  radar.  These  constraints  are  particularly 
severe  for  the  high  latitude  VHF  radars  which  are  unable  to  benefit  from 
ionospheric  refraction  to  bring  the  ray  path,  at  the  interaction  height,  orthogonal  to 
the  local  geomagnetic  field.  Nielsen  and  Schlegel  (1985)  have  also  shown  that  for 
metre  scale  irregularities  in  the  E-region  correction  factors  must  be  applied  to 
accurately  determine  the  drift  magnitude  and  thus  electric  field.  In  the  F-region  this 
does  not,  however,  appear  to  be  a  problem  (Ruohoniemi  et  al,  1987). 

An  alternative  and  complementary  form  of  instrumentation  for  ionospheric  drift 
measurements  is  a  digital,  ground-based  ionosonde  employing  an  antenna  array  and 
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phase  measurements  to  determine  the  arrival  angles  of  the  reflected  signals.  In  this 
study  a  Digisonde  256,  designed  and  built  at  the  University  of  Lowell  Center  for 
Atmospheric  Research  (ULCAR)  is  employed.  Like  the  coherent  radars  the 
Digisonde  is  ideal  for  synoptic  measurements  but  unlike  the  coherent  radars  it  does 
not  rely  on  orthogonality  between  the  ray  direction  and  the  geomagnetic  field. 
Consequently,  the  logistical  constraints  of  location  are  mitigated.  The  Digisonde 
can  determine  the  drift  by  measuring  the  motion  of  isoionic  contours.  This  is  the 
most  likely  condition  at  mid-latitudes  where  the  natural  unevenness  of  these 
surfaces,  accentuated  by  the  transit  of  gravity  waves,  provides  ample  regions  where 
off  vertical  perpendicularity  is  achieved.  In  the  high  latimde  ionosphere,  v/here  the 
measurements  described  here  were  made,  field  aligned  irregularities  (FAIs)  also 
exist  throughout  the  ionosphere  (spread-F).  Here  the  orthogonality  condition  can 
often  be  met  due  to  refraction,  by  the  background  ionisation,  of  the  rather  low 
frequency  waves  used  in  these  measurements  coupled  with  the  large  range  and 
altitude  extent  of  FAIs  (Buchau  and  McNamara,  1986).  The  Digisonde  also  yields 
good  time  resolution  with  a  nominal  value  of  five  minutes  -  under  certain 
favourable  conditions  this  may  be  improved  upon.  The  disadvantage  of  the 
Digisonde  is  that  it  illuminates  a  relatively  small  area  of  the  ionosphere  above  the 
antenna  and  current  data  reduction  techniques  require  uniform  ionospheric  drift 
within  tlie  antenna  beam. 

Reinisch  et  al.  (1987)  showed  good  agreement  between  Digisonde  measurements  of 
F-region  plasma  drift  and  the  major  features  of  plasma  convection  patterns. 
Buchau  et  al.  (1988),  by  examining  various  specific  cases,  showed  that  not  only 
were  the  major  convection  features  obser\'ed  but  also  details  such  as  convection 
reversals  in  response  to  quiet  magnetic  conditions.  These  two  papers  suggested  that 
continuous  monitoring  of  the  large-scale  features  of  the  convection  pattern,  using  a 
spatially  diverse  network  of  relatively  inexpensive  Digisondes,  is  possible. 

The  determination  of  drift  velocities  using  the  Digisonde  is  subject  to  certain 
assumptions  which  are  not  universally  applicable.  As  a  first  step  to  investigating  the 
ability  of  the  technique  for  making  ionospheric  convection  measurements  in  the 
high  latitude  region  this  paper  discusses  Digisonde  drift  measurements,  from  a  site 
close  to  the  corrected  geomagnetic  pole,  as  a  function  of  the  IMF.  At  this  location, 
under  Bz<0  conditions, the  plasma  convection  is  known  to  be  well  ordered  and 
quasi-laminar  over  large  areas  producing  conditiono  for  which  the  currently 
available  Digisonde  drift  analysis  technique  has  been  developed.  We  demonstrate 
the  efficacy  of  the  technique  for  polar  cap  convection  measurements  by  comparing 
our  data  with  '^ther  measurements.  Specifically,  our  statistical  study  will  be 
compared  to  the  extensive  Bz  negative  statistical  models  of  Heppner  and  Maynard 
(1987).  A  number  of  previous  studies  have  demonstrated  broadly  anti-sunward 
flow  during  periods  of  negative  Bz.  The  models  of  Heppner  and  Maynard  1987j 
show  that  when  By  is  negative  the  convection  streamlines,  near  the  northern 
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geomagnetic  pole,  deviate  in  an  anti-clockwise  sense  from  the  anti-sunward. 
Conversely,  during  periods  when  By  is  positive  the  average  convection  direction  is 
rotated  clockwise.  We  also  discuss  our  measurements  of  ionospheric  drift  when  Bz 
is  positive.  A  statistical  validation  of  the  Digisonde  drift  measurement  technique 
for  plasma  convection  studies  when  Bz  is  northward  is,  however,  difficult  since 
there  are  no  generally  acceptable  models  of  the  plasma  convection  pattern  . 

Whilst  the  Digisonde  is  capable  of  measuring  drift  speeds,  as  well  as  directions,  this 
paper  will  deal  only  with  a  comparison  of  the  drift  directions.  The  paper  is  divided 
into  seven  sections.  In  Section  2  a  description  of  the  Digisonde  is  given,  together 
with  a  discussion  of  the  assumptions  built  into  the  measurement  technique,  and  in 
Section  3  the  temporal  evolution  of  the  drift  direction  is  examined  by  example. 
Section  4  describes  the  method  of  selecting  Digisonde  drift  and  IMF  data  with  a 
view  to  examining  the  sensitivity  of  the  measurement  technique  to  variations  in  the 
IMF  components  By  and  Bz.  Section  5  presents  the  results  of  these  analyses  and 
compares  the  results  to  the  satellite  convection  models.  Section  6  deals  with  the 
problem  of  how  to  interpret  data  during  Bz  positive  periods  which  is  more  difficult 
to  do  than  during  periods  when  Bz  is  negative.  In  Section  7  we  present  the 
conclusions  of  the  paper. 


2.  DESCRIPTION  OF  THE  DRIFT  MEASUREMENT  TECHNIQUE 
AND  HARDWARE  AT  QAANAAQ 

The  Digisonde  (Bibl  and  Reinisch,  1978)  drift  technique  uses  spectral  analysis  of 
ionospheric  echoes  to  locate  the  source  of  reflected  signals  (Pfister,  1974a,  b);  this 
is  an  important  aspect  of  the  technique  which  obviates  the  need  for  a  large  antenna 
array.  \^ilst  the  Digisonde  drift  technique  can  be  employed  to  measure  E-region 
drifts  this  paper  is  restricted  to  F-region  measurements. 

The  hardware  consists  of  a  Digisonde  256  digital  ionospheric  sounder  with  one 
transmitter  antenna  and  seven  receiver  antennas. The  transmitter  antenna  has  a 
beamwidth  of  >90°  and  illuminates  a  large  area  of  the  F-region  above  the 
ionosonde.  The  transmitter  delivers  1 0  kW  peak  pulse  power. 

Seven  receiver  antennas  are  deployed  on  a  100  m  equilateral  triangle  (Figure  1) 
although  sometimes  only  the  central  four  antennas  are  used.  The  array  receives 
echo  signals  from  a  frequency  dependent  area  which  is  <  9x10'^  km^  in  the  F- 
region  (equivalent  to  (-3°)-  in  latitude  degrees).  Specifically,  at  10  MHz  the  6  dB 
beamwidth  of  the  array  is  -15°  with  correspondingly  larger  values  at  the  lower 
frequencies.  Echoes  are  obtained  from  all  points  on  a  F-region  isoionic  surface 
where  the  plasma  frequency  is  equal  to  the  sounding  frequency,  and  where  trie 
vector  normal  to  the  surface  points  towards  the  sounder  (neglecting  refraction 
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effects).  In  tlie  polar  cap  F-region  ionosphere  kilometre  scale  irregularities  and 
medium  scale  polar  cap  patches  (eg  Buchau  et  ah,  1985;  Weber  et  al.,1986)  create 
a  number  of  points  off  the  vertical,  where  these  conditions  are  met.  The  ionospheric 
drift  is  determined  by  measuring  the  Doppler  frequencies  of  the  individual  echoes. 

In  the  presence  of  multiple  echoes  a  three  dimensional  drift  velocity  vector  can  be 
derived  from  this  small  array  by  using  the  radial  Doppler  velocity  information 
associated  with  each  echo.  This  may  be  achieved  subject  to  two  conditions.  Firstly, 
the  drift  velocity  must  be  uniform  over  the  illuminated  area.  This  is  a  particularly 
important  restriction  which  we  shall  demonstrate  is  met  in  the  polar  cap  under  some 
IMF  conditions.  Secondly,  the  technique  assumes  that  refraction  effects  can  be 
ignored.  The  validity  of  this  assumption  has  been  demonstrated  by  Bullett  (1990). 

The  Qaanaaq  Digisonde  records  vertical  ionograms  every  5  or  15  minutes  and  the 
interval  between  each  ionogram  is  utilised  to  make  multiple  measurements  of  the 
ionosphere  bulk  velocity.  For  these  velocity  measurements  the  single  ionosonde 
receiver  is  multiplexed  to  either  all  seven  or  just  four  of  the  antennas  at  the  pulse 
repetition  frequency.  The  velocity  measurements  are  made  at  two  frequencies  and 
two  virtual  ranges.  The  operating  frequencies  and  range  gates  are  automatically 
determined  by  a  real  time  software  system  known  as  ARTIST  (Automatic  Real 
Time  Ionogram  Scaler  with  True  height  analysis).  ARTIST  is  used  to  scale  the 
ionogram  immediately  prior  to  the  drift  measurement  and  it  selects  two  lower  F- 
region  frequencies  (~0.5foF2)  and  two  range  gates  which  correspond  to  strong 
signals.  Quadrature  sampling  of  the  signal  received  at  each  antenna  is  performed. 
These  signals  are  Fourier  transfomied,  in  real  time,  to  yield  a  complex  Doppler 
spectmm  for  each  antenna  signal,  (Ai,0i)j,  with  i=l,...64  (64  spectral  lines)  and  j 
=1,...4  or  7  (for  four  or  seven  receiving  antennas),  which  are  recorded  on  magnetic 
tape  for  post  processing. 

Dozois  (1983)  has  described  how  the  echo  data  are  being  processed  to  generate  the 
required  drift  information.  In  the  post  processing  phase  those  Doppler  lines  which 
have  significant  signal  energy  are  identified  as  "sources"  of  reflected  signal.The 
spectral  indices,  i  give  the  Doppler  frequency,  ds  of  each  source: 

dj  =^v.ks  s  =  1,2,.. ..(source  index) 

where  v  is  the  source  drift  velocity  and  k.s  the  vector  from  the  reflection  point,  s  to 
the  receiving  antennas.  From  the  spectral  phases  0ij  across  tlie  antenna  array,  the 
source  wave  vector,  ks  is  determined.  The  intersection  of  tliese  ks  vectors  with  tlie 
horizontal  plane,  at  the  virtual  range  gate,  generates  skymaps  of  the  source 
locations.  Figure  2  shows  a  skymap  sequence  covering  one  and  a  half  hours. 
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The  three-dimensional  velocity  vector  v  is  estimated  from  a  collection  of  source 
locations,  defined  by  ks  and  Doppler  shifts  ds,  by  minimizing  the  least-square 
error  : 


[d^  -  —  v.k  s  =  1,2,3 . (source  index). 

s  ^  71  ^ 

The  teclmique  breaks  down  when  there  are  less  than  three  sources  and  the  errors 
become  large  when  the  sources  are  tightly  clustered. 

The  unaliased  upper  limit  of  measurable  drift  velocities  is  constrained  by  the 
antenna  sampling  frequency.  To  maximize  this  limit  more  recent  Qaanaaq 
Digisonde  drift  measurements  use  only  four  antennas  and  a  200  Hz  pulse  repetition 
frequency  together  with  tw-o  frequencies  and  two  ranges.  Earlier  measurements 
used  a  seven  antenna  configuration  which  decreased  the  sampling  frequency.  A  ±32 
point  spectrum  is  generated  every  5.12  s  with  a  spectral  resolution  of  0.1953  Hz. 
In  the  four  antenna  configuration  the  frequency  range  for  unaliased  signals  is 
±6.25  Hz.  This  frequency  range  corresponds  to  a  line  of  sight  velocity  of  312.5  ms- 
1  at  3.0  MHz,  a  typical  measurement  frequency.  Most  sources  are,  however,  found 
nearer  the  zenith  than  tlie  horizon  and  consequently  the  maximum  horizontal 
velocity  is  significantly  greater  than  312  ms'k  For  a  source  located  at  19°  to  the 
zenith  the  equivalent  maximum  horizontal  velocity  is  1000  ms-^  with  a 
corresponding  smaller  cone  appropriate  to  higher  velocities.  During  each  5  minute 
period  40  measurements  were  made,  10  at  each  frequency  and  range.  Tlie  combined 
sources  (typically  ten  to  several  hundred)  from  each  frequency  were  combined, 
assuming  no  variation  in  the  bulk  velocity  over  five  minutes  and  no  N’ariation  in 
height,  to  yield  two  ionospheric  velocity  vectors. 

3.  DRIFT  DIRECTIONS 

F-region  ionospheric  drift  velocity  vectors,  recorded  by  a  Digisonde  256  located  at 
Qaanaaq,  Greenland  (77. 5N,  290. 8E  geographic;  86. 8N  corrected  geomagnetic 
latitude)  have  been  determined  for  the  days  given  in  Table  1.  The  period  studied 
covers  data  from  a  period  of  about  two  years^and  includes  1986,  1987  and  1988. 
Interplanetary  magnetic  field  data  from  the  IMP  8  satellite  have  been  employed  for 
comparison  purposes.  These  days  w'ere  chosen  from  a  larger  data  set  such  that  both 
the  IMP8  satellite  was  correctly  placed  (.see  Section  4)  and  such  that  5  minute 
resolution  ionospheric  drift  vectors  were  available.The  corrected  geomagnetic 
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coordinate  system  (Hultqvist  (1958a,]  958b)  has  been  used  to  reference  directions 
using  the  coordinates  computed  by  Gustafsson  (1970). 

As  an  example  of  the  ionospheric  drift  measurements  Figure  3  presents  the  drift 
directions,  in  the  corrected  geomagnetic  coordinate  system,  for  the  days  040-042 
(9-11  Februar>')  inclusive  during  1988.  Each  point  is  the  median  drift  direction 
vector  calculated  from  measurements  within  frequency  bins  0.5  MHz  wide  during  a 
5  minute  period.  For  approximately  40%  of  the  5  minute  intervals  two  median  drift 
directions  were  obtained  because  the  automatically  selected  operating  frequencies 
fell  in  different  frequency  bins.  The  diagonal  lines  denote  the  anti-sunward 
direction  along  the  corrected  geomagnetic  noon-midnight  meridian.  On  all  three 
days  Bz  (Figure  4)  showed  variability  including  periods  of  both  significant 
negative  and  positive  values  (5-10  nT).  Day  040  exhibited  long  periods  when  By 
was  negative  (Figure  5)  while  conversely.  By  was  positive  through  much  of  days 
041  and  042. 

As  expected  from  the  results  of  numerous  other  experiments  a  comparison  of  the 
drift  directions  (Figure  3)  and  Bz  (Figure  4)  reveals  that  when  Bz  is  predominantly 
negative  there  is  a  strong  antisunward  trend  in  the  ionospheric  drift  direction,  see 
day  040.  Wlien  Bz  is  south  the  measured  drift  directions  also  vary  in  a  systematic 
way  from  one  measurement  to  the  next.  When  Bz  swings  northward  the  drift 
direction  variation  is  sometimes  well  ordered  but  sometimes  badly  ordered  and 
there  are  also  long  periods  when  the  drift  direction  departs  significantly  from  the 
anti-sunward.  At  2130  UT  on  day  041  the  Bz  component  can  be  identified  as 
directed  northwards  and  except  for  a  brief  interlude  between  00  and  01  UT  on  day 
042  it  remains  so  until  ~  0630  UT.  Throughout  the  period  00  to  0630  UT  on  day 
042  the  drift  directions  are  widely  distributed  over  about  90°  in  azimuth  and 
consistent  ordering  of  the  drift  direction  is  not  easily  discemable.  This  should  be 
contrasted  with  the  period  between  16  UT  and  2330  UT  on  day  042  when  Bz 
oscillates  about  zero  remaining  either  positive  or  negative  for  up  to  2  hours.  When 
Bz  is  positive  during  this  period  the  drift  direction  variation  is  well  ordered  (and 
quasi-periodic).  Here  also  the  drift  direction  deviates  significantly  from  the  anti¬ 
sunward. 

The  systematic  or  ordered  variation  in  drift  direction  from  one  five  minute 
estimate  to  the  next  is  an  important  characteristic  of  nearly  all  of  the  data 
corresponding  to  Bz  negative  conditions  together  with  much  of  the  data 
corresponding  to  Bz  positive  conditions 

During  the  period  between  16  and  23  UT  (Day  042)  visual  inspection  of  Figures  3, 
4  and  5  reveals  a  well  correlated  and  systematic  variation  of  tlie  drift  direction  in 
response  to  the  changes  in  the  IMF.  There  is  also  some  indication  that,  during  the 
one  to  two  hour  excursions  in  Bz,  the  drift  direction  fails  to  establish  a  quiescent 
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direction  applicable  to  that  Bz  polarity.  The  delay  of  the  ionospheric  signature  to 
the  IMF  change  is  short  in  this  example. 


4.  COMPARISON  OF  DRIFT  WITH  IMF  -TECHNIQUE 

Vo  ..'..Lermine  the  relationship  between  the  IMF,  measured  in  geocentric  solar 
magj-*..(.ospheric  coordinates  (GSM),  and  the  drift  direction  it  is  necessary  to  know 
the  Direction  of  the  IMF  when  the  drift  measurement  was  made.  IMF  data  from  the 
IMP  8  satellite  have  been  used  for  this  purpose.  IMP  8  is  in  an  ecliptic  plane  orbit 
which  extends  ~40  Re  (earth  radii)  in  the  GSE  (geocentric  solar  ecliptic)  x- 
directlon  and  ~30  Re  in  the  GSE  y  direction.  In  these  smdies  it  is  assumed  that  the 
magnetopause  is  located  at  x=10  Re  and  only  IMP  8  data  were  considered  when  the 
satellite  position  was  located  at  x>10  Re.  Variations  in  the  IMF,  measured  onTMP 
8,  tlien  manifest  themselves  at  the  magnetopause  and  in  the  ionospheric  F-region 
some  time  later.  Like  Heppner  and  Maynard  (1987)  we  have  defined  a  window 
preceding  the  drift  measurement  during  which  tlie  IMF  must  remain  stable  if  we 
are  to  knovv  unequivocally  the  direction  of  the  IMF  at  the  time  of  the  measurement. 
By  stable  we  mean  that  each  IMF  component  (in  GSM  coordinates)  must  remain 
eitlier  positive  or  negative  during  the  whole  of  the  window.  This  window  may  be 
divided  into  two  parts.  The  first  is  the  time  delay  for  IMF  disturbances  to  propagate 
between  satellite  and  magnetopause.  The  second  is  tlie  time  delay  of  the  ionospheric 
response  due  to  a  variation  in  the  IMF  at  the  magnetopause.  The  length  of  the  first 
part  was  determined  by  considering  both  the  corotation  and  the  convective  time 
delay  approximation  (Couzens  and  King,  1986)  of  IMF  disturbances  between  the 
satellite  and  magnetopause  at  the  sub-solar  point.  In  order  to  estimate  this  delay  we 
have  assumed  that  the  solar  wind  speed  v/as  500  kms'h  that  the  earth's  speed  was  30 
kms'l  and  the  solar  rotation  speed  was  428  kms'k  The  maximum  corotation  delay 
for  satellite  positions  corresponding  to  x  >  10  Re  is  then  ~13  minutes,  the  maximum 
convective  delay  ~8  minutes  and  the  maximum  differential  delay  is  ~10  minutes. 
Russell  et  al.  (1980)  have  shown  that  neither  approximation  can  generally  be 
considered  better  than  tlie  other  and,  therefore,  we  have  set  die  first  component  of 
the  window  to  be  15  minutes  in  order  to  include  the  worst  case  propagation  delay 
from  satellite  to  magnetopause.  With  regard  to  the  second  part  of  the  window  Todd 
et  al.(1988)  and  Etemadi  et  al.  (1988)  have  demonstrated  that  the  ionospheric 
response  time  is  >  2  minutes  but  <15  minutes.  The  total  window  requiring  stable 
IMF  components  was,  therefore,  set  at  30  minutes. 

GSM  coordinates  rather  than  the  GSE  coordinates  were  selected  because  IMF 
statistical  correlations  with  other  geophysical  parameters  such  as  Kp  and  AE  are 
better  in  the  GSM  coordinate  system  (eg  Burch, 1974).  This  is,  however,  in  conflict 
with  Heppner  and  Maynard  (1987)  who  used  die  GSE  coordinate  system.  Another 
difference  between  our  analysis  procedure  and  that  of  Heppner  and  Ma^mard 
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(1987)  is  that  their  model  is  applicable  to  a  range  of  Kp  values  between  3+  and  4*. 
We  explored  filtering  our  data  according  to  these  criteria  but  were  unable  to  detect 
any  statistically  significant  variations  in  average  drift  direction.  All  the  results  in 
this  paper,  therefore,  reflect  the  full  range  of  Kp  values  applicable  to  the  data  set. 

5.  COMPARISON  OF  DRIFT  DIRECTIONS  WITH  IMF  -  RESULTS 

Tlie  aforementioned  criteria  have  been  applied  to  the  thirty-two  days  of  data  to 
determine  how  the  drift  direction  varies  when  Bz  is  negative  and  By  is  negative  and 
also  how  it  varies  when  Bz  is  negative  but  By  is  positive.  No  account  has  been  taken 
of  the  Bx  condition  except  for  tlie  initial  requirement  of  stability.  A  similar  analysis 
has  been  carried  out  when  Bz  is  positive.  Drift  direction  data  have  been  binned  in 
4°  direction  intervals.  In  order  to  ensure  that  the  magnetosphere  is  responding 
strongly  to  the  various  IMF  conditions  an  additional  requirement  of  a  1  nT  guard 
band  has  been  imposed.  That  is  any  drift  vector  which  corresponds  to  a  30  minute 
average  IMF  component  such  that  IByl  or  IBzl  <  1  nT  is  ignored.  Tne  number  of 
points  selected  from  each  day  are  given  in  Figure  6.  No  attempt  was  made  to  select 
periods  containing  specific  IMF  conditions  and  there  are  relatively  few  drift  data 
when  Bz  is  positive  and  By  is  negative. 

Bz<0.,  By<0 

Four  hundred  and  twelve  drift  vectors  were  selected  meeting  these  various  criteria 
and  they  are  well  distributed  throughout  the  period  13  to  04  UT  but  with  relatively 
few  points  at  other  times  (Figure  7,  top  panel).  Their  deviation  from  the  anti- 
sunward  direction  (Figure  7,  bottom  panel)  shows  the  convection  direction  centred 
on  -r2‘’ (clockwise  rotation)  from  the  anti-sunward.  Simple  inspection  reveals  a 
broad  symmetrical  spread  of  convection  directions  with  the  probability  still  at  50% 
of  its  hiaximum  value  at  ~  ±20°  of  the  modal  drift  direction. 

Heppher  and  Maynard  (1987)  using  satellite  data  presented  two  models  for  Bz<0 
and  By<0,  corresponding  to  two  polar  cap  potential  distributions,  with  average 
convection  directions  at  Qaanaaq  of  0°  for  the  A  model  and  -15°  for  tlie  DE  model. 
Maynard  (private  communication,! 989  )  has  indicated,  however,  that  the  A  model 
is  transitional  between  the  DE  and  BC  models  for  negative  By  and  that  its 
occurrence  probability  is  less  than  the  DE  model.  As  a  consequence,  it  may  not  be 
surprising  that  the  Digisonde  measurements  do  not  reveal  a  secondary  peak  at  0°. 
The  ground-based  Digisonde  average  of  -12°  is  in  good  agreement  with  the  satellite 

pwororro  rvP  .1 

Associated  with  cac/i  individual  five  minute  resolution  drift  direction  estimate  is  a 
measurement  inaccuracy  (ie  error).  The  magnitude  of  this  error  is  controlled  by  a 
number  of  factors  including  the  number  and  spatial  distribution  of  reflection  points 
on  the  skymap,  the  temporal  stability  of  the  drift  direction  within  the  antenna  beam 
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and  the  spatial  uniformity  of  the  drift  within  the  antenna  beam.  An  examination  of 
the  error  magnitudes  consequently  provides  an  indication  of  the  validity  of  the 
Digisonde  drift  measurement  technique,  at  least  using  the  current  data  processing 
algorithm,  for  F-region  polar  cap  drift  measurements.  The  difference  between  the 
median  and  quartile  drift  direction  estimates  was  used  as  a  measure  of  this  error. 
For  each  of  the  drift  directions  meeting  the  IMF  selection  criteria  the  difference 
between  the  median  and  both  the  upper  and  lower  quartile  drift  direction  estimate 
was  evaluated.  Figure  8a  shows  the  mean  of  these  two  differences  as  a  density  plot. 
The  large  majority  of  the  errors  (75  %)  are  less  than  5°.  This,  in  association  with 
the  ordered  variation  of  drift  direction  usually  seen  during  these  IMF  conditions,  is 
indicative  of  good  data  and  processing  reliability.  Consequently,  the  distribution  of 
drift  directions  about  the  mode  (Figure  7)  appears  to  be  principally  due  to  the 
physical  variability  of  drift  directions  and  does  not  appear  to  be  an  artifact  of  the 
equipment  and  data  processing. 


Bz<0,  By>0 

Four  hundred  and  three  drift  vectors  met  these  criteria  and  their  time  distribution  is 
weir  spread  over  the  entire  day  (Figure  9,  top  panel).  The  deviation  from  the  anti- 
sunward  direction  peaks  at  24°  although  a  smooth  cur\'e  is  centred  around  36°. 
Again  there  is  a  broad  spread  in  convection  directions. The  Heppner  and  Maynard 
(1987)  model  places  the  average  convection  direction  at  Qaanaaq  under  these  IMF 
conditions  at  45“  which  is  in  good  agreement  witli  tlie  Digisonde  measurement  of 
36°. 

The  individual  drift  direction  errors  associated  with  tliis  IMF  configuration  (Figure 
8b)  lead  to  identical  conclusions  to  those  drawn  above  for  Bz<0  and  By<0.  The 
number  of  points  with  errors  less  than  5°  is  72%  which  is  very  close  to  the  value 
obtained  for  the  Bz<0,  By<0  condition. 

Bz>0,  By<0 

A  statistical  examination  of  Digisonde  data  under  positive  Bz  conditions  cannot 
hope  to  validate  the  measurement  technique,  because  there  is  no  commonly  held 
view  on  the  average  pattern.  These  data  have,  however,  been  examined  in  order  to 
provide  insight  into  tlie  Digisonde  measurement  results  for  Bz  nortli. 

Figure  10  presents  all  the  data  when  Bz  was  positive  and  By  negative. 
Unfortunately,  only  175  points  were  available;  but  surprisingly  for  this  limited  data 
set  the  average  convection  direction  (bottom  panel)  is  still  relatively  well  defined 
and  centred  around  -16°  from  the  anti-sunward  direction.  Sixty  eight  of  the  points 
selected  for  inclusion  in  this  figure  were,  however,  derived  from  a  relatively  .short 
period,  between  0230  and  05  UT  (approximately  2330  MLT  to  02  MLT),  on  day 
344, 1988  when  the  convection  direction  remained  close  to  the  anti-sunward.  This 
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period  is  contained  within  a  longer  time  frame  from  0130  to  07  UT  during  which 
Bz  was  positive.  (We  are  unfortunately  unable  to  comment  on  the  period 
immediately  after  07  UT  because  IMPS  data  were  not  available.)  This  Bz  positive 
frame  was  preceded  by  a  Bz  negative  frame  lasting  from  1930  UT  on  the  previous 
day.  'Fhe  latter  frame  was  characterised  by  strong  antisunward  convection  and  this 
continued  into  the  Bz  positive  period  until  05  UT.  A  characteristic  of  the  anti¬ 
sunward  drift  is  well  ordered  data  with  a  change  of  only  a  few  degrees  between 
successive  five  minute  estimates  of  tlie  drift  direction.  After  05  UT  the  measured 
drift  directions  became  highly  variable  and  disordered.  Likewise  another  34  points 
which  exhibit  a  slowly  varying  and  ordered  drift  direction  stem  from  the  period 
between  20  UT  and  2330  UT  (approximately  17  MLT  to  2030  MLT)  on  day  294, 
1987.  These  two  days  of  data  severely  bias  the  graph  and  without  them  no  well 
defined  average  drift  direction  would  be  discemable. 

The  slowly  varying  ordered  data  also  exhibit  small  errors  associated  with  the 
individual  drift  measurements  (Figure  8c)  but  the  remaining  selected  points  are 
associated  with  a  wider  distribution  of  errors  ranging  from  tlie  small  to  the  large. 
For  the  entire  category  of  data  when  Bz>0  and  By<0  only  48%  of  the  points 
exhibited  errors  less  than  5°  -  a  performance  some  25%  inferior  to  the  Bz  negative 
data  set. 


Heppner  and  Maynard  (1987)  argued  in  favour  of  a  distorted  two  cell  model  when 
Bz  was  positive.  This  distortion  places  Qaanaaq,  at  87°  CG  latitude,  well  into  the 
morning  cell  and  a  complicated  diurnal  variation  of  convection  direction  results 
ranging  from  sunward  through  anti-sunward.  The  ordered  data  presented  in  Figure 
10  are  clustered  in  the  period  20  -  2330  UT  and  0230  -  5  UT  which  corresponds  to 
convection  directions,  scaled  from  the  weak  Bz  positive  model  of  Heppner  and 
Maynard  (1987),  of  between  0°  and  -15°  and  -15°  respectively.  These  data  are, 
therefore,  not  in  disagreement  with  this  model;  they  are,  however,  in  disagreement 
with  the  various  three  and  four  cell  models  which  sugcest  anti-sunward  flow  durina 
positive  Bz  periods.  The  significance  of  such  a  small  data  subset  is,  however,  in 
question.  How  to  interpret  the  remaining  group  of  selected  data  points  showing  a 
wide  scatter  in  azimutli,  is  discussed  below  in  Section  6. 


B2>0,By  >0 

Four  hundred  drift  vectors  met  this  criterion.  No  clear  statistically  stable 
convection  direction,  which  applies  to  the  whole  day,  is  evident  in  these  data 
(Figure  11,  bottoiii  panel).  During  certain  periods  the  convection  direction  is  anti¬ 
sunward  but  at  other  times  it  is  sunward.  Figure  1 1  (top  panel)  fails  to  show  any 
clear  pattern  which  might  be  compared  to  cither  the  distorted  two  cell  pattern  of 
Heppner  and  Maynard  (1987)  or  a  three  or  four  cell  model  (eg  Potemra  et  al, 
1984). 
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In  this  IMF  category  only  42%  of  the  drift  directions  estimates  had  an  associated 
error  of  <5°  (Figure  8d);  a  value  similar  to  that  obtained  for  the  Bz>0,  By<0  data 
set. 


6.  THE  INTERPRETATION  OF  POLAR  CAP  DIGISONDE  DRIFT 

VECTORS  WHEN  Bz  >  0 

When  Bz  is  negative  we  have  demonstrated  that,  at  least  on  a  statistical  level,  there  is 
good  agreement  between  the  Digisonde  drift  directions  and  model  convection 
directions.  When  Bz  is  positive  we  also  propose  to  tentatively  interpret  well 
ordered  variations  in  the  drift  direction  as  a  measure  of  the  convection  direction. 
Wliilst  we  are  unable  to  support  tliis  assumption  by  comparison  with  measurements 
using  other  techniques  we  contend  that  there  is  no  reason  not  to  do  this. 

Wlien  Bz>0  the  ordered  convection  direction  measurements,  selected  according  to 
tlie  IMF  criteria,  divide  into  two  groups.  Some  vaiy'  slowly  with  time  and  examples 
of  these  occurred  on  Day  294,  1987  and  Day  344,  1988  when  By<0.  We  designate 
those  data  which  exhibit  this  behaviour  OS  (ordered  and  slowly  varying).  On  other 
occasions  the  drift  direction  variation  is  also  well  ordered  but  the  variation  in 
direction  with  time  is  fast.  An  example  of  such  a  case  was  given  in  the  post  noon 
period  of  day  042  (Figure  3).  Such  points  are  spread  over  a  wide  range  of  azimuths 
and  a  large  number  of  points,  . so  selected  from  many  days,  has  a  tendency  to  appear 
disordered.  These  points  which  originate  from  periods  when  the  azimuthal 
variation  is  ordered  but  quickly  varying  are  designated  OQ. 


In  addition  to  these  periods  of  ordered  variation  in  convection  direction  there  are 
other  periods  when  the  convection  direction  appears  to  vary  in  a  disordered 
manner  from  one  5  minute  estimate  to  the  next.  After  selection  (according  to  IMF 
criteria)  these  individual  drift  directions,  designated  D  (disordered),  are 
indistinguishable  from  the  00  vectors  except  perhaps  by  the  magnitude  of  their 
associated  errors  which  are  expected  to  be  larger  for  D  type  measurements.  An 
example  of  these  D  variations  was  given  in  Figure  3  between  00  and  06  UT,  day 
042.  One  reason  for  these  disordered  Digisonde  drift  vectors  was  proposed  by 
Buchau  et  al.(1988).  The  Digisonde  drift  analysis  algorithm  assumes  uniform  drift 
across  the  sky  and  does  not  allow  for  flow  variations  within  the  antenna  beam;  if 
flow  irregularities  do  exist  the  algorithm  produces  speeds  and  directions  which  are 
not  meaningful.  Carlson  et  al.  (1984)  have  shown  that  velocity  shears  occur  either 
side  of  sun  aligned  arcs  and  the  latter  are  known  to  occur  when  Bz  is  positive 
(Gussenhoven,  1982).  It,  therefore,  seems  likely  that  visual  or  sub-visual  sun 
aligned  arcs  were  contained  within  the  Digisonde  antenna  beam  when  D  type 
measurements  were  made.  A  zeroth  order  measure  of  the  probability  of  making  D 
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measurements  were  made.  A  zeroth  order  measure  of  the  probability  of  making  D 
type  measurements,  rather  than  OS  or  OQ  measurements  can  be  made  by 
examining  the  relative  errors  associated  with  the  individual  measurements  when 
Bz<0  as  opposed  to  when  Bz>0.  When  Bz<0  the  errors  were  <5°  in  azimuth  for 
about  70%  of  the  time  where  as  when  Bz>0  this  decreased  to  only  -45%  of  the 
time.  This  degradation  in  measurement  accuracy  indicates  that  v/hen  Bz>0  the 
Digisonde  measured  significant  velocity  shears  ~25%  more  of  the  time  than  when 
Bz<0. 

Tlie  summed  effect  of  the  OQ  type  variations  and  the  D  variations,  rather  than  the 
D  variations  alone,  are  central  to  the  reason  why  we  au  unable  to  determine  any 
average  convection  direction  when  Bz>0.  With  only  ~20  to  30  %  of  the  selected 
measurements  impaired  by  velocity  shears  a  diurnal  variation  should  be 
discemable  given  a  good  temporal  distribution  of  data  points  over  the  day.  No  such 
variation  is,  however  seen. 

The  existence  of  OQ  variations  also  leads  us  to  suggest  that  the  synthesis  of 
meaningful,  average  models,  categorised  simply  by  the  sign  of  the  IMF  components 
may  not  be  possible  when  Bz  is  positive.  Inspection  of  the  32  days  of  data  studied  in 
this  paper  reveals  that  the  azimuth  of  OQ  variations  may  be  tied,  in  a  complicated 
manner,  to  the  strength  of  the  IMF  components. 

Another  potential  instrumentation  and/or  data  analysis  error  which  is  considered  to 
produce  second  order  effects  for  polar  cap  observations,  derives  from  the  inability 
of  the  Digisonde  to  distinguish  between  electric  field  induced  F  region  movements 
and  movements  induced  by  other  forces.  The  major  inaccuracy  in  this  regard  is 
believed  to  derive  from  gravity  waves  which  will  be  most  obvious  when  the 
induced  gravity  wave  ionospheric  velocities  are  large,  or  comparable  to,  the 
electric  field  induced  drifts.  When  Bz  is  negative  the  polar  cap  electric  field  induced 
drift  velocities  are  generally  a  few  hundred  metres  per  second  or  higher  and  the 
electric  field  induced  movements  are  expected  to  dominate.  When  Bz  is  positive  the 
electric  field  induced  drift  velocity  falls  to  lower  values  and  these  other  effects  may 
disturb  the  data  on  occasion. 


7.  CONCLUSIONS 


Thirty-two  day.s  of  data  have  been  analysed  and  Fi 
results. 


cure  12  cives  a  synopsis  of  the 


WTien  Bz  is  negative  and  when  By  is  positive  and  negative  these  data  indicate  a 
clear  shift  in  the  convection  drift  direction  respectively  positive  and  negative  of  the 
anti-sunward  direction  .  Heppner  and  Maynard  (1987)  propose  two  models  for  tlie 
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negative  By  case  -  one  with  the  stream  lines  flowing  anti-sunward  (ie  0°)  and  the 
other  offset  by  -15°.  Our  measurements  indicate  a  single  median  direction  of  -12°  . 
For  the  positive  By  case  Heppner  and  Maynard  give  an  offset  of  ~45°  positive 
whereas  our  measurements  indicate  a  modal  convection  direction  of  36“ .  We 
conclude  that  when  Bz  is  negative  the  ground  based  Digisonde  has  been  shown  to 
produce  average,  daily  polar  cap  drift  directions  which  are  close  to  those  reported 
by  Heppner  and  Maynard  (1987)  from  satellite  measurements.The  usefulness  of  the 
Digisonde  for  making  polar  cap  convection  measurements,  when  Bz<0,  has, 
therefore,  been  demonstrated  when  data  are  averaged  over  several  hours  to  several 
days.  Generally  there  is  also  consistency  of  drift  azimuth  from  individual 
measurement  to  individual  measurement  which  is  suggestive  that  the  Digisonde 
can  be  used  for  case  studies  (Buchau  et  al.,1988),  however,  this  aspect  has  not  been 
addressed  in  this  paper . 

When  Bz  is  positive  and  By  negative  u,e  data  divides  into  three  groups  -  convection 
directions  which  are  ordered  and  vary  slowly  with  time  (OS),  convection  directions 
which  are  ordered  and  which  vary  quickly  with  time  (OQ)  and  convection 
directions  which  are  disordered  (D)  wiA  time.  The  latter  are  believed  to  be  due  to 
the  breakdown  of  the  measu  sinent  procedure  as  a  result  of  velocity  shears.  Using 
the  analysis  selection  procedures  described  in  this  paper  OQ  and  D  are 
indistinguishable  from  each  other  but  ther  is  hope  that  this  might  be  overcome  by 
examining  the  errors  associated  wiih  ne  individual  5  minute  resolution 
measurements.  For  the  small  sub-group  OS  there  is  a  favoured  convection 
direction  centred  on  between  0°  and  -20“  from  the  anti-sunward  during  the  periods 
2330  to  02  MLT  and  17  to  2030  MLT.  This  flow  direction  is  in  good  agreement 
with  the  distorted  two  cell  (satellite)  data  model  of  Heppner  and  Maynard  (1987) 

For  Bz  positive  and  By  positive  the  original  drift  data  divides  into  two  groups  OQ 
and  D  and  no  diurnal  variation  of  drift  direction  can  be  identified. 

The  existence  of  OQ  variations,  which  represent  significant  variations  in  drift 
directions  over  inten'als  of  5  to  10  minutes,  illustrates  one  of  tlie  problems  of  using 
satellites  which  take  1 0  minutes  or  more  to  traverse  the  polar  cap  and  which  require 
temporal  stability  over  that  period.  The  low  probability  of  OQ  variations  when 
Bz<0  (which  may  also  be  infered  from  the  relatively  narrow  drift  direction  spread 
of  ±20°  from  tlie  daily  modal  direction)  indicates  that  less  severe  problems  are  to  be 
expected  by  temporal  averaging  of  convection  data  during  these  IMF  conditions. 

The  magnitude  and  direction  of  plasma  drift  has  not  been  quantified  on  a  short  time 
scale  witli  respect  to  other  techniques.  In  the  polar  cap  tliis  involves  comparison  of 
the  convection  data  from  Qaanaaq  with  satellite  or  ground  based  optical 
measurements.  In  the  auroral  region  it  awaits  the  analysis  of  data  from  another 
Digisonde  at  Sondre  Stromfjord  and  comparison  with  data  from  the  co-located  ISR. 
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At  mid  latitudes  it  involves  comparative  measurements  at  Millstone  Hill  between  a 
Digisonde  and  the  ISR.  Due  to  the  assumption  that  the  ionospheric  bulk  velocity 
remains  constant  over  the  antenna  field  of  view  and  the  varying  velocity  shears  that 
are  expected  in  each  of  these  regions  a  series  of  measurement  comparisons  and 
improvements  in  the  analysis  technique  will  be  required  to  validate  the 
measurement  technique.  The  Digisonde  drift  measurements  are  also  susceptible  to 
movements  induced  by  other  ionospheric  dynamic  effects,  such  as  travelling 
ionospheric  disturbances  generated  by  neutral  gravity  waves,  which  may  induce  a 
degree  of  variability  on  the  calculated  convection  velocities.  The  percentage  error 
will  be  a  function  of  ih.e  induced  plasma  convection  velocity  and  may  be  expected  to 
be  small,  under  the  expected  high  plasma  velocities  at  auroral  and  polar  cap 
latitudes,  but  relatively  more  important  at  middle  latitudes. 
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Table  1 

Selected  Data  Periods 


Year 

Day  Number 

Date 

1986 

302-305 

29  Oct-1  Nov 

1987 

139-141 

19-21  May 

167-168 

16-17  June 

265-267 

22-24  September 

293-294 

20-21  October 

328-330 

24-26  November 

356 

22  December 

1988 

039-042 

8-11  February 

052-055 

21-24  February 

340-345 

5-10  December 

Figure  1 
Figure  2 


Figure  3 


Figure  4 

4 


Captions 

Receiver  antenna  geometry. 

A  skymap  survey  for  part  of  day  052.  Each  panel 
shows  the  positions  of  oblique 
ionospheric  sources  scaled  in  degrees  with 
zenith  angles  varying  from  0®  to  ±80®.  Under  each 
panel  is  shown  the  universal  time,  the  number  of 
drift  soundings  contributing  to  the  panel  and  the 
number  of  sources  mapped  onto  the  panel. 

Drift  directions  for  days  040  to  042, 1988.  OS 
indicates  ordered  and  slowly  vai^dng  drift  direction 
variations;  00  ordered  and  quickly  varying  and  D 
disordered  variations,  see  Section  6. 

IMF  Bz  component  (nT)  at  IMP8  on  days  040  to  042, 
1988 


Figure  5 


IMF  By  component  (nT)  at  IMPS  on  davs  040  to  042, 
1988 
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Figure  6 


Figure  7 


Figure  8 


Ficure  9 


Figure  10 


Figure  11 


Figure  12 


No  of  drift  vectors  selected  per  day  per  IMF 
category. 

Selected  drift  data  when  Bz<0  and  By<0  :top  panel, 
temporal  distribution  of  selected  data,  bottom 
panel,  deviation  of  convection  direction  from  the 
anti-sunward. 

Diurnal  variation  of  the  quartile  error  of  each  (five 
minute  resolution)  drift  direction  estimate:  a) 

Bz<0,  By<0;  b)  Bz<0,  By>0;  c)  Bz>0,  By<0;  d 
Bz>0,  By>0. 

Selected  drift  data  when  Bz<0  and  By>0:  Top  panel, 
temporal  distribution  of  selected  data;  bottom  panel 
deviation  of  convection  direction  from  tlie  anti- 
sunward. 

Selected  drift  data  when  Bz>0  and  By<0:  Top  panel, 
temporal  distribution  of  selected  data;  bottom  panel 
deviation  of  convection  direction  from  the  anti- 
sunward. 

Selected  drift  data  when  Bz>0  and  By>0:  Top  panel, 
temporal  distribution  of  selected  data;  bottom  panel 
deviation  of  convection  direction  from  the  anti- 
sunward. 

Synopsis  of  measured  convection  directions:  polar 
plot,  shading  indicates  no  clear  convection  direction. 
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Figure  2 
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6.  2 
OQ  Variations 


141 


11/13/89 


2 


'-1B0. 


Universal  ^  Time 


Figure  7 
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Figure  8 


Percentage  Deviation  (degrees) 


Figure  10 
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Percentage  Deviation  (degrees) 


Drift  direction  -  Antisunward  direction  (Degrees) 


Figure  11 
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Figure  12.  Synopsis  of  measured  convection  directions:  polar  view, 
shading  indicates  no  clear  convection  direction 
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Plasma  convection  measurements,  using  Digisonde  ionospheric  sounders,  have  been  conducted  in 
the  central  polar  cap  at  Thule,  Greenland  (86'  CGL)  and  more  recently  at  Qaanaaq,  Greenland  (87° 
CGL),  in  the  auroral  oval  at  the  Air  Force  Geophysics  Laboratory  Goose  Bay  Ionospheric  Observa¬ 
tory  (65°  CGL).  and  at  suboval  latitudes  at  Argentia  NAS  (57°  CGL).  The  plasma  convection  or 
ionospheric  drift  measurements  conducted  at  Thule  and  Qaanaaq  during  campaigns  from  the  1983- 
1984  winter  to  the  present  provide  evidence  that  antisunward  convection  dominates  in  the  polar  cap 
with  velocities  typically  between  200  and  900  m  s"‘.  Velocity  reversals  or  shears  were  observed  in 
association  with  polar  cap  F  layer  auroras  during  quiet  magnetic  conditions.  Observations  of  the 
plasma  drift  at  Goose  Bay  show,  as  expected,  a  drift  reversal  from  westward  to  eastward  around 
midnight  CGLT,  indicating  the  rotation  of  Goose  Bay  from  the  evening  into  the  morning  convection 
cell.  Observations  at  Argentia,  typically  a  suboval/trough  station,  provide  evidence  under  magnetically 
disturbed  conditions  for  the  midnight  reversal  of  the  antisunward  flow  pattern.  However,  the  data  at 
Argentia  are  generally  less  consistent  under  magnetically  quiet  conditions.  This  likely  indicates  the 
proximity  of  Argentia  to  the  boundary  between  -'tating  and  convecting  plasma.  Recent  theoretical 
calculations  of  electron  density  profiles  within  the  iu,jh-latitude/polar  cap  ionosphere  demonstrate  that 
the  diurnal  f,  F,  variation  observed  at  Thule  is  controlled  by  the  plasma  convection  pattern  and  the 
associated  drift  velocities.  The  model  calculations  for  5,  <  0  and  B,ts0  show  factors  of  2  to  3 
differences  in  over  Thule,  supporting  the  stated  importance  of  convection  pattern  and  velocity 
t.-ieasurements  for  the  modeling  of  the  high-latitude  ionosphere. 


1.  INTRODUCTION 

The  high-latitude  ionosphere  exhibits  several 
large-.scale  characteristics,  for  example,  winter  polar 
cap  UT  dependence  [Saro  and  Rourke,  1964;  Buchau 
et  ai,  1985]:  longitudinal  variation  in  the  F  region 
densities  in  the  nightside  auroral  zone  [</e  la  Beau- 
jardiere  et  al.,  1985];  and  the  main  F  layer  trough. 
These  large-scale  characteristics,  at  least  to  a  large 
extent,  are  a  result  of  the  plasma  convection  across 
the  polar  cap  [_Knudsen,  1974;  Knudsen  et  ai,  1977; 
Watkins,  1978;  Sojka  et  al.,  1981]  and  the  displace- 
iTient  of  the  geomagnetic  pole  from  the  geographic 
pole  ISojka  et  ai,  1979]. 

The  existence  of  polar  plasma  convection  was  in¬ 
ferred  from  electric  field  measurements  from  satellites 
ICauffman  and  Gurnett,  1972],  or  by  directly  mea- 
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suring  the  convection  with  satellite-borne  plasma 
drift  meters  IHanson  and  Heelis,  1975],  Ground- 
based  high-latitude  plasma  convection  measurements 
have  been  conducted  during  the  decade-long  mea¬ 
suring  campaign  of  the  Chatanika  incoherent  scatter 
radar  (ISR)  IFoster,  1983]  and  by  the  northward 
looking  Millstone  Hill  ISR  [Feans  et  ai,  1980].  The 
large  body  of  plasma  convection  data  shows  that  in 
its  simplest  form  the  convection  occurs  in  a  two-cell 
pattern.  Models  of  this  two-cell  convection  pattern 
developed  by  Volland  [1975]  and  improved  by 
Heelis  et  al.  [1982]  have  been  used  extensively  to 
model  the  ionosphere  under  the  influence  of  the  con¬ 
vection  (above  references).  Figure  1  [from  Heelis  et 
ai,  1982]  shows  the  two-cell  configuration,  with  anti¬ 
sunward  plasma  flow  across  the  polar  cap  and  sun¬ 
ward  return  flow  along  the  morning  and  evening 
flanks  of  the  auroral  oval. 

Past  observations  and  modeling  work  have  shown 
that  knowledge  of  the  convection  pattern  is  a  pre¬ 
requisite  for  understanding  the  polar  ionosphere.  The 
two-cell  pattern  occurs  most  clearly  when  the  inter¬ 
planetary  magnetic  field  (IMF)  has  a  southward 
component  (B.  <  0),  while  the  exact  flow  geometry  is 
controlled  by  *B,.  [Heelis,  1984],  For  B.  >  0,  dusk  to 
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IHSUNUNEOUS  CLOItL  EATTEKN 


Fig.  1.  Polar  plasma  convection  pattern  [Heelis  ei  at.,  1982]  in 
CGL/CGLT  coordinate  system.  The  locations  of  the  drift  mea¬ 
suring  stations  Qaanaaq  (Q),  Goose  Bay  (G).  and  Argentia  (A)  are 
indicated  by  dashed  circles. 


dawn  polar  cap  electric  fields  have  been  measured, 
which  cause  sunward  plasma  flow  across  the  central 
polar  cap  [Burke  et  al,  1979],  leading  to  the  model 
of  a  four-cell  convection  pattern.  Satellite  magnetom¬ 
eter  measurements  [Zanetti  et  al.,  1984]  suggest  a 
three-  or  four-cell  pattern  depending  on  the  sign  of 
By.  Lassen  [1979]  interpreted  the  polar  cap  auroral 
arcs,  observed  during  magnetically  quiet  periods,  as 
being  the  result  of  a  shear  associated  with  a  third 
central  polar  cap  cell,  with  the  shear  occurring  either 
at  the  edge  of  the  morning  (when  By  <  0)  or  the  eve¬ 
ning  cell  (when  By  >  0).  Recent  analysis  of  Dynamics 
Explorer  (DE)  2  data  [Heppner  and  Maynard,  1987] 
indicates  that  it  is  possible  to  <plain  the  data  for  the 
full  range  of  B.  and  By  with  a  basic  two-cell  pattern. 
For  the  B.  >  0  condition  this  approach  results  in 
substantial  distortion  of  the  two  cells  to  explain  sun¬ 
ward  flow  in  the  central  polar  cap.  However,  the 
resulting  pattern  agrees  with  observations  in  other 
regions  of  the  high-latitude  ionosphere. 

In  this  paper  we  describe  a  ground-based  coherent 
HF  radar  (digital  ionosonde)  technique,  which  allows 
continuous  measurement  of  the  plasma  convection  at 
selected  sites.  Early  observations  [Reinisch  et  al, 
1987]  show  good  agreement  between  the  measure¬ 
ments  and  the  major  features  of  the  plasma  convec¬ 
tion  patterns.  The  analysis  presented  here  indicates 
that  not  only  are  the  major  features  observed,  but 
also  details  such  as  shears  and  drift  reversals  in  re¬ 
sponse  to  quiet  magnetic  conditions  are  detected. 


The  digital  ionosonde  measurements  suggest  that 
continuous  monitoring  of  at  least  the  large-scale  fea¬ 
tures  of  the  convection  pattern  with  a  ground-based 
technique  is  possible  and  useful. 

2.  SPACED  ANTENNA  DOPPLER 
DRIFT  TECHNIQUE 

The  Digisonde  [Bibl  and  Reinisch,  1978]  drift  tech¬ 
nique  uses  spectral  analysis  of  ionospheric  echoes 
[PJister,  1974fl,  b].  This  technique  has  been  refined 
for  automatic  data  processing  [Dozois,  1983]  and  is 
an  integral  part  of  all  modern  Digisondes  [Reinisch, 
1986]. 

The  drift  measurements  are  conducted  using  phase 
coherent  pulse  transmission  on  two  HF  frequencies 
which  are  selected  to  sample  desired  heights  in  the  F 
region.  These  transmitted  HF  waves  typically  illumi¬ 
nate  a  large  overhead  F  region  area  of  several 
hundred  kilometers  diameter.  An  array  of  four  or 
seven  antennas  (depending  on  the  site)  deployed  in  a 
100-m  baseline  triangle  receives  the  signals  reflected 
from  the  ionosphere  (the  technique  is  equally  appli¬ 
cable  to  £  region  studies).  The  individual  antennas  of 
the  receiving  array  are  multiplexed  at  the  pulse  repe¬ 
tition  rate  (typically  2(X)  Hz).  For  each  of  the  two 
frequencies  the  time  series  of  the  echo  received  and 
sampled  at  each  of  the  antennas  is  Fourier  trans¬ 
formed  in  real  time,  resulting  in  four  (or  seven)  com¬ 
plex  spectra  at  the  end  of  each  measurement  period, 
one  for  each  antenna. 

During  a  normal  measurement  campaign,  iono- 
gramrs  are  taken  in  5-min  intervals  to  monitor  proper 
frequency  and  sample  gate  selection  in  a  dynamic 
environment.  The  time  between  ionograms  is  used  to 
collect  a  series  of  (nominally  eight)  drift  measure¬ 
ments  of  typically  20-s  each.  Figure  2  shows  an  ex¬ 
ample  of  such  a  sequence  of  ionogram/drift  measure¬ 
ments.  On  the  left  is  the  representation  of  the  digital 
ionogram  (bottom,  amplitude;  top,  status,  i.e.,  coded 
information  on  Doppler/direction/polarization);  and 
to  the  right,  a  subset  of  the  collected  Doppler  spectra, 
displayed  to  permit  inspection  of  data  quality  during 
data  acquisition.  Data  time  runs  from  left  to  right  as 
indicated.  Spectral  frequencies  are  displayed  verti¬ 
cally,  with  0,  plus,  and  minus  Doppler  as  indicated. 
The  bottom  right  half  of  Figure  2  (labeled  “ampli¬ 
tude”)  gives  the  amplitude  spectrum;  the  top  half  (la¬ 
beled  “phase”)  shows  the  phase  of  each  correspond¬ 
ing  Doppler  line;  jointly  these  two  data  sets  provide 
a  set  of  complex  spectra.  For  a  20-s  sampling  period 
the  Doppler  resolution  is  0.05  Hz,  and  the  Doppler 
range  is  ±1.6  Hz  with  64  Doppler  lines.  These  values 
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FREQUENCY - TIME - ^ 

lONOGRAM  MODE  DRIFT  MODE 

Fig.  2.  A  lypic.ll  example  of  a  sequence  of  ionogram  and  drift  measurements  collected  during  a  drift  campaign  at 
Argenlia,  Newfoundland.  Such  a  data  set  is  collected  every  5  min. 


can  be  changed  over  a  certain  range  to  adjust  to 
changing  ionospheric  conditions  by  selecting  differ¬ 
ent  sampling  rates  (pulse  repetition  frequencies)  and 
sampling  windows. 

Spatial  cross  correlation  IBibl  et  ai,  1975]  of  the 
complex  spectra  (see  Figure  3)  determines  the  angle 
of  arrival  for  each  spectral  component  containing 
significant  energy.  As  a  result  of  this  analysis  one  can 
construct  a  “sky  map"  showing  the  location  of  each 
reflection  point,  or  source  [Dozo/s,  1983],  specified 
by  a  given  Doppler  frequency.  The  Doppler  fre¬ 
quency  provides  the  radial  velocity  component  of  the 
moving  plasma  for  this  source.  The  returns  from  ail 
sources  superimpose  as  a  combined  signal  at  each 


A  iiw  iiwvjuuiivj’  id  29*^^** 


d,  =  (l/ji)  •  V  -  k,  j  =  I,  2,  •  •  ■ 


where  s  is  the  source  index,  v  is  the  drift  velocity,  and 
k  the  wave  vector  from  source  s.  The  sky  map  covers 


the  area  of  the  antenna  array’s  main  lobe,  a  fre¬ 
quency  dependent  area,  nominally  smaller  than  400 
km  along  each  side  for  F  region  reflections  (<45° 
half  beam  width).  In  the  sky  map  of  Figure  3  the 
majority  of  the  sources  are  in  the  NW  quadrant;  the 
letters  indicate  positive  and  the  numbers  negative 
Doppler  shifts.  By  assuming  that  the  observed  Dop¬ 
pler  shifts  are  the  result  of  a  uniform  motion  of  the 
reflecting  plasma,  one  can  determine  the  three- 
dimensional  velocity  vector  v  which  in  the  least 
squares  error  sense  best  represents  the  Doppler  fre¬ 
quencies  d,  measured  at  the  source  locations  specified 
by  k, .  The  drift  direction  for  the  sky  map  in  Figure  3 
is  approximately  to  the  east-northeast. 

The  size  (or  diameter)  of  the  convection  pattern 
and  the  actual  convection  velocities  are  important 
quantities  in  modeling  the  polar  ionosphere.  Three 
Digisonde  stations,  Qaanaaq,  Greenland  (87°  CGL), 
Goose  Bay,  Labrador  (65°  CGL),  and  Argentia, 
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Cross  correlation  P  of  complex  Doppler  spectrum  F  from  4  antennas: 
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yields  angle  of  arrival  of  each  reflected  spectral  line. 
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Fig.  3.  The  identification  of  the  location  of  reflecting  areas  (sources)  for  those  Doppler  components  of  the 
ionospheric  returns  which  contain  energy  above  a  threshold  level  is  made  by  spatial  cross  correlation.  The  result  is 
a  sky  map  (bottom)  of  source  locations.  Letters  indicate  sources  with  positive,  numbers  those  with  negative 
Doppler.  The  three-dimensional  drift  vector  is  obtained  by  least  squares  fit.  assuming  uniform  drift  over  the  sky 
map  area. 


NAS,  Newfoundland  (57°  CGL)  provide  a  first  step 
at  continously  measuring  features  of  the  convection. 
Their  respective  locations  over  a  full  day  are  indicat¬ 
ed  in  Figure  1  by  dashed  circles  marked  Q  (Qaa¬ 
naaq),  G  (Goose  Bay),  and  A  (Argentia).  The  figure 
shows  that  for  the  conditions  represented  by  this  pat¬ 
tern,  Qaanaaq  is  fuiiy  embedded  in  the  antisunward 
flow  in  the  central  polar  cap,  while  Goose  Bay  is 
under  the  influence  of  westward  (sunward)  flow 
before  midnight  and  of  eastward  (sunward)  flow  after 


midnight.  Under  quiet  conditions,  Argentia  is  equa- 
torward  of  the  high-latitude  convection  and  observes 
mid-latitude  corotating  plasma.  Under  more  active 
magnetic  conditions  (in  response  to  B.  <  0)  the  con¬ 
vection  pattern  expands  in  diameter,  bringing  Argen¬ 
tia  under  its  influence  for  part  of  the  night  and  ex¬ 
tending  the  duration  during  which  Goose  Bay  ob¬ 
serves  sunward  return  flow.  Figure  4  presents  diurnal 
variation  of  the  drift  direction  in  a  corrected  geomag¬ 
netic  (CG)  direction, CG  local  time  diagram  which 


152 


BUCHAU  ET  AL.:  POLAR  CAP  MEASUREMENTS  AND  REAL-TIME  MODELING 


525 


TIME  HISTORIES  AT  POLAR  AND  AURORAL  STATIONS 

Fig.  4.  Expected  diurnal  change  of  the  drift  direction  (as¬ 
suming  the  convection  shown  in  Figure  1)  for  Qaanaaq  and 
Goose  Bay  in  a  CG  direction  and  CG  local  time  coordinate 
system. 

would  be  observed  at  Qaanaaq  (top)  and  Goose  Ba> 
(bottom),  assuming  the  convection  pattern  of  Figure 
1.  At  Qaanaaq  the  drift  direction  rotates  linearly 
through  360°,  while  at  Goose  Bay  the  westward  drift 


before  midnight  changes  to  eastward  at  or  near  mid¬ 
night  CGLT. 

3.  DRIFT  MEASUREMENTS 

3.1.  Thule — Magnetically  active  period 

Several  winter  measurement  campaigns  were  con¬ 
ducted  by  the  Airborne  Ionospheric  Observatory 
while  on  the  ground  at  Thule  AB,  Greenland 
(86°  CGL)  to  study  the  structure  and  dynamics  of  the 
polar  cap  ionosphere.  During  these  campaigns  the 
antisunward  convection  of  large  nonauroral  plasma 
patches  was  observed  with  an  all-sky  imaging  pho¬ 
tometer  (ASIP)  [Weber  et  al.,  1984],  and  the  UT 
dependence  of  the  maximum  electron  densities  in 
these  patches  was  established  [Bitchau  ct  al.,  1985]. 
For  the  campaigns  in  December  1983,  January- 
February  1984  and  March  1985  a  four-antenna  array 
was  deployed  near  the  aircraft  to  provide  measure¬ 
ments  of  the  plasma  drift  independent  of  the  optical 
observations.  The  patches  observed  by  the  ASIP 
under  active  conditions  in  the  1200  to  00{X)  UT  time 
window  permit  an  independent  measurement  of  the 
convection  velocities.  Between  0000  and  1200  UT  the 
patches  are  generally  not  detectable  by  the  ASIP 
owing  to  the  low  plasma  densities  [Buchau  et  al., 
1985].  Drift  results  (Figure  5)  for  measurements 
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Fig.  5.  F  region  drift  direction  (bottom)  and  velocity  (top,  showing  mean  velocity  and  range  of  observations) 
taken  during  four  moderately  disturbed  days  in  January  and  February  1984  at  Thule  AB.  Greenland.  The  anti¬ 
sunward  direction  is  indicated  for  reference. 
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taken  at  Thule  during  four  disturbed  to  moderately 
disturbed  days  between  January  30  and  February  4, 
1984  {ZKp  for  January  30  and  31  and  February  1-4 
was  34.  29  —  ,  23 -f,  22-,  23,  and  41-1-,  respectively) 
show  a  predominantly  antisuiiward  drift  (the  anti- 
sunward  direction,  rotating  through  360  in  24  hours, 
IS  indicated).  This  agrees  with  the  analysis  of  simulta¬ 
neous  optical  observations  for  February  3-4,  1984 
[HVfrer  et  al..  1986]  and  earlier  results  derived  from 
optical  data  of  antisunward  plasma  transport  under 
disturbed  conditions  [Weber  et  ai,  1984],  The  ob¬ 
served  velocities  shown  in  the  top  panel  of  Figure  5 
(mean  velocity  and  range)  fluctuate  between  150  and 
900  m  s"',  well  within  the  range  of  previous  satellite 
measurements  and  in  general  agreement  with  the  op¬ 
tical  observations.  With  the  confidence  in  the  tech¬ 
nique  gained  from  these  measurements,  permanent 
digital  sounding  capabilities  were  established  at  the 
Qaanaaq  Geophysical  Observatory  of  the  Danish 
Meteorological  Institute,  in  northern  Greenland 
(87*'  CGL),  and  at  a  suboval  station  (Argentia,  New¬ 
foundland.  58'  CGL).  Together  with  the  Air  Force 
Geophysics  Laboratory  (AFGL)  Goose  Bay  Iono¬ 
spheric  Observatory  (65‘  CGL)  these  stations  form  a 
meridian  chain  along  approximately  20°E  CG  lon¬ 
gitude. 

3.2.  Thule — Miuinetkally  quiet  period 

Earlier  campaigns  had  established  that  sun-aligned 
polar  cap  arcs,  especially  subvisual  F  layer  arcs, 
occur  in  the  polar  cap  under  magnetically  quiet  con¬ 
ditions  [Weber  and  Bucliau,  1981;  Buchau  et  al., 
1983].  Reiff  et  al.  [1978]  showed  that  current  conti¬ 
nuity  at  plasma  flow  shears  requires  upward  (or 
downward)  Birkcland  currents,  therefore  producing 
subvisual  and  visual  signatures,  i.e.,  auroras  in  the 
region  of  upward  currents.  Lassen  [1979]  concluded 
that  polar  cap  arcs  arc  optical  signatures  of  shears 
and  flow  irregularities  in  the  magnctosphcric  convec¬ 
tion.  The  analysis  of  Sondrestrom  ISR  data  in  the 
context  of  ASIP  arc  observations  has  confirmed  that 
arcs  delineate  plasma  shears  [Carlson  et  al.,  19o4]. 

The  Thule  drift  data  taken  on  December  9-10, 
1983  cover  a  9-hour  period  at  the  end  of  a  27-hour- 
long  quiet  period  (December  9  "EKp  -  5  +  ).  During 
the  quiet  period  from  19(X)  to  2300  UT  on  December. 
9,  1983,  subvisual  and  visual  arcs  typical  for  quiet 
periods  were  observed  in  the  ASIP  data  to  occupy 
major  parts  of  the  sky  and  to  drift  from  dawn  to 
dusk.  The  concurrent  drift  measurements  are  shown 
in  Figure  6  (bottom).  The  drift  direction  oscillates 


dramatically  throughout  the  three  magnetically  quiet 
3-hour  periods  at  the  beginning  of  the  data  sample 
{Kp  is  indicated  below  the  abscissa),  suggesting  the 
transit  of  several  shear  boundaries  through  the  sta¬ 
tion  zenith.  As  magnetic  activity  increases,  the  con¬ 
vection  direction  becomes  more  orderly,  with  a  gen¬ 
eral  direction  30‘’-45'’  clockwise  from  antisunward. 
The  drift  velocities  are  low  (150-300  ms"')  during 
the  quiet  period  and  increase  somewhat  (150-600  m 
s" ')  during  the  disturbed  period. 

Also  of  significance  here  is  the  first  comparison  of 
drift  data  with  IMF  data,  shown  in  the  top  part  of 
Figure  6.  The  large  fluctuations  in  drift  direction, 
indicating  the  presence  of  shears,  as  well  as  the  arc 
observations  occur  during  times  when  >  0,  B,,  « 
0,  and  B^  >  0.  As  B^  decreases  and  B,.  increases  with 
still  (strongly)  positive  B.,  the  plasma  drift  becomes 
more  orderly,  but  skewed  by  60”-90°  clockwise  from 
antisunward.  For  the  steady  B.  <  0  conditions  from 
1100  to  2100  UT,  the  flow  is  steady  antisunward 
(skewed  30“  clockwise).  As  B.  becomes  positive  again 
after  2100  UT,  shears  start  to  dominate  again.  No 
clear  cflects  of  B,.  and  B,  on  the  drift  direction  (under 
B.  <  0  conditions)  can  be  derived  from  this  short 
measurement. 

For  the  drift  analysis  the  convection  is  assumed  to 
be  uniform  across  the  sky  map.  This  clearly  is  not  the 
case  during  the  9-hour  period  o<‘  arc  observation,  and 
the  resulting  drift  directions  and  the  calculated  veloc¬ 
ities  are  not  meaningful  but  merely  give  an  indication 
of  the  existence  of  shears.  A  preliminary  attempt  to 
separate  sky  maps  into  regions  within  and  on  either 
side  of  arcs,  using  ASIP  images  as  a  guide,  shows 
evidence  of  shears.  A  detailed  discussion  of  this 
analysis  is  beyond  the  scope  of  this  paper. 

3.3.  Qaanaaq — First  extended  drift  measurements 

Drift  measurements  until  early  1986  had  been 
made  by  manually  selecting  frequencies  and  setting 
sampling  gales,  thereby  requiring  continuous  man¬ 
ning  of  the  Digisonde.  Since  1986  this  procedure  has 
been  automated  for  all  Digisonde  256  stations.  The 
automation  is  based  on  the  ARTIST  trace  identifi¬ 
cation  software  for  the  Digisonde  ionograms  [Re/n- 
isc/i  and  Huang,  1982;  Reinisch  et  al.,  1982]  that  re¬ 
sides  in  an  IBM  AT  microcomputer  which  is  part  of 
the  Digisonde  256. 

The  data  (Figure  7)  cover  a  continuous  period  of 
almost  6  days  of  drift  measurements,  from  October 
29  to  November  3.  1986.  A  cursory  inspection  sug¬ 
gests  at  least  for  the  first  four  days  steady  anti- 
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r  REGION  DRIFT 

OIQISONDE  OBSERVATIONS  AT  TNULE,  GREENLAND 
AND 

MAGNETIC  FIELD  PARAMETERS 
•/10  DEC  t3  1S00  -  214S  UT 


Fig.  6.  First  comparison  of  F  region  drift  measurements,  taken  on  December  9-10,  1983  at  Thule,  and  inter¬ 
planetary  magnetic  field  (IMF)  daia  (for  discussion  see  text).  Also  shown  are  the  Kp  values  for  the  measurement 
period. 


sunward  drift  direction.  However,  cleai  deviations 
are  seen  for  the  period  1500-1800  UT  on  October  29, 
and  1200-1600  UT  on  November  1.  Significant  oscil¬ 
lations  in  direction,  similar  to  those  discussed  in  sec¬ 
tion  3.2,  indicate  the  presence  of  shears/arcs  starting 
1700  UT  on  November  1  and  continuing  more  or 


less  unabated  throughout  the  remainder  of  the 
measurement  period. 

Comparison  of  magnetic  activity  (Kp)  and  drift 
measurements  in  Figure  7  indicates  antisunward  con¬ 
vection  at  velocities  greater  than  300  ms"'  during 
more  disturbed  periods  (Kp  >  3)  on  October  29  and 
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Pig  7.  ConUnuous  F  region  dnfi  tneasurcmems  laken  at  Qaanaaq.  Greenland,  from  October  28,  1986  (2100 
UT)  to  November  3,  1986  (15(X)  UT),  Kp  values  are  shown  between  corresponding  3-hour  tick  marks.  Superim¬ 
posed  on  the  presentation  of  the  drift  veloutj  for  comparison  is  an  analog  representation  of  the  Kp  time  variation. 


early  October  30.  1986.  This  active  pet  tod  was  inter¬ 
rupted  by  a  3-liour  quiet  period  (Kp  =  1,  1500  to 
1800  UT  on  October  29)  during  which  the  drift 
turned  sunward,  in  agreement  with  the  previously 
mentioned  satellite  observation  for  B.  >  0.  Sunward 
drift  during  quiet  times  following  disturbed  con¬ 
ditions  was  again  observed  on  November  1,  1986. 
Magnetically  active  conditions  (Kp  =  4,  4  +  ,  3  +  , 
2  +  )  during  the  first  half  of  the  day  were  interrupted 
by  quiet  (Kp  =  1  — )  conditions  from  1200  to  1500 
UT,  during  which  period  well-organized  antisunward 
drift  with  velocities  greater  than  600  ms"'  was  mea¬ 
sured.  From  here  to  the  end  of  the  measurements  the 
magnetic  activity  remained  generally  quiet,  and  clear 
evidence  of  shears  was  observed  throughout.  A  short 
period  of  moderate  activity  (Kp  =  3  +  .  3  —  )  within 
this  quiet  period  tended  to  show  somewhat  more  an¬ 
tisunward  organization  of  the  drift.  While  the  pre¬ 
viously  discussed  measurements  from  Figure  7  agree 
with  generally  accepted  behavior  of  convection 


during  various  levels  of  magnetic  activity,  this  is  not 
the  case  for  October  31,  1986,  when  despite  extended 
quiet  magnetic  conditions  the  drift  is  consistently  an¬ 
tisunward,  with  velocities  decaying  from  ~450  to  75 
m  s" '  by  the  end  of  the  quiet  period. 

The  (horizontal)  drift  velocities  are  shown  above 
the  respective  drift  direction  panels.  Superimposed  is 
the  variation  of  Kp  as  a  function  of  time  in  analog 
presentation  (continuous  line).  There  is  evidence  that 
higher  Kp  values  correlate  with  higher  velocities  and 
lower  Kp  values  with  lower  velocities  (for  example, 
2100  UT  on  October  31  to  1500  UT  on  November 
1).  The  correlation  breaks  down  during  the  extended 
period  of  shear  observation  (second  half  of  Novem¬ 
ber  1.  and  all  of  November  2-3,  1986). 

3.4.  Qaanaaq — Drift  nwusun'incnls  under  i,unlil 
conditions 

A  limited  set  of  drift  measurements  taken  at  Qaa¬ 
naaq  under  fully  sunlit  conditions  corroborates  the 
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F  REGION  DRIFT 

DIGISONDE  OBSERVATIONS  AT  QANAQ.  GREEMUND 
7/8  APRIL  1986  1900  TO  1900  CGLT 

Fig.  8.  F  region  drift  observations  from  Qaanaaq.  Greenland,  taken  on  April  7-8,  1986  taken  under  magnetically 

quiet  conditions. 


observations  in  darkness.  Drift  measurements  con¬ 
ducted  at  Qaanaaq  on  July  29-31,  1986  (data  not 
presented  here)  show  consistent  antisunward  convec¬ 
tion  with  velocities  from  150  to  450  m  s"‘  for 
the  prevailing  moderately  active  conditions  (LKp  = 
26—,  20,  20 -f,  respectively,  or  Kp>2).  For  mea¬ 
surements  taken  under  magnetically  quiet  conditions, 
in  full  sunlight  on  April  8,  1986  (Figure  8),  with 
Kp  =  1  — ,  1  -F ,  1  + ,  2  -F  for  the  period  of  measure¬ 
ments  (0000  to  1200  UT  on  April  8,  1986  or  2200  to 
0900  CGLT)  the  convection  is  initially  sunward,  as 
expected  for  quiet  or  (inferred)  B,  >  0  conditions.  As 
magnetic  activity  increases  through  Kp  =  1  +  to  2+, 
the  drift  direction  turns  antisunward,  skewed  clock¬ 
wise  by  30°-60“. 

It  is  clear  from  the  set  of  Qaanaaq  observations 
discussed  so  far  that  Kp  is  at  best  a  coarse  ordering 
parameter  for  the  drift  observations,  especially  since 
B.  >  0  does  not  necessarily  correspond  to  quiet  mag¬ 
netic  conditions  as  Figure  6  indicates.  There  a  jump 
from  B.  =  -F  3  nT  to  -F  6  nT  corresponds  to  a  change 
in  Kp  from  Oq  to  3o ,  contrary  to  an  expectation  of 
continued  quiet  conditions. 


3.5.  Goose  Bay  observations 

As  Figures  1  and  4  indicate.  Goose  Bay,  for  the 
selected  average  convection  pattern,  is  under  the  in¬ 
fluence  of  the  sunward  return  plasma  flow  for  +6 
hours  around  CG  midnight,  with  westerly  flow  ex¬ 
pected  prior  to  and  easterly  flow  past  midnight.  At 
other  times.  Goose  Bay  is  under  the  influence  of  co¬ 
rotating  plasma,  measuring  the  lower  drift  velocities 
of  the  mid-latitude  ionosphere  [Kent  and  Wright, 
1968]. 

A  typical  example  of  Goose  Bay  drift  measure¬ 
ments  are  the  data  from  January  14-15,  1983  (Figure 
9)  taken  under  moderately  disturbed  conditions 
(LKp  =  16-f  and  28+  for  January  14  and  15,  re¬ 
spectively).  The  convection  in  the  premidnight  sector 
is  consistent  and  toward  the  west  (sunward).  After 
CG  midnight  the  convection  switches  toward  the 
east,  moving  smoothly  through  the  expected  south¬ 
erly  direction.  A  valid  measurement  indicating  a 
quick  switch  back  to  westerly  drift  at  02(X)  AST  in¬ 
terrupts  the  smooth  transition  toward  easterly  flow. 
It  should  be  noted  that  the  velocities  are  consider- 


157 


530 


BUCHAU  ET  AU.:  POLAR  CAP  MEASUREMENTS  AND  REAL-TIME  MODELING 


AST  '  a'o  '  22  '  00  '  0*2  '  o'o  '  o'o  '  o'o 

I 

F  REOION  DRIFT 

DIQISONDE  OBSERVATIONS  AT  QOOSE  BAY,  LABRADOR 
.  14/18  JAN  os  18:00  TO  0000  AST 

Fig.  9.  F  region  drift  observations  at  Goose  Bay  on  January  14-15.  1983  show  the  expected  switch  of  the  sunward 
convection  return  flow  from  westward  to  eastward  near  midnight. 


ably  lower  than  those  observed  in  the  polar  cap,  with 
premidnight  velocities  measured  at  250  m  s“‘  and 
postmidnight  velocities  at  100  m  s“‘. 

A  substantial  data  set  (38  full  or  partial  days)  col¬ 
lected  between  1981  and  1985  during  the  devel¬ 
opment  of  the  drift  technique  is  not  of  the  same  qual¬ 
ity  and  consistency  as  the  polar  cap  data  set.  This  is 
due  to  (1)  manual  data  collection,  (2)  interference 
with  the  F  layer  measurements  by  blanketing  £, 
during  magnetically  disturbed  periods,  and  (3)  lack  of 
solid  echoes  during  low-/^  Fj  (trough)  conditions,  co¬ 
inciding  with  quiet  magnetic  activity  and  small  oval 
and  convection  pattern  diameter.  But  throughout  the 
measurement  set,  consistent  premidnight  westward/ 
postmidnight  eastward  drift  direction  is  observed. 
Short  15  min  to  1  hour  drift  reversals  are,  however, 
not  uncommon,  especially  in  the  postmidnight  sector. 
The  drift  velocities  shown  in  Figure  9  are  typical  of 
the  higher  velocities  observed  in  the  full  data  set; 
however,  the  high  velocities  may  be  observed  after 
midnight  as  well  as -before.  Low  velocities  (<50  m 
s"*)  are  more  routinely  observed  and  are  more  typi¬ 


cal  for  times  away  from  midnight.  A  thorough  analy¬ 
sis  in  the  context  of  available  IMF  data  is  planned, 
focusing  specifically  on  the  time  of  the  midnight 
switch  from  west  to  cast,  the  relation  of  short-lived 
drift  reversals  to  IMF  conditions,  and  the  daily  dura¬ 
tion  of  distinct  convection  conditions.  This  latter 
measurement  may  provide  information  on  the  diam¬ 
eter  of  the  convection  pattern. 

3.6.  Argentia  observations 

Our  southernmost  measurements  have  been  taken 
at  Argentia  NAS,  Newfoundland  (57°  CGL).  Under 
quiet  and  moderate  magnetic  conditions,  Argentia  is 
likely  to  remain  under  the  influence  of  the  corotating 
ionosphere,  at  or  south  of  the  plasmapause  (Figure 
1).  For  low  magnetic  activity  {ZKp<  25),  observed 
drift  directions  to  date  do  not  exhibit  any  recog¬ 
nizable  patterns.  Measured  velocities  are  in  the 
50-150  m  s"*  range.  Figure  10,  however  provides 
evidence  that  for  the  magnetically  active  period  of 
February  9-12,  1986  the  diameter  of  the  convection 


158 


BUCHAU  ET  AL.:  POLAR  CAP  MEASUREMENTS  AND  REAL-TIME  MODELING 


531 


— F  HORIZONTAL  ^  . 

.  velocity  : 

■ '  _ _ '-I 

■  . • . ^ 


HORIZONTAL  ■ 
DIRECTION  ' 


.Veb'  oo'oo3f  ■  I  '•  'io7i 


lES  12b0  UT 


,00-  ..;;:;:::::;;..«oflizoNTM.  -«i 

««-  . . A  ■  ■  Ae-A'«-VA . r  - 

.  v|  1/  i  .......  .  ■  . . .  yy 

§;*'  .  :  "  r.’,  ,::  horizontal  direction  ;; 


§;*|  .  ;  horizontal  direction,. 

P.==r:Ep»^^ 

10  FEB  0000  UT  1  t1  FEB  liOO  UT 

•000  COLT 

4}0- . . .  . . .  . . . .  ,  .4$ 

,00- 1  A  ' ;  .  -thK-'-'-x  '  horizontal 

,50-  ■■■■;■  kJk  -/TlV  ’■  ••V,  VELOCITY 

‘••.^..^■‘.1: ,  V./J ..  V..  ...  ;L  fVERT.VEL. 

0*  . . .  .■..O...r.  ....  . .  01  .11..  —ft 


lY  vHfe  horizontal  direction  ^ 


KpN  4::  f— 4?  no  I  I  I  I  i  I  ^ 

'  11  FEB  0000  UT  I  12  FEB  1200  UT 

0000  COLT 

F  REGION  DRIFT 

DIGISONDE  OBSERVATIONS  AT  ARGENTIA,  NEWFOUNDLAND 
9/12  FEBRUARY  1986 

Fig.  10.  F  region  drill  observations  at  Argenlia,  Newfoundland  on  February  9-12,  1986  during  moderately 

disturbed  conditions. 


pattern  expands  to  bring  Argentia  into  the  area  of  Inspection  of  Figure  10  shows  evidence  for  west- 
sunward  return  flow.  For  this  station,  magnetic  mid-  ward  (sunward)  drift  at  75-150  m  s"*  for  Kp^4 
night  occurs  at  0400  UT.  For  each  of  the  3-hour  prior  to  local  magnetic  midnight,  with  occasional  ve- 
periods  marked  on  the  time  axes  of  the  three  panels,  locity  increases  to  450  m  s"‘.  During  less  active  con- 
the  Kp  values  are  indicated.  ditions  (2  <  Kp  <  4-)  the  data  for  February  10-11 
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show  switching  between  sunward  (westward)  and  an- 
tisunward  (eastward)  flow  prior  to  magnetic  mid¬ 
night.  This  suggests  a  sequence  of  convection  pattern 
contractions  and  expansions,  alternatingly  shifting 
the  convection  pattern  or  the  corotating  plasma  over 
the  station.  The  postmidnight  (CGLT)  data  on  Feb¬ 
ruary  10  and  11  indicate  generally  east  to  southeast 
(sunward)  drift  with  velocities  between  50  and  100  m 
s"'  for  quiet  conditions  (February  10)  and  between 
100  and  300  ms"’  for  more  active  conditions  (Feb¬ 
ruary  11).  Short  deviations  from  the  east  drift  toward 
northerly  were  observed  in  both  morning  periods. 

4.  DISCUSSION  OF  DRIFT  MEASUREMENTS 

A  newly  developed  digital  ionosonde  ground- 
based  technique  to  measure  plasma  drift  at  high  lati¬ 
tudes  has  produced  data  which  agree  with  previous 
observations  from  satellites  and  from  incoherent 
scatter  radars; 

1.  In  the  central  polar  cap  we  find  antisunward 
convection  for  B.  <  0  or  for  Kp  >  3. 

2.  During  quiet  conditions,  especially  when  fol¬ 
lowing  an  active  period,  sunward  convection  is  ob¬ 
served  in  the  polar  cap. 

3.  For  extended  quiet  periods  and  B.  >  0,  shear 
conditions  are  observed  typical  of  the  Sun-aligned 
polar  cap  arcs. 

4.  At  oval/suboval  latitudes  we  observe  sunward 
flow  (westward  before  midnight  and  eastward  after 
midnight)  under  magnetically  disturbed  conditions. 

5.  Argentia  (57®  CGL)  and  to  some  extent  also 
Goose  Bay  (65®  CGL)  show  random  switches  from 
westward  to  eastward  and  back  under  moderate 
{Kp  ^  2)  conditions.  This  suggests  that  the  station(s) 
move  alternatingly  under  the  influence  of  convecting 
or  corotating  plasmas. 

Important  parameters  of  the  plasma  convection, 
namely  the  antisunward  convection  velocity  and  the 
approximate  diameter  of  the  convection  pattern,  can 
be  measured  for  part  of  the  day. 

5.  DRIFT  MEASURE.MENTS  AND  MODELING 
OF  THE  HIGH-LATITUDE  IONOSPHERE 

The  importance  of  polar  convection  measurements 
lies  in  the  relation  of  the  convection  pattern  and  the 
convection  velocities  (magnitude  and  direction)  to 
the  condition  of  the  interplanetary  magnetic  field 
(IMF).  The  latter  controls  the  polar  cap  potential 
through  the  solar  wind  to  magnetosphere  energy 
coupling  function  c  [Akasofu,  1981,  1984;  Reiffet  at.. 


1981].  Ionospheric  modeling  has  shown  that  polar 
cap  and  auroral  oval  ionization  and  the  development 
of  the  F  layer  trough  to  the  south  of  the  oval  are 
controlled  (location  and  intensity)  by  the  magnitude 
of  the  polar  cap  potential  and  are,  therefore,  related 
to  the  polar  cap  convection  velocities.  We  expect  that 
future  studies  will  lead  to  the  development  of  either  c 
(IMF)  or  (IMF)  dependent  high-latitude  iono¬ 
spheric  models  which  can,  in  the  absence  of  measure¬ 
ments  for  either  parameter,  be  driven  by  polar  cap 
plasma  convection  measurements. 

As  an  example,  we  have  investigated  the  effects  of 
B.  on  the  electron  densities  in  the  central  polar  cap. 
Using  solar  production  for  the  time  dependent  mod¬ 
eling  of  the  ionization  buildup  (and  decay)  within  a 
flux  tube  (for  details  see  Anderson  et  al  [this  issue]) 
as  it  travels  along  a  trajectory  described  by  a  B, 
dependent  convection  model  [Hee/is,  1984],  the  time 
history  of  the  maximum  electron  density  (W„,j)  in 
flux  tubes  which  cross  Thule  at  selectd  times  has 
been  computed.  The  effects  due  to  particle  precipi¬ 
tation  have  been  neglected.  The  contribution  to 
by  the  cusp  precipitation  is  small  (~2.4  x  10“  cm"^), 
as  the  flux  tubes  that  typically  travel  at  500-1000  m 
s"'  spend  little  time  traversing  the  cusp  [Knudsen  et 
a/.,  1977]. 

Figure  1 1  shows  the  trajectories  of  the  flux  tubes 
which  traverse  Thule  at  1400  CGLT  (squares)  and 
2000  CGLT  (triangles)  in  a  CGL/CGLT  coordinate 
system  for  B.  <0  and  B.aO  condition,  respectively. 
Other  trajectories  arc  also  shown  but  are  not  dis¬ 
cussed  further.  The  locations  of  the  flux  tubes  are 
indicated  in  At  =  30  min  increments;  the  starting 
point  of  each  trajectory  is  6  hours  before  the  time  of 
crossing  over  Thule.  The  figure  illustrates  the  large 
difference  in  convection  pattern  diameter  and  in  the 
convection  velocities  (v  is  proportional  to  the  dis¬ 
tance  between  half-hour  marks)  for  the  two  selected 
B.  conditions.  Of  special  importance  is  the  exposure 
of  the  relevant  flux  tubes  to  solar  UV.  The  relation  of 
the  flux  tube  trajectories  to  the  terminator  for  De¬ 
cember  is  shown  in  Figure  12.  For  B.  <  0  the  flux 
lubes  arriving  at  Thule  al  1400  and  2000  CGLT 
spend  at  least  4.5  hours  in  sunlight,  much  of  that 
time  at  solar  zenith  angles  of  <80®.  Because  of  the 
high  convection  velocities,  it  takes  the  flux  tubes  only 
between  15  and  30  min  to  reach  Thule,  after  crossing 
the  terminator  into  darkness,  a  time  short  in  com¬ 
parison  to  the  recombination  time  of  the  F  layer 
ionization.  In  contrast,  for  B.  a  0  the  respective  flux 
tubes  experience  sunlight  for  only  1. 5-3.5  hours  and 
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CORRECTED  GEOMAGNETC  LOCAL  TIME  CORRECTED  GEOMAGNETIC  LOCAL  TIME 

CORRECTED  GEOMAGNETIC  LATITUDE  CORRECTED  GEO.MAGNET1C  LATITUDE 


SELECTED  CONVECTION  TRAJECTORIES 
TERMINATING  AT  THULE  (86“  CGL) 

At  =  30  min. 

Fig.  11.  Convection  trajectories  for  flux  tubes  crossing  Thule  at  0200  CGLT  (crosses),  0800  COLT  (circles). 
1400  CGLT  (squares)  and  2000  CGLT  (triangles)  for  B.  <  0  and  B,  a:  0  conditions,  in  CGL/CGLT  coordinates 
(Al  =  30  min).  The  2000  COLT  trajectory  has  not  been  shown  in  the  left  diagram,  for  clarity,  since  it  fcllow.s  closely 
the  1400  CGLT  trajectory. 


only  with  solar  zenith  angles  of  >  G5°.  As  a  conse¬ 
quence  of  the  low  convection  velocities  it  takes  the 
flux  tubes  between  45  min.  and  1.5  hours  to  move 
from  the  terminator  to  Thule.  Figure  13  shows  the 
results  of  the  modeling  of  the  ionization  time  histor¬ 
ies,  as  the  flux  tubes  travel  along  the  respective  tra¬ 
jectories.  The  results  differ  by  a  factor  of  3  (2)  in  Af„,„ 
for  ionization  observed  at  Thule  at  1400  (2000) 
CGLT,  for  the  two  considered  conditions,  D.  -  0 
and  B.  <  0,  respectively. 

6.  CONCLUSION 

A  chain  of  Digisondes,  with  automatic  plasma  drift 
measurement  capabilities,  is  able  to  continuously 
measure  features  of  the  polar  cap  convection  pat¬ 
terns,  which  previously  had  been  derived  from  many 
“snapshot”  measurements  by  polar  satellites,  and 


from  a  limited  (in  location  and  time)  set  of  incoherent 
scatter  radar  measuroments. 

Digisonde  drift  measurements  in  the  central  polar 
cap  identified  antisunward  convection  under  dis¬ 
turbed  conditions,  and  sunward  convection  when 
magnetic  conditions  switched  from  active  to  quiet 
Under  extended  quiet  conditions  the  existence  of 
shears  has  been  documented,  which  is  expected  to 
accompany  the  simultaneously  occurring  polar  cap 
auroral  E  and/or  F  layer  arcs.  The  observed  veloci¬ 
ties  in  the  polar  cap  tended  to  increase  with  magnetic 
activity,  corresponding  to  the  positive  correlation  be¬ 
tween  magnetic  activity  and  the  polar  cap  potential 
drop  \_Heppner  and  Maynafd.  19871[;  however,  a  clear 
overall  correlation  between  activity  and  velocity 
could  not  be  established.  This  may  possibly  be  due  to 
incorrect  urift  velocities  ootained  during  arc  or  quiet 
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^  L0C4L  sans  TIME  LOCAL  SOLAR  TIME 

GEOGRAPHICAL  LATITUDE  GEOGRAPHICAL  LATITUDE 


SELECTED  COUVECTION  TRAJECTORIES 
TERMINATING  AT  THULE  (86»  C6U 
At  *  30  min. 

Fig.  12,  The  siime  conveciion  trajectories  as  in  Figme  1 1,  here  in  geographic  latitude,  local  solar  time  coordinates. 

Also  ihowii  is  the  ierminmor(x  =  90°). 


conditions,  when  ihe  assuinplion  of  uniform  drift 
over  the  shy  map  area  docs  not  hold. 

At  the  auroral  and  .subauroral  stations  Goose  Bay 
and  Argentia  the  sunward  return  flow  is  observable 
for  several  hours  near  magnetic  local  midnight,  in- 
cluditig  the  expected  drift  reversal  from  westward  to 
eastward,  as  long  as  the  activity  is  high  enough  to 
bring  the  convection  patlc.’n  over  tnc  statiori(s). 
During  more  quiet  or  changing  magnetic  conditions 
the  switching  of  a  station  from  inside  the  convecting 
region  into  the  corotating  region  has  been  observed 
at  the  southernmost  station  (57°  CGL).  Thus  such 
drift  measurements  may  allow  oetcrminalion  of  the 
dimensions  of  the  region  under  the  influence  of  con¬ 
vection  for  at  least  part  of  the  day. 

Our  previously  discussed  measurements  and  other 
studies  of  polar  cap  convection  under  B.>  0  con¬ 
dition  iBurke  c!  al.,  1979;  Fran',  c-i  u!..  1986.  Heppner 
and  Maynard..  1987]  suggest  that  the  highly  struc¬ 
tured,  often  turbulent  convection,  represented  in 
multiple-cell  or  contorted  flow  patterns,  cannot  be 


deterniined  with  any  degree  of  accuracy  by  a  limited 
(even  large)  set  of  ionosonde  drift  measurements. 
However,  the  data  discussed  here,  backed  up  by  data 
from  a  full  year  of  Regular  World  Day  measurements 
(Qaanaaq,  1986),  provide  evidence  that  the  drift 
measurements  can  readily  identify  this  turbulent  con¬ 
dition. 

In  general,  the  slow  convection  under  B^>  0  con¬ 
dition  is  not  as  important  to  the  maintenance  or  en¬ 
hancement  of  the  polar  cap  (and  auroral  oval)  F 
layer  plasma  densities  as  the  higher  convection  veloc¬ 
ities  for  B.  <0  (increasing  magnetic  activity).  The 
previous  ionospheric  model  study  supports  this  as¬ 
sumption,  where  3  limes  higher  plasma  densities  over 
Thule  resulted  from  the  stronger  convection  for  the 
<  0,  in  comparison  to  those  for  the  B^  ~  0.  It  is 
expected  that  the  continous  drift  measurements  will 
define  the  nasic  parameters  of  the  convection  pattern 
at  least  for  periods  when  the  more  regular  two-cell 
configuration  prevails,  i.e.,  for  B,  £  0. 

But  even  for  B.  <  0  it  will  not  be  possible  to  speci- 
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Calculattd  NMAX  VakMS  Along  ConvacUon 
Tta)«cto!l«»  Crossing  Ihuls  at  14  and  20  MLT 

Fig.  13.  Calculated  values  along  convection  trajectories  crossing  Thule  at  1400  and  2000  CGLT. 


fy  the  details  of  the  convection  pattern  for  tl"  “’'hole 
polar  region  with  the  limited  number  of  polar  sta¬ 
tions  currently  envisioned  (five  to  10).  We  do  expect, 
however,  on  the  basis  of  the  data  presented  in  this 
paper,  that  a  detailed  analysis  of  the  existing  and 
growing  drift  data  base  will  lead  to  a  technique  of 
selecting  (for  at  least  part  of  the  day)  the  convection 
model  and  the  polar  cap  potential  drop  which  best 
represent  the  observations.  The  improving  convec¬ 
tion  models,  the  addition  of  Digisonde  drift  measure¬ 
ments  at  locations  significant  to  the  determination  of 
the  convection  pattern,  and  the  combined  use  of  ion 
drift  data  from  satellites  (Defense  Meteorological 
Satellite  Program)  and  time  continuous  Digisonde 
drift  measurements  hold  the  promise  for  a  convection 
pattern  specification  capability.  A  first  step  in  this 
direction  is  the  planned  deployment  of  a  Digisonde 
at  Sondrestrom  to  measure  the  skewing  of  the  trans¬ 
cusp  plasma  flow  under  the  influence  of  By  [_Heelis, 
1984]  with  the  goal  of  providing  in  the  skewing  an¬ 
other  parameter  for  the  selection  of  the  best  fitting 
model. 
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